D-A189  715  MODEL  SELECTION  FOR  THE  MULTIPLE  MODEL  AOAPTIVE 

ALGORITHM  FOR  IN-FLIGHT  SIMULATIONS)  AIR  FORCE  INST  OF 
TECH  HRIGHT-PATTERSON  AFB  OH  J  R  MATHES  DEC  87 
UNCLASSIFIED  AFIT/GE/ENG/870-48  F/G  1/1 


H 


DEPARTMENT  OF  THE  AIR  FORCE 
AIR  UNIVERSITY 

AIR  FORCE  INSTITUTE  OF  TECHNOLOGY 


Wright-Patterson  Air  Force  Base,  Ohio 
DISTRIBUTION  STATEMENT  A. 

Approval  for  public  ralocms; 


AFIT/GE/ENG/87D-40 


MODEL  SELECTION  FOR  THE  MULTIPLE 
MODEL  ADAPTIVE  ALGORITHM  FOR 
IN-FLIGHT  SIMULATION 

THESIS 


James  R.  Matthes,  Jr. 
Captain,  USAF 

AFIT/GE/ENG/87D-40 


DTIC 


S^LECTEI 
MAR  0  7  1988 

c^H 


Approved  for  public  release;  distribution  unlimited 


AFIT/GE/ENG/87D-40 


MODEL  SELECTION  FOR  THE  MULTIPLE 
MODEL  ADAPTIVE  ALGORITHM  FOR 
IN-FLIGHT  SIMULATION 

THESIS 


Presented  to  the  Faculty  of  the  School  of  Engineering 
of  the  Air  Force  Institute  of  Technology 
Air  University 

In  Partial  Fulfillment  of  the 
Requirements  for  the  Degree  of 
Master  of  Science  in  Electrical  Engineering 


James  R.  Matthes,  Jr.,  B.S. 
Captain,  USAF 

December  1987 


Approved  for  public  release;  distribution  unlimited 


e 


Acknowledgments 


The  thesis  represents  the  culmination  of  my  AFIT 
experience.  This  effort  was  accomplished  to  a  large  degree 
by  depriving  myself  of  one  of  my  favorite  recreational 
pursuits  --  sleep.  Being  cognizant  of  my  limited  ability,  I 
realize  that  without  the  guidance  and  suggestions  of  my 
peers  and  professors  this  effort  would  have  not  been 
possible . 

My  special  thanks  are  extended  to  the  AFIT  faculty 
members  who  patiently  and  expertly  guided  me  through  this 
research  effort.  The  technical  knowledge,  suggestions,  and 
dedication  of  my  thesis  advisor  Professor  John  J.  D’Azzo 
were  instrumental  in  the  success  of  this  thesis.  I  also 
wish  to  express  thanks  to  Professor  Peter  S.  Maybeck  for 
whom  I  have  the  utmost  respect.  His  tireless  reviews  of  the 
thesis  drafts  provided  many  constructive  technical  improve¬ 
ments  to  my  thesis  and  enhanced  my  understanding  of  the 


*<N 


*  V 

O,  O 


material  presented  in  this  document. 

Gratitude  is  also  due  to  Captain  Julio  Velez  and 
Lieutenant  Tom  Berens.  The  many  stimulating  and 


•jslon  For 

GRAJtl 

TAB 

Junced 

'lcatlon_ 


enlightening  discussions  we  had  were  enjoyable  and  were  an 


invaluable  addition  to  this  thesis 


'Distribution/ 

j _ Availability  Codes 

■  Avail  and/or 

Dlst  Special 

Vx  ! 


□  □ 


c 


Heartfelt  thanks  are  also  due  to  my  wonderful  wife, 
Helene,  and  my  two  precious  children,  Noelle  and  Aaron. 
Helene  constantly  provided  me  with  the  love,  encouragement, 
and  understanding  that  seemed  to  ease  the  pain  associated 
with  the  AFIT  experience  while  Noelle  and  Aaron  were  always 
there  to  keep  my  focus  on  the  important  matters  of  life. 

Finally,  I  wish  to  express  my  love  and  appreciation  to 
the  Lord  Jesus  Christ  for  the  strength  and  guidance  that  He 
is  continually  willing  to  provide.  I  realize  that  "without 
Him  I  can  do  nothing..."  John  15:5. 

--  Captain  James  R.  Matthes,  Jr. 


Q 

G 

9 

9 


i  i  i 


wm 

C 

Table  of  Contents 

• 

Page 

Acknowledgments  . . 

ii 

• 

List 

of  Figures  . 

vi  i 

List 

of  Tables  . 

xxi  i 

Abstract  . 

xxiv 

i . 

• 

1  . 

Introduction  . 

1 

» 

1  .  1 

Overview  . 

1 

* 

1.2 

Background  . 

2 

f 

’* 

1 . 3 

Problem  Statement  . 

3 

r 

4 

e 

1 . 4 

Summary  of  Current  Knowledge  . . 

4 

1.4.1  In-Flight  Simulation  . 

4 

.V 

1.4.2  Adaptive  Control  . 

6 

1.4.3  MMAE  . 

11 

1  .  5 

Assumptions  . . . 

12 

• 

1.6 

Approach  . . . 

13 

1.7 

Overview  . . . 

15 

2. 

Aircraft  Description  and  Models  . 

16 

• 

2  .  1 

AFTI/F-16  Description  . 

16 

2.2 

Aircraft  Models  . 

20 

2.2.1  State  Space  Model  . 

2.2.2  Autoregressive  Difference 

20 

* 

Equation  Model  . 

26 

j 

c 

2.3 

Actuator  Model  . 

28 

•«< 

2.4 

Summary  . . . 

29 

,» 

f/, 

■: 

3. 

Control  Law  Algorithm  . 

30 

5 

c 

3.1 

Introduction  . 

30 

*• 

3.2 

Model  Following  Techniques  . 

30 

i* 

3.3 

Porter’s  Control  Law  . 

34 

i* 

3.3.1  Fixed  Gain  Controller  Matrices  . 

37 

i* 

c 

3.3.2  Adaptive  Controller  Matrices  . 

43 

m 

)' 

3.4 

Summary  . 

45 

t’ 

tr 

»' 

» 

» 

t 

t 

0 

iv 

‘  *  „  *  .  » 

r.v.  r 

4.  Parameter  Estimation  .  47 

4.1  Introduction  .  47 

4.2  Multiple  Model  Algorithm  .  48 

4.2.1  Prediction  Error  .  53 

4.2.2  Prediction  Error  Covariance  .  54 

4.3  Multiple  Model  Parameter  Estimator  .  55 

4.4  Summary  .  59 

5.  Model  Selection  Criteria  .  60 

5.1  Introduction  .  60 

5.2  Performance  Boundaries  .  61 

5.3  Sensor  Noise  Effects  .  67 

5.4  Summary  .  71 

6.  Simulation  Setup  and  Results  .  73 

6.1  Simulation  Setup  .  73 

6.2  Simulation  Results  .  77 

6.2.1  Single  Model  Analysis  .  77 

6.2. 1.1  Performance  Boundary- 

Evaluation  . 77 

6.2. 1.2  Control  Surface  Activity  ...  101 

6.2.2  Two-Model  Configurations  .  106 

6.2.3  Three-Model  Configurations  .  123 

6.2.4  Control  Law  Gain  Adjustment  .  150 

6.2.5  Sensor  Noise  Effects  .  175 

6.3  Summary  . 215 

7.  Conclusions  and  Recommendations  .  217 

7.1  Introduction  . . 217 

7.2  Conclusions  .  217 

7.3  Recommendations  for  Further  Study  .  221 

Appendix  A  .  223 

Appendix  B  .  242 

Appendix  C  . 247 

Appendix  D  .  251 


v 


Bibliography  . . . . .  260 

Vita  .  263 


o 


a 


vi 


6 


& 


6 


List  of  Figures 


Figure  Page 

1- 1  Parameter  Adaptive  Control  System  .  9 

2- 1  The  AFTI/F-16  .  17 

2-2  AFTI/F-16  Control  Surfaces  and  Definition  of 

Positive  Deflections  .  19 

2-3  Positive  Directions  of  Forces,  Moments,  and 

Angles  in  the  Body  Reference  Frame  .  22 

2- 4  Flight  Envelope  of  Interest  .  23 

3- 1  Explicit  Model  Following  System  .  32 

3-2  Improved  Explicit  Model  Following  System  .  33 

3- 3  System  Block  Diagram  -  Discrete  Case  .  38 

4- 1  Multiple  Model  Parameter  Estimator  .  56 

5- 1  Region  of  Performance  Robustness  .  63 

5-2  Intersection  of  Performance  Boundaries  .  67 

5- 3  Noise  Effects  on  Model  Selection  .  72 

6- 1  Top  Level  Description  of  Parameter 

Adaptive  System  Simulation  Package  .  74 

6-2  Flight  Path  Angle  Reference  Tracking  Signal  ...  76 

6-3  Pitch  Rate  Reference  Tracking  Signal  .  76 

6-4  Pitch  Rate  Response 

Nominal  Model:  18K  0.45M 

Operating  Point:  18K  0.45M  .  79 


vi  i 


6 


Pitch  Rate  Performance  Criterion 
Nominal  Model:  1_K  0.45M 
Operating  Point:  18K  0.45M  . 


6-6  Flaperon  Position 

Nominal  Model:  18K  0.45M 
Operating  Point:  18K  0.45M 

6-7  Elevator  Position 

Nominal  Model:  18K  0.45M 
Operating  Point:  18K  0.45M 

6-8  Flaperon  Rate 

Nominal  Model:  18K  0.45M 
Operating  Point:  18K  0.45M 

6-9  Elevator  Rate 

Nominal  Model:  18K  0.45M 
Operating  Point:  18K  0.45M 

6-10  Pitch  Rate  Response 

Nominal  Model:  18K  0.45M 
Operating  Point:  18K  0.8M  . 


6-11  Pitch  Rate  Performance  Criterion 
Nominal  Model:  18K  0.45M 
Operating  Point:  18K  0.8M  . 


6-12  Flaperon  Position 

Nominal  Model:  18K  0.45M 
Operating  Point:  18K  0.8M 

6-13  Elevator  Position  ’ 

Nominal  Model:  18K  0.45M 
Operating  Point:  18K  0 . 8M 

6-14  Flaperon  Rate 

Nominal  Model:  18K  0.45M 
Operating  Point:  18K  0 . 8M 

6-15  Elevator  Rate 

Nominal  Model:  18K  0.45M 
Operating  Point:  18K  0 . 8M 


VI  1  1 


V  . 


7T77TF —  '*»,V4V  V  V  V,  V 


y: v:y: y: ^ 'y v 'v  y y.  v v yyyyyy 


# 


6-16  Performance  Boundary  for  Nominal  Model 

at  18K  0.45M  .  88 

6-17  Performance  Boundary  for  Nominal  Model 

at  38K  0.7M  .  89 

6-18  Pitch  Rate  Response 

Nominal  Model:  18K  0.45M 

Operating  Point:  26K  0.45M  .  93 

6-19  Pitch  Rate  Performance  Criterion 
Nominal  Model:  18K  0.45M 

Operating  Point:  26K  0.45M  .  93 

6-20  Flaperon  Position 

Nominal  Model:  18K  0.45M 

Operating  Point:  26K  0.8M  .  94 

6-21  Elevator  Position 

Nominal  Model:  18K  0.45M 

Operating  Point:  26K  0.45M  .  94 

6-22  Flaperon  Rate 

Nominal  Model:  18K  0.45M 

Operating  Point:  26K  0.45M  .  95 

6-23  Elevator  Rate 

Nominal  Model:  18K  0.45M 

Operating  Point:  26K  0.45M  .  95 

6-24  Pitch  Rate  Response 

Nominal  Model:  38K  0.7M 

Operating  Point:  26K  0.45M  .  96 

6-25  Pitch  Rate  Performance  Criterion 
Nominal  Model:  38K  0.7M 

Operating  Point:  26K  0.45M  .  96 


6-26  Flaperon  Position 

Nominal  Model:  38K  0 . 7M 
Operating  Point:  26K  0.45M 


97 


6 

6-27 

Elevator  Position 

Nominal  Model:  38K  0.7M 

# 

Operating  Point:  26K  0.45M  . 

97 

6-28 

Flaperon  Rate 

Nominal  Model:  38K  0.7M 

Operating  Point:  26K  0.45M  . 

98 

* 

6-29 

Elevator  Rate 

Nominal  Model:  38K  0.7M 

Operating  Point:  26K  0.45M  . 

98 

• 

6-30 

Lift  Curve  . 

101 

6-31 

Performance  Boundaries  for  Configuration  1  .... 

108 

6-32 

Performance  Boundaries  for  Configuration  2  .... 

109 

C 

6-33 

Performance  Boundaries  for  Configuration  3  .... 

110 

6-34 

Performance  Boundaries  for  Configuration  4  .... 

111 

• 

6-35 

Performance  Boundaries  for  Configuration  5  .... 

112 

6-36 

Pitch  Rate  Response/Configuration  2 

Operating  Point:  14K  0.35M  . 

114 

6-37 

Pitch  Rate  Performnace  Criterion 

o 

Configuration  2 

Operating  Point:  14K  0.35M  . 

114 

1 

i 

i 

6-38 

Model  Probability  Weightings 

\ 

f 

Configuration  2 

Operating  Point:  14K  0.35M  . 

115 

6-39 

Pitch  Rate  Response/Configuration  4 

Operating  Point:  14K  0 . 9M  . 

116 

c 

1 

f 

1 

» 

t 

1 

6-40 

Pitch  Rate  Performance  Criterion 

Configuration  4 

Operating  Point:  14K  0 . 9M  . 

116 

. 

'c 

X 

\ 

L . 

V  Vi 

W.  W  W  W  '"  ' 
y.v,  v,v4\*  v  v 

€ 


m 


e 


f 


c 


* 


6-41  Model  Probability  Weightings 
Configuration  4 

Operating  Point:  14K  0 . 9M  .  117 

6-42  Pitch  Rate  Response/Configuration  3 

Operating  Point:  10K  0.9M  .  118 

6-43  Pitch  Rate  Performance  Criterion 
Configuration  3 

Operating  Point:  10K  0.9M  . . .  118 

6-44  Model  Probability  Weightings 
Configuration  3 

Operating  Point:  10K  0 . 9M  .  119 

6-45  Performance  Boundary  Overlap  .  122 

6-46  Three  Model  Configurations  .  125 

6-47  Performance  Boundaries  for  Configuration  6  ....  126 

6-48  Performance  Boundaries  for  Configuration  7  ....  127 


6-49  Performance  Boundaries  for  Configuration  8  ....  128 


6-50  Model  Selection  Data  for  Configuration  6  .  142 

6-51  Model  Selection  Data  for  Configuration  7  .  143 

6-52  Model  Selection  Data  for  Configuration  8  .  144 

6-53  Model  Probability  Weightings 

Configuration  6 

Operating  Point:  26K  0.45M  . .  145 

6-54  Model  Probability  Weightings 
Configuration  8 

Operating  Point:  26K  0.65M  .  147 

6-55  Pitch  Rate  Response 

Nominal  Model:  22K  0.65M 
Operating  Point:  22K  0.65M 
SIG1 =0 . 4  SIG2=0 . 7  RHO=0.8 


153 


6-56  Pitch  Rate  Performance  Criterion 
Nominal  Model:  22K  0.65M 
Operating  Point:  22K  0.65M 

SIG1 =0 . 4  SIG2=0 . 7  RHO=0.8  .  153 

6-57  Flaperon  Position 

Nominal  Model:  22K  0.65M 
Operating  Point:  22K  0.65M 

SIG1 =0 . 4  SIG2=0 . 7  RHO=0.8  .  154 

6-58  Elevator  Position 

Nominal  Model:  22K  0.65M 
Operating  Point:  22K  0.65M 

S IG1 =0 . 4  SIG2=0 . 7  RHO=0.8  .  154 

6-59  Flaperon  Rate 

Nominal  Model:  22K  0.65M 
Operating  Point:  22K  0.65M 

SIG1 =0 . 4  SIG2=C . 7  RHO=0.8  .  155 

6-60  Elevator  Rate 

Nominal  Model:  22K  0.65M 
Operating  Point:  22K  0.65M 

SIG1 =0 . 4  SIG2=0.7  RHO=0.8  .  155 

6-61  Pitch  Rate  Response 

Nominal  Model:  22K  0.65M 
Operating  Point:  22K  0.65M 

SIG1 =0 . 4  S IG2=0 . 7  RHO=0.4  .  158 

6-62  Pitch  Rate  Performance  Criterion 
Nominal  Model:  22K  0.65M 
Operating  Point:  22K  0.65M 

SIG1 =0 . 4  SIG2=0 . 7  RHO=0.4  .  158 

6-63  Flaperon  Position 

Nominal  Model:  22K  0.65M 
Operating  Point:  22K  0.65M 

SIG1 =0 . 4  SIG2=0 . 7  RHO=0.4  .  159 

6-64  Elevator  Position 

Nominal  Model:  22K  0.65M 
Operating  Point:  22K  0.65M 
SIG1 =0 . 4  S IG2  =  0 . 7  RHO=0.4 


c 


159 


c 


c 


6-65  Flaperon  Rate 

Nominal  Model:  22K  0.65M 
Operating  Point:  22K  0.65M 
SIG1=0 . 4  S IG2=0 . 7  RHO=0.4  . 

6-66  Elevator  Rate 

Nominal  Model:  22K  0.65M 
Operating  Point:  22K  0.65M 
SIG1 =0 . 4  S IG2=0 . 7  RHO=0.4  . 

6-67  Pitch  Rate  Response 

Nominal  Model:  22K  0.65M 
Operating  Point:  22K  0.65M 
SIG1 =0 . 8  SIG2=0 . 7  RHOsO.8  . 

6-68  Pitch  Rate  Performance  Criterion 
Nominal  Model:  22K  0.65M 
Operating  Point:  22K  0.65M 
SIG1 =0 . 8  SIG2=0 . 7  RHO=0.8  . 

6-69  Flaperon  Position 

Nominal  Model:  22K  0.65M 
Operating  Point:  22K  0.65M 
SIG1 =0 . 8  SIG2=0 . 7  RHO=0.8  . 

6-70  Elevator  Position 

Nominal  Model:  22K  0.65M 
Operating  Point:  22K  0.65M 
SIG1 =0 . 8  SIG2=0 . 7  RHO=0.8  . 

6-71  Flaperon  Rate 

Nominal  Model:  22K  0.65M 
Operating  Point:  22K  0.65M 
SIG1=0 . 8  SIG2=0 . 7  RHO=0.8  . 

6-72  Elevator  Rate 

Nominal  Model:  22K  0.65M 
Operating  Point:  22K  0.65M 
SIG1 =0 . 8  SIG2=0 . 7  RHO=0.8  . 


160 


160 


161 


161 


162 


162 


163 


163 


e 


xi  i  i 


c 


c 


Q 


6-73  Pitch  Rate  Response 

Nominal  Model:  22K  0.65M 
Operating  Point:  22K  0.65M 

SIG1 =0.1  SIG2=0 . 7  RHO=0 . 8  .  164 

6-74  Pitch  Rate  Performance  Criterion 
Nominal  Model:  22K  0.65M 
Operating  Point:  22K  0.65M 

SIG1 =0 . 1  S IG2=0 . 7  RHO=0 . 8  .  164 

6-75  Flaperon  Position 

Nominal  Model:  22K  0.65M 
Operating  Point:  22K  0.65M 

SIGU0.1  SIG2  =  0 . 7  RHO=0 . 8  .  165 

6-76  Elevator  Position 

Nominal  Model:  22K  0.65M 
Operating  Point:  22K  0.65M 

SIG1 =0 . 1  S IG2=0 . 7  RHO=0 . 8  .  165 

6-77  Flaperon  Rate 

Nominal  Model:  22K  0.65M 
Operating  Point:  22K  0.65M 

SIG1 =0.1  SIG2=0 . 7  RHOsO.8  .  166 

6-78  Elevator  Rate 

Nominal  Model:  22K  0.65M 
Operating  Point:  22K  0.65M 

SIG1 =0 . 1  SIG2=0 . 7  RHO=0 . 8  .  166 

6-79  Pitch  Rate  Response 

Nominal  Model:  22K  0.65M 
Operating  Point:  22K  0.65M 

SIG1 =0 . 4  SIG2= 1 . 5  RHO=0 . 8  .  168 

6-80  Pitch  Rate  Performance  Criterion 
Nominal  Model:  22K  0.65M 
Operating  Point:  22K  0.65M 

S IG1 =0 . 4  SIG2= 1 . 5  RHO=0 . 8  .  168 

6-81  Flaperon  Position 

Nominal  Model:  22K  0.65M 
Operating  Point:  22K  0.65M 

SIG1 =0 . 4  SIG2=1.5  RHO=0 . 8  .  169 


6-82  Elevator  Position 

Nominal  Model:  22K  0.65M 
Operating  Point:  22K  0.65M 

SIG1=0 . 4  SIG2= 1 . 5  RHO=0 . 8  .  169 

6-83  Flaperon  Rate 

Nominal  Model:  22K  0.65M 
Operating  Point:  22K  0.65M 

SIG1 =0 . 4  SIG2=1.5  RHO=0 . 8  .  170 

6-84  Elevator  Rate 

Nominal  Model:  22K  0.65M 
Operating  Point:  22K  0.65M 

SIG1=0 . 4  SIG2= 1 . 5  RHO=0 . 8  .  170 

6-85  Pitch  Rate  Response 

Nominal  Model:  22K  0.65M 
Operating  Point:  22K  0.65M 

SIG1=0.4  SIG2=0 . 1  RHO=0 . 8  .  171 

6-86  Pitch  Rate  Performance  Criterion 
Nominal  Model:  22K  0.65M 
Operating  Point:  22K  0.65M 

SIG1 =0 . 4  S IG2=0 . 1  RHO=0.8  .  171 

6-87  Flaperon  Position 

Nominal  Model:  22K  0.65M 
Operating  Point:  22K  0.65M 

SIG1 =0 . 4  SIG2=0 . 1  RHO=0 . 8  .  172 

6-88  Elevator  Position 

Nominal  Model:  22K  0.65M 
Operating  Point:  22K  0.65M 

SIG1 =0 . 4  SIG2=0 . 1  RHO=0 . 8  .  172 

6-89  Flaperon  Rate 

Nominal  Model:  22K  0.65M 
Operating  Point:  22K  0.65M 

SIG1 =0 . 4  SIG2=0 . 1  RHO=0 . 8  .  173 


6-90  Elevator  Rate 

Nominal  Model:  22K  0.65M 
Operating  Point:  22K  0.65M 
SIG1 =0 . 4  SIG2=0 . 1  RHO=0 . 8  . 

6-91  Pitch  Rate  Response 
Configuration  8 
Operating  Point:  10K  0 . 9M 
SIG1 =0 . 4  SIG2=0.1  RHO=0 . 8  . 

6-92  Pitch  Rate  Performance  Criterion 
Configuration  8 
Operating  Point:  10K  0 . 9M 
SIG1 =0 . 4  SIG2=0 . 1  RHO  =  0 . 8  . 


6-93  Flaperon  Position 
Configuration  8 
Operating  Point:  10K  0.9M 
SIG1 =0 . 4  SIG2=0 . 1  RHO=0 . 8 


6-94  Elevator  Position 
Configuration  8 
Operating  Point:  10K  0 . 9M 
SIG1 =0 . 4  SIG2=0 . 1  RHOsO . 8 


6-95  Flaperon  Rate 

Configuration  8 
Operating  Point:  10K  0 . 9M 
SIG1 =0 . 4  SIG2=0.1  RHO=0 . 8 


6-96  Elevator  Rate 

Configuration  8 
Operating  Point:  10K  0 . 9M 
SIG1 =0 . 4  SIG2=0.1  RHO=0 . 8 


6-97  Pitch  Rate  Response 
Configuration  8 
Operating  Point:  38K  0 . 9M 
SIG 1 =0 . 4  SIG2  =  0.1  RHO  =  0 . 8 


6-98  Pitch  Rate  Performance  Criterion 
Configuration  8 
Operating  Point:  38K  0 . 9M 
SIG1 =0 . 4  SIG2=0.1  RHO=0 . 8  . 


9 


e 


• 

6-99 

Flaperon  Position 

Configuration  8 

Operating  Point:  38K  0 . 9M 

SIG1 =0 . 4  SIG2=0.1  RHO=0 . 8  . 

•  «  •  t 

.  180 

• 

6-100 

Elevator  Position 

Configuration  8 

Operating  Point:  38K  0.9M 

SIG1 =0 . 4  SIG2=0 . 1  RHO=0 . 8  . 

•  •  •  •  . 

.  180 

• 

6-101 

Flaperon  Rate 

Configuration  8 

Operating  Point:  38K  0 . 9M 

SIG1 =0 . 4  S IG2=0 . 1  RHO=0 . 8  . 

•  •  •  •  > 

.  181 

• 

6-102 

Elevator  Rate 

Configuration  8 

Operating  Point:  38K  0.9M 

SIG1 =0 . 4  SIG2=0 . 1  RHO=0 . 8  . 

•  •  •  •  i 

.  181 

• 

6-103 

Model  Probability  Weightings 
Operating  Point:  10K  0.35M  .... 

•  •  •  •  i 

.  183 

6-104 

Prediction  Error  Variance  (1,1) 
Operating  Point:  10K  0.35M  .... 

for 

•  •  •  •  < 

Model 

1 

.  184 

• 

6-105 

Prediction  Error  Variance  (2,2) 
Operating  Point:  10K  0.35M  .... 

f  or 

•  •  •  •  i 

Model 

1 

.  184 

6-106 

Prediction  Error  Variance  (1,1) 
Operating  Point:  10K  0.35M  .... 

for 

t  •  •  •  i 

Model 

2 

.  185 

£ 

6-107 

Prediction  Error  Variance  (2,2) 
Operating  Point:  10K  0.35M  .... 

f  or 

>  •  •  •  i 

Model 

2 

»  •  •  •  • 

.  185 

£ 

6-108 

Model  Probability  Weightings 
Operating  Point:  10K  0.35M 
Sensor  Noise  of  Table  6-15  ..... 

»  •  •  •  i 

»  *  •  •  « 

.  187 

6-109 

Prediction  Error  Variance  (1,1) 
Operating  Point:  10K  0.35M 
Sensor  Noise  of  Table  6-15  .... 

for 

•  •  •  •  i 

Model 

1 

»  •  •  •  • 

.  188 

0 


6-110  Prediction  Error  Variance  (2,2)  for  Model  1 
Operating  Point:  10K  0.35M 

Sensor  Noise  of  Table  6-15  .  188 

6-111  Prediction  Error  Variance  (1,1)  for  Model  2 
Operating  Point:  10K  0.35M 

Sensor  Noise  of  Table  6-15  .  189 

6-112  Prediction  Error  Variance  (2,2)  for  Model  2 
Operating  Point:  10K  0.35M 

Sensor  Noise  of  Table  6-15  .  189 

6-113  Model  Probability  Weightings 
Operating  Point:  10K  0.35M 

Sensor  Noise  *  100  . .  190 

6-114  Prediction  Error  Variance  (1,1)  for  Model  1 
Operating  Point:  10K  0.35M 

Sensor  Noise  *  100  .  191 

6-115  Prediction  Error  Variance  (2,2)  for  Model  1 
Operating  Point:  10K  0.35M 

Sensor  Noise  *  100  .  191 

6-116  Prediction  Error  Variance  (1,1)  for  Model  2 
Operating  Point:  10E  0.35M 

Sensor  Noise  *  100  .  192 

6-117  Prediction  Error  Variance  (2,2)  for  Model  2 
Operating  Point:  10K  0.35M 

Sensor  Noise  *  100  .  192 

6-118  Model  Probability  Weightings 

Operating  Point:  10K  0.35M  .  195 

6-119  Prediction  Error  Variance  (1,1)  for  Model  1 

Operating  Point:  10K  0.35M  .  196 

6-120  Prediction  Error  Variance  (2,2)  for  Model  1 

Operating  Point:  10K  0.35M  . .  196 

6-121  Prediction  Error  Variance  (1,1)  for  Model  2 

Operating  Point:  10K  0.35M  .  197 


xv  i  i  i 


e 


6-122 

Prediction  Error  Variance  (2,2) 
Operating  Point:  10K  0.35M  .... 

for 

*  •  •  •  i 

Model 

2 

i  •  •  •  • 

.  197 

# 

6-123 

Model  Probability  Weightings 
Operating  Point:  10K  0.35M 
Sensor  Noise  *  100  . 

.  198 

• 

6-124 

Prediction  Error  Variance  (1,1) 
Operating  Point:  10K  0.35M 
Sensor  Noise  *  100  . 

for 

•  •  •  •  i 

Model 

1 

.  199 

6-125 

Prediction  Error  Variance  (2,2) 
Operating  Point:  10K  0.35M 
Sensor  Noise  *  100  . . 

f  or 

Model 

1 

>  •  •  •  • 

.  199 

• 

6-126 

Prediction  Error  Variance  (1,1) 
Operating  Point:  10K  0.35M 
Sensor  Noise  *  100  . 

for 

•  •  •  •  i 

Model 

2 

i  «  •  •  • 

.  200 

1 

s 

* 

I 

6-127 

Prediction  Error  Variance  (2,2) 
Operating  Point:  10K  0.35M 
Sensor  Noise  *  100  . 

for 

Model 

2 

.  200 

; 

• 

6-128 

Model  Probability  Weightings 
Configuration  8 

Operating  Point:  26K  0 . 9M  . . 

•  •  •  •  \ 

i  •  •  •  • 

.  203 

• 

6-129 

Prediction  Error  Variance  (1,1) 
Configuration  8 

Operating  Point:  26K  0.9M  . . 

for 

>  •  •  •  i 

Model 

1 

i  •  •  •  • 

.  204 

» 

> 

1, 

1 

! 

6-130 

Prediction  Error  Variance  (2,2) 
Configuration  8 

Operating  Point:  26K  0.9M  . 

for 

Model 

1 

i  t  *  t  * 

.  204 

i 

i 

> 

i 

‘  r- 

\  ' 

6-131 

Prediction  Error  Variance  (1,1) 
Configuration  8 

Operating  Point:  26K  0.9M  . 

for 

Model 

2 

.  205 

\ 

6-132 

Prediction  Error  Variance  (2,2) 
Configuration  8 

Operating  Point:  26K  0 . 9M  . 

f  or 

•  •  «  •  i 

Model 

2 

.  205 

e 


m 


« 


i 

» 


4 

4 

4 

i 


t 


6-133  Prediction  Error  Variance  (1,1)  for  Model  3 
Configuration  8 

Operating  Point:  26K  0 . 9M  .  206 

6-134  Prediction  Error  Variance  (2,2)  for  Model  3 
Configuration  8 

Operating  Point:  26K  0.9M  .  206 


6-135  Model  Probability  Weightings 
Configuration  8 
Operating  Point:  26K  0 . 9M 

Sensor  Noise  of  Table  6-15  .  207 

6-136  Model  Probability  Weightings 
Configuration  8 
Operating  Point:  26K  0.9M 

Sensor  Noise  *  100  . . .  208 

6-137  Prediction  Error  Variance  (1,1)  for  Model  1 
Configuration  8 
Operating  Point:  26K  0.9M 

Sensor  Noise  *100  . .  209 

6-138  Prediction  Error  Variance  (2,2),  for  Model  1 
Configuration  8 
Operating  Point:  26K  0.9M 

Sensor  Noise  *  100  .  209 

6-139  Prediction  Error  Variance  (1,1)  for  Model  2 
Configuration  8 
Operating  Point:  26K  0.9M 

Sensor  Noise  *  100  .  210 

6-140  Prediction  Error  Variance  (2,2)  for  Model  2 
Configuration  8 
Operating  Point:  26K  0 . 9M 

Sensor  Noise  *  100  .  210 

6-141  Prediction  Error  Variance  (1,1)  for  Model  3 
Configuration  8 
Operating  Point:  26K  0.9M 

Sensor  Noise  *  100  .  211 


xx 


C 


6-142  Prediction  Error  Variance  (2,2)  for  Model  3 
Configuration  8 
Operating  Point:  26K  0 . 9M 

Sensor  Noise  *  100  . 


6-143  Model  Probability  Weightings 
Configuration  8 

Operating  Point:  38K  0.6M  . .  213 

6-144  Model  Probability  Weightings 
Configuration  8 
Operating  Point:  38K  0 . 6M 

Sensor  Noise  of  Table  6-15  .  213 

6-145  Model  Probability  Weightings 
Configuration  8 
Operating  Point:  38K  0 . 6M 

Sensor  Noise  *  10  . . .  214 


6-146  Model  Probability  Weightings 
Configuration  8 
Operating  Point:  38K  0 . 6M 
Sensor  Noise  *  100  . 


List  of  Tables 


Table  Page 

2-1  Flight  Conditions  of  Interest  .  24 

2-2  Control  Surface  Position  and  Rate  Limits  .  29 

6-1  Flight  Conditions  Failing  Performance  Criteria 

Nominal  Flight  Condition=18 , 000  feet  0.45  Mach  .  86 

6-2  Flight  Conditions  Failing  Performance  Criteria 

Nominal  Flight  Condition=38 , 000  feet  0.7  Mach  ..  87 

6-3  Flight  Conditions  of  Nominal  Models  with 

Similar  Dynamic  Pressures  .  90 

6-4  Eigenvalues  for  the  A  matrix  of  State 

Space  Model  . 99 

6-5  Angle  of  Attack  for  Nominal  Flight  Conditions  . .  100 

6-6  Aerodynamic  Gains  for  Transfer  Functions  of 

Equations  (6-3)  through  (6-6)  104 

6-7  Aerodynamic  Gains  for  Transfer  Functions  of 

Equations  (6-8)  through  (6-11)  105 

6-8  Two  Model  Configurations  .  107 

6-9  Three  Model  Configurations  . 124 

6-10  Model  Selection  Data  for  Configuration  6  .  130 

6-11  Model  Selection  Data  for  Configuration  7  .  134 

6-12  Model  Selection  Data  for  Configuration  8  .  138 

6-13  Performance  Index  Comparison  of  Selected 

Flight  Conditions  .  149 

6-14  Design  Parameter  Values  for  Simulations 

at  22,000  ft  0.65  Mach  . . . .  156 


15  Sensor  Noise  Data  . 

16  Two  Model  Configuration  for  Noise 

Consideration  . 

Two  Model  Configuration  for  Noise 
Consideration  . 


xx  i  i  i 


AFIT/GE/ENG/87D-40 


ABSTRACT 

'This  thesis  extends  the  research  accomplished  by  Capt 
Pineiro  and  Lt  Berens  in  the  area  of  adaptive  algorithm 
implementation.  Specifically,  this  thesis  explores  the 
performance  characteristics  of  the  multiple  model  estimation 
algorithm  and  how  they  influence  the  selection  of  aircraft 
models  to  allow  the  parameter  adaptive  control  system  to 
maintain  tracking  performance  over  a  desired  portion  of  the 
flight  envelope.  The  aircraft  dynamic  equations  used  are 
those  of  the  AFTI/F-16  and  the  control  law  design  is  based 
on  the  method  developed  by  Professor  Porter. 

Numerous  fixed  gain  simulations  are  presented  in  order 
to  determine  the  performance  robustness  to  plant  variations 
of  selected  flight  conditions.  From  these  simulations, 
performance  or  robustness  boundaries  for  each  nominal  flight 
condition  are  determined.  Data  is  then  generated  to  ascer¬ 
tain  the  effect  on  the  control  surfaces’  responses  of  re¬ 
placing  the  nominal  model  of  a  performance  boundary  with 
another  model  within  the  same  boundary.  That  leads  to  the 
determination  of  the  critical  factors  to  be  considered  in 


choosing  an  aircraft  model  to  represent  a  specific  set  of 
performance  boundaries. 

Following  the  determination  of  what  factors  are 
important  in  establishing  performance  bounds  for  a  fixed 
gain  system,  model  selection  for  the  multiple  model  algo¬ 
rithm  is  evaluated.  In  attempting  to  select  a  set  of  air¬ 
craft  models  that  would  ensure  an  acceptable  level  of  track¬ 
ing  performance  over  the  desired  flight  envelope,  several 
two-,  three,  and  four-model  configurations  are  examined. 

Each  configuration  is  evaluated  to  determine  the  amount  of 
overlap  of  performance  bounds  required  for  proper  tracking 
performance,  the  amount  of  the  flight  envelope  that  is 
covered,  and  control  surface  performance. 

After  selecting  a  set  of  aircraft  models  that  results 
in  the  best  overall  system  response,  the  effect  of  adjusting 
the  control  law  gains  on  the  performance  of  the  multiple 
model  estimation  algorithm  is  evaluated.  By  assuming  that 
all  states  are  accessible,  sensor  noise  is  then  added  to 
each  of  the  longitudinal  states  to  study  how  noise  impacts 
model  selection.  A  set  of  models  that  produces  acceptable 
tracking  performance  over  the  desired  flight  envelope  and 
the  most  immunity  to  sensor  noise  is  then  selected. 


XXV 


MODEL  SELECTION  FOR  THE  MULTIPLE 
MODEL  ADAPTIVE  ALGORITHM  FOR 
IN-FLIGHT  SIMULATION 


1  .  Introduction 


1 . 1  Overview 

There  are  a  large  number  of  estimation  and  control 
problems  where  algorithm  design  and  implementation  are  de¬ 
pendent  on  the  modeling  and  evaluation  of  uncertain  para¬ 
meters.  One  such  problem  was  analyzed  by  Capt  Luis  Pineiro 
in  his  thesis,  Parameter  Adaptive  Model  Following  For  In 
Flight  Simulation  (1).  Pineiro's  thesis  documented  the  fact 
that  parameter  adaptive  control  techniques  can  be  implemen¬ 
ted  to  allow  a  specific  system  to  maintain  a  desired  level 
of  tracking  performance  in  the  presence  of  plant  parameter 
changes.  Lt  Thomas  Berens  continued  Pineiro’s  research  by 
implementing  a  multiple  model  algorithm  (2)  as  the  primary 
parameter  estimation  technique.  This  estimation  algorithm 
allows  the  system  to  maintain  tracking  performance  when  the 
plant  has  deviated  considerably  from  the  nominal  condition. 
This  thesis  extends  Beren’s  research  by  investigating  the 


performance  characteristics  of  the  multiple  model  algorithm 


and  how  they  affect  the  selection  of  aircraft  models  that 


will  yield  a  desired  level  of  tracking  performance  over  a 
specific  portion  of  the  flight  envelope. 

1.2  Background 

Using  the  techniques  developed  by  Professor  Brian  Porter 
(3;4;5;6;7;8),  Pineiro  developed  an  adaptive,  fast  sampling 
control  law  to  compensate  for  changing  aircraft  parameters 
of  the  Advanced  Fighter  Technology  Integration  F-16  (AFTI/F- 
16)  operating  in  a  model  following  configuration.  This 
technique  involved  the  use  of  on-line,  recursive,  step  re¬ 
sponse  matrix  identifiers  to  update  the  control  law  gains  as 
needed  to  account  for  plant  parameter  variations.  Pineiro’s 
research  was  limited  to  the  linearized,  longitudinal,  rigid 
body  dynamics  of  the  AFTI/F-16  using  perturbation  equations 
of  motion  at  a  nominal  flight  condition  of  10,000  feet 
altitude  and  Mach  0.9  (1). 

In  order  to  keep  the  computational  loading  and  conver¬ 
gence  time  to  a  minimum,  Pineiro  allowed  only  a  small  number 
of  parameters  to  be  estimated  (only  four  parameters  were 
adjustable  out  of  20)  while  retaining  an  accurate  represen¬ 
tation  of  the  plant.  The  remaining  parameters  that  were  not 
among  those  to  be  estimated  were  fixed  at  values  represents- 
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tive  of  the  current  flight  condition.  The  results  of  Pinei¬ 
ro’s  research  documented  the  fact  that  parameter-adaptive 
control  techniques  can  be  implemented  to  allow  a  specific 
system  to  maintain  a  desired  level  of  tracking  performance 
in  the  presence  of  plant  parameter  changes  (produced  by 
large  changes  in  the  stability  derivatives )( 1  )  . 

Capt  Pineiro  recommended  that  analysis  of  this  research 
area  continue  with  the  parameter  adaptive  algorithm  extended 
to  include  a  multiple  model  estimator  (1).  In  his  thesis, 

Lt  Berens  implemented  such  an  estimation  algorithm.  The 
multiple  model  algorithm  implementation  can  be  approached  by 
assuming  that  the  uncertain  parameters  take  on  discrete 
values  within  a  possible  range  of  parameter  values.  For 
each  flight  condition  a  complete  system  model  is  designed 
which  produces  an  estimate  of  the  selected  parameters.  Fol¬ 
lowing  the  design  of  a  bank  of  system  models,  each  model  is 
weighted  by  the  probability  of  that  particular  model  being 
correct  and  averaged  with  the  outputs  of  all  other  models 
(9)  . 

1.3  Problem  Statement 

The  purpose  of  this  thesis  is  to  extend  the  research 


accomplished  by  Capt  Pineiro  and  Lt  Berena  in  the  area  of 
adaptive  algorithm  implementation.  Specifically,  this  the¬ 
sis  explores  the  performance  characteristics  of  the  multiple 
model  estimation  algorithm,  how  they  influence  the  selec¬ 
tion  of  aircraft  models,  and  how  well  they  allow  the  system 
to  maintain  tracking  performance  over  a  desired  portion  of 
the  flight  envelope. 

1.4  Summary  of  Current  Knowledge 

This  section  summarizes  the  current  knowledge  that  is 
pertinent  to  this  thesis  effort.  The  areas  of  in-flight 
simulators,  adaptive  control,  and  multiple  model  adaptive 
estimation  are  examined  in  the  following  paragraphs. 

1.4.1  In-Flight  Simulators.  An  in-flight  simulator  is  an 
aircraft  that  is  able  to  fly  with  the  same  characteristics 
as  another  aircraft.  While  flying  the  simulator,  the  pilot 
actually  experiences  the  flying  qualities  and  feel  of  the 
aircraft  that  is  being  simulated  (10). 

The  in-flight  simulator  has  become  an  extremely  impor¬ 
tant  tool  in  the  aircraft  research  and  development  process. 
As  the  cost  of  developing  aircraft  has  increased,  the  bene¬ 
fits  of  in-flight  simulation  have  become  more  pronounced. 


In-flight  simulation  provides  the  capability  to  fly  aircraft 
systems  prior  to  committing  to  the  full  scale  or  production 
phases  of  the  acquisition  process.  Another  area  in  which 
the  simulator  has  been  extremely  helpful  is  the  area  of 
experimental  testing  of  aircraft.  Investigations  into 
flight  control  systems,  aircraft  flying  qualities,  and  sys¬ 
tem  integration  have  been  addressed  in  detail  with  the 
advent  of  the  in-flight  simulator  (11). 

An  advanced  in-flight  simulator  designated  as  the  Va- 
r iable-Stabi lty  In-Flight  Simulator  Test  Aircraft  (VISTA) 
will  be  developed  by  the  USAF .  VISTA  will  be  a  high  perfor¬ 
mance,  supersonic,  six  degree  of  freedom  simulator  based  on 
a  two  place  fighter  aircraft.  By  having  a  variable  stabili¬ 
ty  system,  VISTA  will  be  able  to  support  a  diverse  range  of 
aeronautical  research  in  the  fields  of  flight  controls, 
handling  qualities,  and  flying  characteristics.  It  will 
also  have  the  ability  to  simulate  current,  modified,  or  new 
aircraft  (11). 

Because  VISTA  will  be  able  to  operate  over  a  large  range 
of  flight  conditions,  wide  variations  in  the  aircraft  sta¬ 
bility  derivatives  will  be  possible.  Jones  and  Porter  have 
shown  that  degradation  in  closed  loop-system  behavior  will 
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inevitably  occur  in  the  case  of  large  plant  parameter  chan¬ 
ges  and  a  fixed  control  law  (3).  To  ensure  that  the  simula¬ 
tion  is  not  degraded  substantially  due  to  inaccurate  stabil¬ 
ity  derivatives,  a  model-following  mode  of  operation  may  be 
used.  The  term  model- f ol lowing  means  that  the  response 
characteristics  of  the  host  aircraft  are  tailored  to  'fol¬ 
low’  the  responses  of  the  simulated  aircraft  as  generated  by 
a  computer  model  (12). 

1.4.2  Adaptive  Control.  Hartman  and  Krebs  (13)  have  shown 
that  adaptive  control  is  a  method  of  maintaining  desired 
tracking  performance  throughout  the  flight  envelope.  In  his 
research,  Pineiro’s  main  objective  was  to  implement  a  para¬ 
meter  adaptive  control  law  that  would  compensate  for  the 
degraded  simulation  performance  that  occurs  when  the  sta¬ 
bility  derivatives  of  the  host  aircraft  take  on  values  which 
have  deviated  considerably  from  the  expected  nominal  values 
(  1  )  . 

The  area  of  adaptive  control  has  become  an  area  of  great 
current  interest  in  automatic  control  research.  The  reason 
for  this  is  that  adaptive  control  is  a  technique  that  pro¬ 
poses  solutions  when  the  system  varies  with  time  in  a  manner 
that  is  not  totally  specifiable  a-priori.  This  is  in 
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contrast  to  the  manner  in  which  most  automatic  control 


techniques  afford  solutions  to  time  invariant  systems  (14). 

The  main  benefit  of  adaptive  control  is  that  the  system 

parameters  can  be  adjusted  on-line  as  the  system  environment 

changes.  Hagglund  (14)  stated  in  his  doctoral  thesis,  New 

Estimation  Techniques  For  Adaptive  Control, 

The  idea  of  adaptive  control  is  to  automate  the 
analysis  and  design.  Hence,  an  adaptive  controller 
does  not  only  control  the  process,  but  it  also 
collects  information  about  the  process  behaviour 
from  the  control  signals  and  the  measured  output 
signals.  Based  on  this  information,  or  analysis, 
the  controller  parameters  are  adjusted  on-line. 


There  are  numerous  schemes  reported  in  the  litera¬ 
ture  for  parameter  adaptive  control,  three  of  which  are  gain 
scheduling,  model  reference  control,  and  self  tuning  regula¬ 
tors  (15). 

The  actual  concept  of  gain  scheduling  started  with  the 
development  of  flight  control  systems.  In  fact,  it  is  the 
predominant  method  used  in  current  aircraft  to  handle  varia¬ 
tions  of  parameters  in  flight  control  systems.  Mach  number 
and  dynamic  pressure  are  measured  and  used  as  scheduling 
variables.  When  the  scheduling  variables  have  been  chosen, 


the  control  parameters  are  evaluated  at  a  number  of  opera- 


ting  conditions  using  an  acceptable  design  method.  The  draw¬ 
backs  of  gain  scheduling  are  that  it  is  an  open  loop  compen¬ 
sation  method  (no  feedback  is  provided  to  correct  for  an 
incorrect  schedule,  resulting  in  gains  being  implemented 
which  could  degrade  system  performance  or  even  result  in 
system  instability)  and  the  design  is  time  consuming  (15). 

Model  reference  adaptive  controllers  try  to  obtain  a 
control  behavior  close  to  that  of  a  reference  model  for  a 
given  input  signal.  The  disadvantage  of  this  type  of  adap¬ 
tive  control  is  that  if  a  fixed  reference  model  is  used,  the 
system  will  approach  the  behavior  of  the  reference  model  as 
opposed  to  the  behavior  of  a  proper  model  at  the  current 
operating  condition.  The  advantage  of  the  model  reference 
system  is  that,  when  the  operating  condition  is  near  the 
fixed  reference  model  and  the  input  is  one  for  which  the 
system  has  been  designed,  quick  adaptation  is  possible 
(  16)  . 

Another  method  of  adaptively  adjusting  the  parameters  of 
a  system  is  to  use  a  self  tuning  regulator  (15),  also  refer¬ 
red  to  as  a  parameter  adaptive  controller.  It  should  be 
noted  that  self  tuning  regulators  are  one  form  of  a  parame¬ 
ter  adaptive  controller  and  are  not  equivalent  to  that 


entire  class.  A  block  diagram  of  the  self  tuning  regulator 
used  by  Pineiro  is  shown  in  Figure  1-1  (1).  The  adaptive 
controller  has  three  different  parts,  the  first  of  which 
(see  Figure  1-1)  is  the  ordinary  feedback  loop  formed  by  the 
control  algorithm  (regulator)  and  the  aircraft  dynamics 
(process).  The  input  and  output  of  the  process  are  fed  into 
a  parameter  estimator  which  performs  on-line  estimation  of 


Figure  1-1.  Parameter  Adaptive  Control  System  (1) 


the  parameters  of  the  process.  The  estimated  parameters  are 
then  used  to  modify  the  controller  gains  by  a  design  calcu¬ 
lation.  The  term  adaptation  mechanism  refers  to  the  parame¬ 
ter  estimator  and  the  design  calculation  (14).  The  host 
aircraft  dynamics  for  the  adaptive  controller  designed  by 
Pineiro  were  the  linearized,  longtitudinal ,  rigid  body  dyna- 


mics  of  the  Advanced  Fighter  Technology  Integration  F-16 
(AFTI/F-16).  The  signals  that  were  used  as  tracking  refer¬ 
ences  (Y  model  in  Figure  1-1)  were  the  elevators  and  flaper- 
ons  and  the  adaptive  system  outputs  (Y  host  in  Figure  1-1) 
were  flight  path  angle  and  pitch  rate  (1). 

The  control  algorithm  and  controller  design  that  were 
implemented  by  Pineiro  were  the  techniques  developed  by 
Professor  Brian  Porter  ( 3 ; 4  ;  5  ;  6  ;  7  ; 8  )  .  Porter  and  Jones  have 
shown  that  on-line  recursive  identifiers  can  be  used  to 
provide  updated  step  response  matrices.  Inclusion  of  these 
updated  step  response  matrices  in  a  digital  proportional 
plus  integral  control  law,  results  in  highly  effective  adap¬ 
tive  digital  set  point  tracking  controllers  for  multivaria¬ 
ble  plants  (3).  Using  these  design  techniques,  Pineiro 
developed  an  adaptive  control  law  to  compensate  for  changing 
aircraft  parameters  of  the  AFTI/F-16  operating  in  a  model 
following  configuration.  This  technique  involved  the  use  of 
on-line,  recursive,  step  response  matrix  identifiers  to 
update  the  control  law  gains  as  needed  to  account  for  plant 
parameter  variations. 

As  shown  in  Figure  1,  the  parameter  adaptive  controller 


also  contains  a  recursive  parameter  estimator.  Many  estima- 
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tion  schemes  have  been  used,  including  stochastic  approxima¬ 
tion,  least  squares,  extended  and  generalized  least  squares, 
instrumental  variables,  extended  Kalman  filtering,  and  the 
maximum  likelihood  method  (15).  A  summary  of  the  different 
estimation  methods  is  given  in  Ljung  and  Sodestrom  (17). 

The  parameter  identification  algorithm  used  by  Pineiro  was  a 
modified  recursive  least  squares  algorithm  developed  by 
Hagglund  ( 1 ) . 

1.4.3  Multiple  Model  Adaptive  Estimation  ( MMAE ) .  In  his 
research,  Berens  continued  Pineiro’s  effort  by  extending  the 
parameter  adaptive  algorithm  to  include  a  multiple  model 
estimator  (2).  MMAE  is  especially  beneficial  when  the  ap¬ 
plication  has  parameters  that  can  undergo  large  changes  such 
as  is  the  case  with  the  stability  derivatives  of  the  host 
aircraft  of  an  in-flight  simulator  (18).  Multiple  model 
adaptation  has  been  successfully  applied  to  numerous  prob¬ 
lems  including  flight  control  problems  (19). 

Multiple  model  implementation  can  be  approached  by  as¬ 
suming  that  the  uncertain  parameters  can  take  on  discrete 
values  within  a  possible  range  of  parameter  values.  A 
complete  system  model  is  designed  for  each  flight  condition, 
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resulting  in  a  bank  of  n  system  models.  Based  upon  the 
characteristics  of  the  residuals  generated  for  each  of  the  n 
system  models,  the  correctness  of  each  of  the  estimated 
parameters  is  calculated  iteratively  from  the  conditional 
probability  for  each  of  the  n  system  models  being  the  best 
model  to  use.  The  conditional  probability  is  determined 
using  the  past  measurement  history.  The  closer  the  value  of 
the  conditional  probability  is  to  one,  the  more  weight  is 
placed  on  that  given  model’s  parameter  estimates  being  cor¬ 
rect.  The  parameter  estimate  of  each  system  model  is  weight¬ 
ed  by  the  value  of  its  probability,  and  the  adaptive  parame¬ 
ter  estimate  is  determined  by  the  probabilistically  weighted 
average  of  all  the  system  model  outputs  {9;  18). 


# 

1.5  Assumptions 

The  assumptions  in  this  thesis  are  those  typically 

e 

adopted  when  using  linearized  longitudinal  perturbation 
equations.  These  assumptions  are  those  used  by  Pineiro  and 
Berens  and  are  listed  below. 


*  The  aircraft  is  a  rigid  body  with  constant  mass. 

*  The  atmosphere  is  fixed  with  respect  to  the  earth. 
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*  The  earth* 8  surface  is  an  inertial  reference  frame. 

*  Linearization  about  an  operating  condition  is  valid 
for  preliminary  aircraft  models. 

*  Aerodynamics  are  constant  for  Mach  and  altitude. 

1 . 6  Approach 

Using  the  control  law  design  parameters  and  commanded 
inputs  of  flight  path  angle  and  pitch  rate  as  given  by 
Pineiro,  numerous  fixed  gain  simulations  were  accomplished 
in  order  to  determine  the  performance  robustness  to  plant 
variations,  of  selected  flight  conditions.  From  these  simu¬ 
lations,  performance  or  robustness  boundaries  for  each  nomi¬ 
nal  flight  condition  were  determined.  Data  was  then  gener¬ 
ated  to  ascertain  the  effect  on  the  control  surfaces’  re¬ 
sponses  of  replacing  the  nominal  model  of  a  performance 
boundary  with  another  model  within  the  same  boundary  on  the 
control  surfaces  responses.  That  led  to  the  determination  of 
the  critical  factors  to  be  considered  in  choosing  an  air¬ 
craft  model  to  represent  a  specific  set  of  performance 
bounds . 

Following  the  determination  as  to  what  factors  were 
important  in  establishing  performance  bounds  for  a  fixed 


gain  system,  model  selection  for  the  multiple  model  algo- 


rithm  was  evaluated.  In  attempting  to  select  a  set  of 
aircraft  models  that  would  ensure  an  acceptable  level  of 
tracking  performance  over  the  desired  flight  envelope,  seve¬ 
ral  two-,  three-,  and  four-model  configurations  were  exam¬ 
ined.  (The  term  model  configuration  refers  to  the  number  of 
models  that  were  explicitly  implemented  in  the  multiple 
model  algorithms).  Each  configuration  was  evaluated  to 
determine  the  amount  of  overlap  of  performance  bounds  re¬ 
quired  for  proper  tracking  performance,  the  region  of  the 
flight  envelope  that  was  covered,  and  control  surface  per¬ 
formance  . 

After  selecting  a  set  of  aircraft  models  that  resulted 
in  the  best  overall  system  response  (discussed  in  a  later 
section),  the  effect  of  adjusting  the  control  law  gains  on 
the  performance  of  the  multiple  model  estimation  algorithm 
was  evaluated.  By  assuming  that  all  states  were  accessible, 
sensor  noise  was  then  added  to  each  of  the  longitudinal 
states  to  study  how  noise  impacts  model  selection.  A  set  of 
models  that  produced  acceptable  tracking  performance  over 
the  desired  flight  envelope  and  the  most  immunity  to  sensor 
noise  was  then  selected. 
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1  . 7  Overview 

The  material  in  this  thesis  begins  in  Chapter  2  with  a 
description  of  the  AFTI/F-16  and  a  presentation  of  the 
mathematical  aircraft  model.  The  actuator  model  is  also 
discussed  in  Chapter  2.  Chapter  3  contains  a  brief  descrip¬ 
tion  of  the  multivariable  design  theory  developed  by  Profes¬ 
sor  Brian  Porter  of  the  University  of  Salford,  England.  The 
multiple  model  estimation  algorithm  is  discussed  in  Chapter 
4  and  model  selection  criteria  are  presented  in  Chapter  5. 
Chapter  6  presents  the  simulation  results,  and  a  discussion 
of  conclusions  and  recommendations  for  future  study  is  found 
in  Chapter  7 . 
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2.  Aircraft  Description  and  Models 

2.1  AFTI/F-16  Description 

Because  the  VISTA  is  still  in  the  process  of  being 
developed  and  the  aircraft  dynamic  equations  are  not  avail¬ 
able,  the  Advanced  Fighter  Technology  Integration  (AFTI)/F- 
16  was  used  as  the  host  aircraft  for  the  parameter  adaptive 
control  system.  The  reason  for  this  choice  is  the  antici¬ 
pated  similarity  of  the  AFTI  and  proposed  VISTA.  The 
AFTI/F-16  aircraft,  shown  in  Figure  2-1,  is  an  F-16A  air 
superiority  fighter  modified  to  be  a  testbed  for  evaluating 
new  aircraft  technologies  (20).  These  modifications  include 
the  addition  of  two  vertical  canards  which  are  mounted  on 
the  engine  inlet,  control  surfaces  that  allow  independent 
motion  of  the  trailing  edge  flaps  and  independent  motion  of 
the  horizontal  tail  halves,  and  a  redundant,  digital  fly-by- 
wire  flight  control  system.  The  AFTI/F-16's  mission  is  very 
important  in  that  the  technologies  tested  on  it  may  be 
incorporated  in  the  Air  Force’s  future  fighter  aircraft. 

The  unaugmented  AFTI/F-16  aircraft  is  statically  unsta¬ 
ble  in  the  longitudinal  axis  for  subsonic  flight.  The 
reason  this  condition  exists  is  that  the  center  of  gravity 
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Figure  2-1.  The  AFTI/F-16 


is  located  behind  the  aircraft's  center  of  pressure.  The 
instability,  manifested  by  an  unstable  short  period  mode, 
allows  the  aircraft  to  withstand  higher  load  factors  and 
reduces  drag, thereby  enhancing  maneuverability.  Additional 
ly,  the  aircraft  has  a  lightly  damped  dutch  roll  mode. 
Therefore,  the  flight  control  system  has  a  twofold  purpose: 
to  stabilize  the  aircraft  longitudinally  and  to  improve  the 
dutch  roll  damping  (20). 

In  addition  to  being  able  to  perform  conventional  ma¬ 
neuvers  such  as  pitching  longitudinally,  rolling  laterally, 
and  turning  with  zero  sideslip,  the  AFTI/F-16  can  also 
perform  unconventional  maneuvers  that  require  decoupling  of 
the  aircraft’s  modes  of  responses.  These  unconventional 
maneuvers  include  pitch-pointing,  yaw-pointing,  and  lateral 


and  longitudinal  translation.  Conventional  and  decoupled 
maneuvers  can  also  be  blended,  which  gives  the  AFTI/F-16  an 
even  wider  array  of  possible  maneuvers  (20). 

The  control  surfaces  that  are  used  to  maneuver  the 
AFTI/F-16,  as  well  as  the  definitions  of  positive  surface 
deflections,  are  shown  in  Figure  2-2.  The  surfaces  include 
left  and  right  horizontal  tail  halves,  left  and  right  flape- 
rons,  left  and  right  canards,  and  the  rudder.  The  horizon¬ 
tal  tail  halves  have  a  dual  function  in  that  they  can  be 
deflected  symmetrically  as  elevators  to  pitch  the  aircraft, 
or  they  can  be  deflected  asymmetrically  to  augment  rolling. 
Likewise,  the  flaperons  can  be  deflected  symmetrically  to 
function  as  conventional  flaps  or  asymmetrically  to  function 
as  ailerons  to  control  roll.  The  canards  can  be  used  in  a 
snow  Pi  ow  configuration  to  function  as  a  speedbrake  or  they 
can  be  used  to  provide  sideforces  on  the  aircraft.  As  in  a 
conventional  aircraft,  the  rudder  is  used  for  yawing.  The 
surfaces  used  in  this  thesis  are  shown  in  heavy  black  in 
Figure  2-2  and  include  the  left  and  right  tail  halves  and 
the  left  and  right  flaperons. 


2.2  Aircraft  Models 

2.2.1  State  Space  Model.  The  aircraft  is  modeled  by  a  set 
of  first  order  differential  equations  in  the  state  space 
form 

x=Ax+Bu  (2-1) 

The  state  equations  are  perturbation  equations  about  a  trim 
condition  and  are  derived  from  the  aircraft  equations  of 
motion  which  consist  of  forces  at  the  center  of  gravity  and 
moments  with  reference  to  the  body  axis.  Detailed  deriva¬ 
tions  of  the  longitudinal  state  perturbation  equations  used 
in  this  thesis  are  presented  in  a  thesis  by  Barfield  (21). 

Equation  (2-2)  shows  the  longitudinal  axis  state  space 
model  for  the  AFTI/F-16.  The  primed  terms  in  Equation  (2-2) 
are  dimensional ized  derivatives  (Appendix  A)  in  the  body 


e 

0 

0 

0 

l 

0 

0 

0 

u 

xe 

x' 

u 

xa 

x' 

q 

U 

+ 

X6e 

x6f 

a 

ze 

z' 

u 

Zoc 

z' 

q 

a 

Z6e 

z8f 

q 

Me 

M' 

u 

M' 

a 

M' 

q 

q 

M«e 

M6f 

8e 
8  f 


(2-2 


where 


6  is  pitch  angle 


u  is  forward  velocity 
&  is  angle  of  attack 
q  is  pitch  rate 

5  is  the  total  elevator  where  both  horizontal 
fail  halves  deflect  symmetrically 

5^,  is  the  total  flap  where  both  flaperons  deflect 
symmetrically 

axis.  The  definition  of  the  body  axis  along  with  the  posi¬ 
tive  directions  of  forces,  moments,  and  angles  are  found  in 
Figure  2-3. 

A  Flight  Dynamics  Laboratory  aerodynamic  data  package 
for  the  AFTI/F-16  was  used  to  obtain  linearized  aerodynamic 
coefficients  and  dimensionalized  derivatives  by  trimming  the 
aircraft  at  various  flight  conditions.  Aircraft  weight  and 
inertias  were  common  to  all  flight  conditions  and  were  also 
provided  by  the  program  (see  Appendix  A). 

This  thesis  is  concerned  with  that  portion  of  the 
flight  envelope  which  ranges  in  altitude  from  10,000  to 
38,000  feet  and  from  0.25  to  0.9  in  Mach  number.  This 
portion  of  the  flight  envelope  is  shown  graphically  in 
Figure  2-4.  The  specific  flight  conditions  considered  in 
the  envelope  of  interest  are  listed  in  Table  2-1.  The 
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Flight  Conditions  of  Interest 


It  ( kf t ) /Mach  No 


Alt  ( kf t ) /Mach  No, 


38 

0.45 

26 

0.40 

38 

0.50 

26 

0.45 

38 

0.55 

26 

0.50 

38 

0.60 

26 

0.55 

38 

0.65 

26 

0.60 

38 

0.70 

26 

0.65 

38 

0.75 

26 

0.70 

38 

0.80 

26 

0.75 

38 

0.85 

26 

0.80 

38 

0.90 

26 

0.85 

34 

0.40 

26 

0.90 

34 

0.45 

22 

0.30 

34 

0.50 

22 

0.35 

34 

0.55 

22 

0.40 

34 

0.60 

22 

0.45 

34 

0.65 

22 

0.50 

34 

0.70 

22 

0.55 

34 

0.75 

22 

0.60 

34 

0.80 

22 

0.65 

34 

0.85 

22 

0.70 

34 

0.90 

22 

0.75 

30 

0.35 

22 

0.80 

30 

0.40 

22 

0.85 

30 

0.45 

22 

0.90 

30 

0.50 

18 

0.30 

30 

0.55 

18 

0.35 

30 

0.60 

18 

0.40 

30 

0.65 

18 

0.45 

30 

0.70 

18 

0.50 

30 

0.75 

18 

0.55 

30 

0.80 

18 

0.60 

30 

0.85 

18 

0.65 

30 

0.90 

18 

0.70 

26 

0.30 

18 

0.75 

26 

0.35 

18 

0.80 

Alt  (kft)/Mach  No, 


0.25 

0.30 

0.35 

0.40 

0.45 

0.50 

0.55 

0.60 

0.65 

0.70 

0.75 

0.80 

0.85 

0.90 

0.25 

0.30 

0.35 

0.40 

0.45 

0.50 

0.55 

0.60 

0.65 

0.70 

0.75 

0.80 

0.85 

0.90 
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longitudinal  d imens ional i zed  derivatives  for  each  of  the 
flight  conditions  listed  in  Table  2-1  have  been  obtained 
using  the  Flight  Dynamics  Laboratory’s  AFTI/F-16  aerodynamic 
data  package.  A  state  space  model  for  each  flight  condition 
is  formed  by  inserting  the  proper  dimensionalized  deriva¬ 
tives  into  Equation  (2-1).  Appendix  A  presents  the  state 
space  models  for  each  of  the  flight  conditions. 

An  output  vector,  y,  is  defined  for  the  model  in  Equa¬ 
tion  (2-1)  in  order  to  command  desired  maneuvers.  The 
output  vector  is  added  to  the  state  space  aircraft  model 
with  an  output  matrix,  C.  The  aircraft  model  state  and 
output  equations  then  take  the  form  of 


x  =  Ax  +  Bu 


y  =  Cx 


(2-3) 


The  output  vector  used  in  this  study  is  the  same  vector  as 
that  used  by  Pineiro  (1)  and  Berens  (2)  in  which  flight  path 
angle  and  pitch  rate  are  commanded.  Flight  path  angle  is 
given  by 


7  =  ©  -  a 


(  2-4  ) 


Because  pitch  angle,  pitch  rate,  and  angle  of  attack 


are  state  variables,  flight  path  angle  and  pitch  rate  are 


easily  included  in  the  model  by  placing  unit  elements  with 
the  proper  signs  in  the  appropriate  positions  in  the  C 
matrix.  This  results  in  an  output  equation 
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(2-5) 


2 . 2 .2 _ Autoregressive  Difference  Equation  Model.  The 

AFTI/F-16  plant  state  and  output  equations  of  Equation 
(2-3)  may  be  discretized  for  a  sampling  period  T  by  using 
the  relationships 
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The  resulting  sampled  data  state  and  output  equations  become 
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y(kT)  =  [c 


x  ( kT ) 
x  (kT) 


(2-9) 


Pineiro  showed  (Appendix  B)  that  the  input-output  relation¬ 
ship  of  Equations  (2-8)  and  (2-9)  can  be  expressed  in  terms 
of  an  Nth  order  autoregressive  vector  difference  equation  of 
the  form 


y(kT)  =  BjUflk-lIT)  -  A^Mk-DT)  +  .  (2-10) 

+  B  u  {  (  k-N  )T)  -  A  y  (  <  k-N  )  T  )  +  e(kT) 

n  n 

or  equivalently 

y  (  kT )  =Te(kT)  +  e(kT)  (2-11) 

where  the  equation  error  vector  e(kT)  is  assumed  to  be  a 

zero-mean  Gaussian  white-noise  vector  with  covariance  A  , 

_e 

Tfk'1’)  6  RmX^  is  a  matrix  of  past  values  of  y  and  u,  the 
mxm  mxm 

matrices  A .  6  R  (1=1,2, B,  6  R  (  i  =  1 , 2 . ,N) 
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and  the  vector  6  €  R  are  the  parameters  of  the  Nth  order 


difference  equation,  m  is  equal  to  the  number  of  inputs  and 


(  is  equal  to  the  number  of  parameters  in  the  difference 
equation.  For  the  case  AFTI/F-16  longitudinal  aircraft 
model  of  Equation  (2-3),  the  parameter  vector  ©  consists 
of  a  column  vector  of  20  elements  which  are  unique  for  a 
given  flight  condition. 


2.3  Actuator  Model 

Consideration  of  control  surface  position  and  rate 
limits  is  important  in  the  design  of  a  realistic  control 
system.  The  actuator  transfer  function  model  for  both  the 
elevator  and  the  flaperon  surfaces  are  first  order  models  of 
the  form 


T 

actuator 


(  s  ) 


44 


(  s  +  44  ) 


(2-12) 


The  control  surface  position  and  rate  limits  used  in  this 
study  are  those  given  in  Pineiro’s  thesis  and  are  found  in 
Table  2-2.  The  position  limits  assume  that  the  elevator  and 
flaperon  have  a  trimmed  position  of  zero  degrees.  Due  to 
the  fact  that  each  flight  condition  has  a  different  nonzero 
trim  value  of  elevator  and  flaperon  position,  the  surface 


position  limits  are  adjusted  accordingly  for  each  point. 


Table  2-2. 
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Control  Surface  Position  and  Rate  Limits 


Surface 

Position  Limit 

Rate  Limit 

(deg  ) 

( deg/sec ) 

Elevators 

+  25 

90 

Flaperons 

+  20 

78 

-23 

The  position  limits  for  each  flight  condition  are  found  in 
Appendix  C. 

2 . 4  Summary 

This  section  presents  the  mathematical  development 
of  the  equations  describing  the  AFTI/F-16.  The  AFTI/F-16 
is  chosen  as  the  plant  due  to  the  non-availibili ty  of 
similar  equations  for  the  VISTA.  The  autoregressive  differ¬ 
ence  equation  representing  an  equivalent  state  space  repre¬ 
sentation  of  the  aircraft  is  also  presented,  along  with  the 
flight  conditions  considered  in  this  thesis. 
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3 . 1  Introduction 


This  chapter  describes  the  control  algorithm  that  is 
implemented  in  the  parameter  adaptive  control  system  as 
shown  in  Figure  1-1.  The  control  law  should  be  such  that 
the  plant  is  able  to  follow  the  time  response  generated  by  a 
computer  model  of  the  vehicle  being  simulated.  The  multi- 
variable  control  algorithm  used  in  this  thesis  is  that 
developed  by  Professor  Brian  Porter  of  the  University  of 
Salford,  England.  Porter’s  method  employs  output  feedback 
with  high-gain  error-actuated  controllers.  Output  feedback 
may  be  advantageous  since  state  variables  may  be  difficult 
to  measure. 

3.2  Model  Following  Technique  (1) 

One  method  of  achieving  in-flight  simulation  is  through 
the  implementation  of  the  model  following  concept  (12).  A 
model  following  in-flight  simulation  scheme  uses  the  signals 
from  the  evaluation  pilot’s  cockpit  as  inputs  to  an  aircraft 
computer  which  simulates  the  equations  of  motion  of  the 
desired  aircraft.  The  output  of  the  same  aircraft  computer 


is  comprised  of  the  time  histories  that  describe  the  respon¬ 
ses  of  the  simulated  aircraft  to  the  inputs  from  the  evalua¬ 
tion  pilot.  The  variable  stability  system  has  the  responsi¬ 
bility  of  cianuevering  the  control  surfaces  of  the  host 
aircraft  so  that  the  aircraft  response  at  the  pilot’s  sta¬ 
tion  is  the  same  as  those  of  the  emulated  airplane.  Obvi¬ 
ously,  achieving  satisfactory  model  following  performance 
depends  on  having  an  accurate  knowledge  of  the  host  air¬ 
craft’s  stability  derivatives. 

The  dynamics  of  the  host  vehicle  can  be  represented  by 
a  differential  equation  of  the  form 

• 

x^(t)  =  A^x^ft)  +  B^u^  ( t )  (3-1) 

and  the  dynamics  of  the  vehicle  to  be  simulated  can  be 
represented  by 

*  <t)  "  +  ^ny»(t)  <3~2) 

m  in  in  mm 

where  A,  ,  A  ,  B.  ,  B  are  matrices  that  contain  the  stability 
“h  ~m  “h  ~m 

derivatives,  x.  and  x  are  the  states  of  the  respective 
h  "Tn 

differential  equations,  u  and  u  are  the  required  control 

n  — m 

surface  deflection,  and  the  subscripts  h  and  m  signify  the 
host  and  model  aircraft,  respectively. 


c 


By  substituting  the  states  and  the  rate  of  change  of 
the  states  of  the  model  into  Equation  (3-1),  a  control  law 


for  exact  model  following  can  be  obtained  by  solving  for  the 


control  input  u,  ( t )  : 

h 


x  (t)  =  A^x  (t)  +  JLjl  (t) 

Ti  riTn  n  h 


(3-3) 


Assuming  B,  ie  invertible  and  solving  for  \y  (t)  leads  to 
~h  “ h 


u  (t)  =  B  1x  (t)  -  B  1A  x  (t) 
h  h  m  h  h  m 

u  (t)  =  Kx  (t)  -  K  x  (t) 
h  1  m  ~2  m 


(3-4  ) 
(3-5) 


From  Equations  (3-4)  and  (3-5)  it  can  be  seen  that  the 
control  inputs  to  the  host  aircraft  are  defined  in  terms  of 
the  states  of  the  model  and  the  first  derivative  of  those 
states.  Also,  Equations  (3-4)  and  (3-5)  show  that  accurate 
stability  derivatives  are  required  to  have  good  model  fol¬ 
lowing  performance.  As  shown  in  Figure  3-1,  the  states  and 


rates  of  change  of  the  states  are  available  from 


the  variable  stability  system’s  computer  that  contains  the 
model  aircraft  dynamics. 

Figure  3-2  shows  an  improved  explicit  model  following 
system.  Notice  that  there  is  a  feedback  loop  containing  a 
gain  matrix  around  the  plant.  The  purpose  of  this  feed- 


Figure  3-2.  Improved  Explicit  Model  Following  System  (1) 

back  loop  is  to  reduce  the  sensitivity  of  the  explicit  model 
following  system  (see  Figure  3-1)  to  plant  parameter  varia¬ 
tions  . 

With  the  inclusion  of  the  feedback  loop  the  feedforward 


gains  become 


c 


* 


c 


% 


If 


Kl  >> 

— h  “h  — h 


(3-8) 


then  the  feedback  loop  becomes  increasingly  tighter  which 
makes  the  system  less  dependent  on  (less  sensitive  to 
plant  parameter  variations). 


3.3  Porter's  Control  Law  (1) 


The  multivariable  control  law  used  in  this  thesis  to 
achieve  model  following  was  developed  by  Professor  Brian 
Porter.  Porter  developed  a  control  law  that  can  be  used  to 
make  the  plant  follow  the  responses  generated  by  a  computer 
model  of  the  vehicle  to  be  simulated.  This  section  presents 
Porter’s  control  law. 

The  host  aircraft  in  this  thesis  is  completely  con¬ 
trollable  and  observable,  and  is  described  by  the  continuous 
time  state  space  model 

x(t)  =  A  x(t)  +  B  u(t)  (3-9) 

Z(t)  =  C  x ( t )  (3-10) 

where 

A  =  continuous  plant  matrix  (n  x  n) 


<• 
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B  =  continuous  input  control  matrix  (n  x  m) 
having  a  rank  "m" 

C  =  continuous  output  matrix  (  £  x  n) 
y  =  output  vector  with  (  outputs 
u  =  input  vector  with  m  inputs 
x  =  state  variable  vector  with  n  states 

For  the  linearized,  longitudinal  dynamics  of  the  AFTI/F-16, 
the  state  variable  vector  has  4  states  as  described  in 
Chapter  2: 

6  =  pitch 

u  =  forward  velocity 
O'  =  angle  of  attack 
q  =  pitch  rate 

The  input  vector  consists  of  2  inputs 

A  =  elevator 
e 

6^.  =  flaperons 

The  2  elements  of  the  output  vector  are 

y  =  flight  path  angle 
q  =  pitch  rate 
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In  using  Porter's  design  method,  it  is  desirable  to 


partition  the  state  equations  so  that  B  and  C  are  square 

fa  fa« 

(m  x  m)  and  (  f  x  f )  matrices,  respectively.  This  process 
yields 


and 


x  ( t ) 


x2(t) 


- 

- 

- 

—  “ 

-12 

x1  ( t ) 

+ 

0 

“2  1 
_ 

— 2  2 

x2(t) 

h 

( 3-m 


y(  t ) 


x  ( t ) 

V11 


(  3-12  ) 


where  it  is  required  that  the  number  of  inputs  to  the  system 
equal  the  number  of  outputs,  m=£.  This  requirement  ensures 
that  the  column  dimension  of  B  equals  the  row  dimension  of 
C.  In  equation  (3-12)  the  elements  associated  with  x ^ ( t ) 
and  x^ ( t )  are 

x  (t)  =  {  e  ]  (3-13) 

x  (t)  =  fu  Qfq]  (3-14) 


The  plant  of  Equations  (3-11)  and  (3-12)  is  defined  as 
regular  if  the  matrix,  CB,  has  full  rank  "m" .  The  input 
output  configuration  of  this  thesis  provides  for  a  regular 
design.  For  regular  plants  with  stable  transmission  zeros, 


c 


© 


© 


© 


c 
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the  control  law  is  a  discrete  proportional-plus-integral 
(PI)  output  feedback  control  law  (see  Figure  3-3)  expressed 

as 

u(kT)  =  (1/T)  [KJ  e ( kT )  +  K2z(kT)]  (3-15) 

and 

u(t)  =  u(kT)  for  t  6  [kT,  (k+l)T]  (3-16) 

where 

r(kT)  is  the  sampled  reference  tracking  signal 
k  is  an  integer 

T=  (1/f)  is  the  sampling  period 

Kj  and  are  the  (m  x  m )  controller  matrices 

e(kT)  is  the  error  vector  [ e ( kT ) =r ( kT ) -£( kT ) ] 

z(kT)  is  the  digital  integral  of  the  error  vector 
[z( (k+1 )T)=z(kT)  +  Te(kT)] 

When  CB  does  not  have  full  rank,  the  plant  is 

labeled  as  irregular  and  a  proportional-plus  integral-plus 

derivative  (PID)  controller  structure  may  be  implemented 

(  6  ;  7  )  . 

3.3.1  Fixed  Gain  Controller  Matrices. 

The  continuous  time  state  equations  given  in  Equations 
(3-11)  and  (3-12)  can  be  discretized  for  the  sampling  period 
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System  Block  Diagram 
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The  resulting  sampled  data  state  and  output  equations  for 
the  plant  are 


Xjllk  +  1)T) 

r*  *  1 

*12 

Xj (kT) 

x2((k  +  1 )T) 

§  $ 
L21  -22J 

x2(kT) 

¥ 

— 2 


and 


Z(kT)  =  [C 


V 


[x  (kT  j] 

[vkT>J 


u(kT) 

(3-19) 


(3-20) 


v.. 


The  augmented  closed-loop  state  and  output  equations 
for  the  control  law  of  equation  (3-15)  are  (8) 


z  { (  k  + 1  )  T  f 

I  -TC  -TC 

"  z ( kT ) 

x  (  (k  +  1 )T) 

= 

x  (kT) 

x2( (k+l)T) 

Hh  i-'VA  i22-fW2J 

x2(kT) 

"TI 

fi?i 


r  ( kT )  (3-21) 
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z(kT) 


X(kT)  = 


[0  c  Cl  x  (kT) 
x‘(KT, 


( 3-22  ) 


A  transformation  block  diagonalizes  equations  (3-21) 
and  ( 3-22  )  as 


where 


x i { { k+ 1 ) T }  A 


x  {(k+l)T)  0 


^  (kT) 


r(  kT  > 


A  x  (kT) 

“4  2 


B  (3-23) 
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X(kT)  =  ( 


C  ]  Tx  (kT) 
x*  ( kT  ) 


( 3-24  ) 


x^kT) 

=  z(kT) 

( 3-25  ) 

x  (kT) 

x2<kT) 

=  x  (kT) 

( 3-26  ) 

Ci  =  CK, 

_1K2  ,  o  1 

( 3-27  ) 

h 

=  £2 

( 3-28 ) 
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( 3-29  ) 
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As  the  sampling  frequency  is  increased,  the  closed-loop 


transfer  function  takes  on  the  asymptotic  form  (23) 


G( z)  =  G  <  z)  +  G  (z )  ( 3-33  ) 


where  z  is  the  discrete  transform  operator  and 


G  (z)  =  C  (zl  -I  -TA  )  TB  (3-34) 

~1  —  1  n  n  ~o  “o 

G  ( z  )  =  C0  {  z I  -I  -A  ) ~ 1 B  (3-35) 

“2  ~2  mm  “4  _4 

with 


Kl‘‘K2 

*12S2’  K2 


Aji-A^C-1  C, 


( 3-36  ) 


B 


4!2£2 


a,  =(-82B,c2] 


( 3-37  ) 


( 3-38 ) 


In  Equations  (3-33),  (3-34),  and  (3-35),  G  ( z )  and 

G^  (  z  )  are  the  slow  and  fast  mode  transfer  functions  respec¬ 


tively.  The  slow  modes  contain  two  sets  of  eigenvalues  Z 
and  which  are  given  by  the  equations 


1 


Z,  =  { z  €  C :  det { zl  -I  +TK~1K»  =  0) 
1  “m  ~m  ~1  “2 


( 3-39  ) 


-1 


Zo  =  { z  €  C :  det { zl  -I  -TA,  +  TA,  C  C, )  =  0)  (3-40/ 

2  n-m  “n-m  “11  ~12“2  ~1 
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The  fast  mode  eigenvalues  are  given  by 


Z0  =  (zf  C:  det  { z  I  -I  +C  R  K  )  =  0}  (3-41) 

3  ~m  ~m  2  2  1 

Because  of  the  form  of  A  ,  B  ,  and  C, ,  the  eigenvalues 

o  ~  o  “1 

of  A  are  uncontrollable  or  unobservable.  Therefore,  as  the 

o 

sampling  frequency  increases,  the  slow  transfer  function 
asymptotically  approaches  zero  and  the  overall  system 
transfer  function  contains  only  the  fast  modes,  as  given  in 
G^(z),  which  can  be  expressed  in  the  equivalent  form 

G(z)  =  G  ( z )  =  (zl  -I  +  C  B  K ,)_1C  B  k  (3-42) 

Z  mm  2  2  1  2  2  1 

The  desired  fixed  gain  controller  matrices  and 
are  then  given  by 


K,  =  (W‘S 


”2  = 


( 3-43  ) 


( 3-44  ) 


where  p  is  a  positive  scalar 
the  ratio  of  proportional  to 
diagonal  tuning  matrix.  The 
mine  the  weighting  effect  of 


greater  than  zero 
integral  control, 
diagonal  elements 
a  particular  error 


and  assigns 
and  £  is  a 
of  4  deter- 
signal  on 


the  respective  control  input.  Hence,  selection  of  the  di- 


agonal  weighting  matrix  elements  can  be  altered  to  achieve 


€ 


<5 


e 


<0 


<0 


desired  tracking  characteristics. 


3.3.2  Adaptive  Controller  Matrices.  To  avoid  performance 
degradation  in  the  presence  of  large  plant  parameter  varia¬ 
tions,  Porter  has  shown  <3;4;5;6;7)  that  by  incorporating 
fast  on-line  recursive  identifiers  which  provide  updated 
step-response  matrices,  satisfactory  response  characteris¬ 
tics  can  be  acheived.  The  significance  of  using  the  step 
response  matrix  is  that  it  can  be  obtained  from  real  time 
input-output  measurements  to  reflect  the  current  charac¬ 
teristics  of  the  plant,  thus  forming  the  basis  for  an  adap¬ 
tive  system. 

The  discrete  step  response  matrix  is  defined  as  (4) 

T 

H ( T )  =  j  C  exp(At)  B  dt  (3-45) 

0 

For  small  sampling  periods  H(T)  ’=  TCB,  and  the  control  law 
design  equation  becomes 


u(kT)  =  K  e(kT)  +  LzIkTI  (3-46) 

X  fa 

where 

R  =  H_1(T)|  (3-47) 

K2  =  pgj  (3-48) 
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Equations  (3-47)  and  (3-48)  show  that  for  the  case  of 
plants  with  variable  parameters ,  updated  step  response  ma¬ 
trices  need  to  be  provided  to  the  control-law  design  equa¬ 
tions.  This  requirement  can  be  achieved  due  to  the  fact 
that  the  step  response  matrix  can  be  obtained  from  real  time 
input-output  data  reflecting  the  current  plant.  Thus,  by 
using  the  step  response  matrix,  an  adaptive  control  system 
can  be  designed  to  track  plants  with  variable  parameters. 

The  elements  of  the  step  response  matrix  for  the  control  law 
equation  (3-47)  can  be  obtained  by  expressing  the  state  and 
output  equations  ((3-19)  and  (3-20))  as  an  Nth  order  autore¬ 
gressive  difference  equation,  as  shown  in  Chapter  2,  of  the 
form 


r  r 


L 


y(kT)  =  B^ulk-DT)  -  A^llk-lIT)  +  ... 

+  B  u ( ( k-N )T )  -  A  y  { ( k-N )T )  +  e ( kT )  (3-49) 

n  ~n  _ 


or 


y(kT)  =  T1  (kT)©  +  e  ( kT ) 


( 3-50) 


where  the  equation  error  vector  e(kT)  is  assumed  to  be  a 

„T  mx  1 

zero-mean  Gaussian  white-noise  vector,  T  (kT)€  R  is  the 


measurement  mat 


rix  consisting  of  past  values  of  (y^ 


)  and 
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mxm  mxm 

(u.  },  the  matrices  A.€  R  ( i  =  1 , 2 , . . • N  )  ,  B.e  R 

~k  ~i  _i 

(i=l,2,..N),  and  the  vector  ©  €  R^  are  the  parameters  of 
the  Nth  order  difference  equation. 

It  can  be  shown  that,  by  using  the  definition  of  the 
step  response  matrix,  the  following  equation  can  be  derived 

H(  T  )  =  11  T  TCB  (  3-51 ) 

Equation  (3-51)  clearly  shows  that  by  identifying  in 
real-time  using  input-output  data,  updated  step-response 
matrix  estimates  can  be  provided  for  the  control  law  design 
equation,  (3-47).  From  Equation  (3-49),  it  can  be  seen  that 
B i  can  be  identified  as  a  set  of  parameters  in  the  parameter 
vector  0. 


3 . 4  Summary 

This  chapter  briefly  presents  the  concept  of  model 
following  in-flight  simulation  and  discusses  a  control  law 
design  technique  that  can  be  used  to  accomplish  model  fol¬ 
lowing.  The  design  technique  discussed  is  the  Porter  me¬ 
thod,  and  both  fixed  gain  and  adaptive  controller  gains 
matrices  are  presented.  Finally,  the  method  of  adaptively 
updating  the  controller  gain  matrices  via  the  parameters  of 
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an  autoregressive  difference  equation  are  presented.  The 
next  chapter  discusses  in  detail  the  parameter  estimation 
algorithm  used  in  this  thesis. 


Parameter  Estimation 


4_. 

4 . 1  Introduction 

This  chapter  details  the  algorithm,  used  in  this  the¬ 


sis,  that  allows  the  process  parameters  to  be  identified. 

Due  to  the  fact  that  the  stability  derivatives  change  as  the 
aircraft  operates  through  varying  flight  conditions,  a  para¬ 
meter  identification  algorithm  is  required  that  estimates 
the  time  varying  parameter  vector  9(t).  In  his  thesis, 
Pineiro  incorporated  a  modified  recursive  least  squares 
algorithm  developed  by  Hagglund  (14)  to  account  for  time 
varying  parameters.  In  order  to  minimize  the  number  of 
parameters  to  be  identified,  Pineiro  fixed  some  parameters 
that  were  not  used  in  the  control  law  and  were  assumed  to  be 
constant  about  a  nominal  flight  condition.  Thus,  the  number 
of  parameters  estimated  by  the  estimation  algorithm  was 


reduced  to  f our . 

In  an  attempt  to  decrease  the  parameter  estimation  con¬ 
vergence  time  required  to  estimate  the  complete  set  of  plant 
model  parameters,  Berens  (2)  implemented  a  multiple  model 


"'■A 


algorithm  parameter  estimation  technique  in  lieu  of  the 


modified  recursive  least  squares  algorithm  used  by  Pineiro. 


»*V 
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The  algorithm  discussed  in  this  chapter  and  used  in  this 
thesis  is  the  multiple  model  parameter  estimation  technique 
which  was  presented  by  Wittenmark  (24)  and  implemented  by 
Berens.  The  manner  in  which  models  are  to  be  selected  for 
the  multiple  model  algorithm,  to  provide  satisfactory  track¬ 
ing  performance  in  the  flight  envelope  of  interest,  is  the 
primary  emphasis  of  this  thesis. 

4.2  Multiple  Model  Algorithm 

The  multiple  model  adaptive  algorithm  is  a  parameter 
estimation  technique  which  can  yield  satisfactory  parameter 
estimates  over  a  desired  parameter  space.  By  selecting  a 
discrete  number  of  parameter  values  that  are  dispersed 
throughout  the  region  of  expected  parameter  space  and  incor¬ 
porating  them  as  a-priori  data,  the  multiple  model  estima¬ 
tion  algorithm  provides  a  method  of  discriminating  which  of 
the  a-priori  models  is  nearest  to  the  true  parameter  value. 

The  multiple  model  algorithm  defines  the  parameter 
vector  as  a  random  variable  whose  conditional  probability 
density  function  can  be  approximated  as  a  weighted  sum  of 


Dirac  delta  functions.  The  weights  that  are  associated  with 


each  of  the  a— priori  models  are  dependent  on  the  fit  of  a 


given  model  to  the  current  real-world  condition  and  are 
updated  from  modeling  errors  generated  for  each  candidate 
model . 

As  is  discussed  earlier,  the  parameter  vector,  G(kT), 
is  composed  of  the  coefficients  of  an  autoregressive  differ¬ 
ence  equation.  Each  flight  condition  has  a  parameter  vector 
associated  with  it  and  hence  becomes  the  source  of  the  a- 
priori  information  required  by  the  multiple  model  algorithm. 
The  candidate  model  parameter  vectors  are  designated  as 

G.  (kT)  i  =  1 ,2 . N  (4-1  ) 

l 

where  N  is  the  number  of  candidate  models  being  considered. 

One  method  of  estimating  G(kT)  is  to  form  a  weighted 
sura  of  a-priori  models  as  given  by  the  equation 

N 

9 ( kT )  =  L  a. (kT)©. (kT)  (4-2) 

icl1  _1 

where  a.(kT)  is  a  weighting  factor  associated  with  a  partic¬ 
ular  candidate  model.  The  conditions  on  the  weighting  fac¬ 
tors  are  given  by 


0  < 


a. (kT) 

i 


<  1 


(4-3) 
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and 


N 

£  ai(kT)  =  1  (4-4) 

From  the  above  equations,  estimating  the  parameter 

vector  0(kT)  becomes  a  matter  of  properly  choosing  the 

weighting  factors,  a.(kT),  associated  with  each  candidate 

1 

model,  which  in  turn  are  dependent  on  the  fit  of  the  asso- 


c 


As  shown  in  this  thesis,  the  multiple  model  algorithm  as 
given  by  Wittenmark  (24)  utilizes  the  variance  of  the  pre¬ 


diction  error  as  the  method  of  determining  the  best  model 

fit. 

The  multiple  model  algorithm  (9,18,19,24)  provides  a 
method  of  blending  a-priori  information  with  parameter  esti¬ 
mation  to  formulate  a  weighting  factor  for  each  candidate 
model.  As  stated  earlier,  this  approach  views  0(kT)  as  a 
random  vector  that  has  a  conditional  probability  density 
function  approximated  as  a  weighted  sum  of  Dirac  delta 
functions.  The  quantities  that  affect  the  shape  of  the 
approximated  probability  density  function  are  the  mean  vec¬ 
tor,  covariance  and  the  weighting  factors  associated  with 
each  of  the  individual  models  of  the  multiple  model  algo- 
r i thm . 


The  mean  of  each  individual  distribution  is  the  a- 


priori  model  parameter  values  themselves.  The  covariance  of 


each  distribution  is  the  corresponding  parameter  covariance 


matrix,  P  (kT),  which,  as  shown  later,  does  not  need  to  be 
“l 


directly  estimated.  The  weighting  factors,  a  ( kT  )  ,  are  se¬ 


lected  via  a  Bayesian  approach  and  are  calculated  based  on 


the  assumption  that  the  residuals  of  the  best  fitting  model 


£ 


to  the  true  model  are  zero  mean  white  Gaussian.  This  obser¬ 


vation  can  be  used  to  test  which  model  best  fits  the  input/ 
output  data.  Assuming  an  a-priori  distribution,  p(0j£(O)), 
for  each  model  (if  there  is  no  a-priori  knowledge  the  as¬ 
sumption  can  be  made  that  all  models  are  equally  probable, 

i.e.  p(0  |  £(0))=1/N  ,  a  conditional  probability  can  be 

i 

calculated  and  updated  following  a  measurement  by  the  equa¬ 
tions  as  given  by  Kittenmark  (24) 

p(0|y(kT))  =  p(0  Iv(kT)  ,y(k-l )T) 

=  p( x(kT)  I  0,z< (k-1 )T) )p(© ly( (k-1 )T) ) 
p(}MkT)  |  £( (k-1 ) T )  ) 

where 

p(y(kT)  |  0,y(  (k-1  )T)  )  =  (  2  7T  A  (  kT  )  )  -  1  /2  *E  (4-7) 

and 

T  - 1 

E  =  exp { - . 5e ,  ( kT  )  ( A .  ( kT  )  )  e.(kT)) 

-i  “le  l 

Therefore  the  weighting  coefficients,  a  (kT) ,  take  the  form 
of 

-1/2  T  - 1 

a  (kT)=Ca. ( (k-1 )T) [A.  (kT)]  exp{-.5e.  ( kT ) ( A  )  e.(kT)) 

l  l  “ie  l  le  l  . 


(  4-5  ) 

(4-6  ) 


c 


where  C  is  a  normalization  factor,  such  that 


E  {a.(kT)}=l 


(4-9) 


6(kT)  =  E  (a. (kT)9. ( kT ) ) 
i=l  1 


(4-10) 


and  e  is  the  residual  generated  by  the  ith  model  and  A 

i  ~ie 

is  the  prediction  error  covariance.  Multiple  model  algo¬ 
rithm  weighting  coefficient  update  equations  similar  to 
those  presented  above  are  given  by  Maybeck  (9)  and  Andersson 


4.2.1 _ Prediction  Error.  To  define  the  error  associated 

with  a  specific  candidate  model  (the  ith  model),  the  predic¬ 
tion  error  or  residual  is  defined  from  Equation  (2-11)  as 


e.  (kT)  =  y(kT)  -  T  (kT)6  (kT) 

l  l 


(4-11) 


The  parameter  error  is  defined  as 


A9.  (kT)  =  0.  (kT)  -  0  (  kT  ) 


( 4-12 ) 


where  ©(kT)  is  the  best  fitting  model  parameter  vector.  By 
incorporating  Equation  (4-12)  into  Equation  (4-11),  the 
prediction  error  becomes 


e.(kT)  =  y ( kT )  -  TT(kT)6(kT>  -  T T ( kT ) Ae . < kT >  (4-13) 

l  i 


'  *-  * 


The  first  two  terms  on  the  right  hand  side  of  Equation 
(4-13)  represent  the  error  introduced  by  actuator  and  sensor 
noise  as  well  as  the  error  due  to  linear  approximation 
errors.  The  third  term  is  an  additional  modeling  error  term 
induced  by  the  model  not  being  the  best  fitting  model. 

4.2.2 _ Prediction  Error  Covariance.  As  can  be  seen  in 

Equation  (4-8),  the  multiple  model  algorithm  as  given  by 
Wittenraark  (24)  utilizes  the  covariance  of  the  prediction 
error  as  the  method  of  determining  the  best  model  fit. 

Berens  (2)  showed  that  an  estimate  of  the  prediction  error 
covariance  is  given  by 

k 

A.  <kT)  ~  (  1/k  > (  E  e  .  ( i i >  *e .  1  (  ii  )  )  (4-14) 

le  . .  i  l 

ii  =  I 

Equation  (4-14)  is  an  ergodic  approximation  to  the  true 

ensemble  average  and  assumes  e.  is  zero-mean. 

'  l 

The  prediction  error  covariance  of  systems  with  gene¬ 
ralized  inputs  is  slowly  time  varying.  Old  values  of 

A  (kT)  are  discarded  exponentially  by  replacing  Equation 
i  e 

(4-14)  with  a  recursive  fading  memory  filter  (9), 

T 

A.  (kT)  —  j  A.  (  ( k  -  1 ) T )  +  (1  - j  )  e .  ( kT ) e .  (kT)  (4-15) 


where 


0  <  j  <  1 

Small  j  corresponds  to  short  memory  and  j  large  corresponds 
to  long  memory. 

The  performance  of  Equation  (4-15)  is  restricted  by  the 

fact  that  the  prediction  error  covariance,  A.  (kT),  is  as- 

le 

sumed  to  vary  slowly  with  time.  When  the  time  varying 

assumption  is  violated,  the  performance  of  Equation  (4-15) 

can  be  enhanced  by  low  pass  filtering  of  the  input,  u(kT), 

and  the  output,  jMkT)  (discussed  in  Section  4.3).  From 

Equations  (4-11)  and  (4-14),  it  can  be  seen  that  A.  ( kT )  is 

"  xe 

a  function  of  T(kT).  It  is  shown  earlier  that  T(kT)  is 
formed  entirely  of  past  values  of  u(kT)  and  y(kT).  There¬ 
fore,  by  reducing  the  high  frequency  variation  in  u(kT)  and 
y(kT)  via  low  pass  filtering,  high  frequency  variation  in 

A  ( kT !  is  also  reduced, 
l  e 

4_._3  Mu  1 1  i  pi  e  _Mode  1_  Parameter  Estimator 

The  multiple  model  parameter  estimator  which  incor¬ 
porates  a-priori  information  for  use  in  in-flight  simulation 
is  shown  in  Figure  4-1  (9).  As  discussed  earlier,  the  a- 


referred  to  as  adaptive  or  composite)  estimate  of  the  best 


fitting  model  parameters.  The  primary  estimates  are  then 


vmv 


used  to  update  the  control  law  gains. 

The  secondary  parameter  estimation  task  can  be  accom¬ 
plished  by  incorporating  one  of  several  different  recursive 
estimation  techniques  such  as  recursive  least  squares  or 
Hagglund's  algorithm.  The  number  of  parameters  to  be  esti¬ 
mated  by  the  recursive  parameter  estimation  routine  is  also 
a  design  option.  The  secondary  estimator  may  be  chosen  to 
estimate  all,  some  (as  in  Pineiro’s  simulation),  or  none  of 
the  parameters. 

A  number  of  problems  become  apparent  if  the  secondary 
parameter  estimators  are  allowed  to  update  all  the  parame¬ 
ters  recursively  at  every  sampling  period.  As  the  number  of 
inputs,  outputs,  and  the  order  of  the  autoregressive  dif¬ 
ference  equation  increases,  the  computational  effort  may 
become  impractical  due  to  long  convergence  and  computation 
time.  Also,  as  the  number  of  parameters  to  be  updated 
increases,  the  system  requires  more  excitation  for  parameter 
convergence  (17). 

To  reduce  the  computational  effort  afforded  by  full 
scale  secondary  parameter  estimation  and  to  decrease  the 
convergence  time,  the  secondary  estimators  may  be  allowed  to 
estimate  part  of  the  parameter  vector  while  the  remaining 


parameters  are  assumed  to  be  fixed  at  a  nominal  flight  con¬ 
dition.  At  times,  this  assumption  may  not  be  satisfactory 
and  in  fact  it  may  preclude  convergence  to  the  proper  model. 

Secondary  parameter  estimation  can  be  avoided  altoge¬ 
ther  by  leaving  all  of  the  parameters  in  each  secondary 
estimator  fixed  at  a  nominal  flight  condition.  This  method 
requires  relatively  little  computational  effort  and  provides 
good  estimates  of  the  parameter  vector,  especially  when  a 
robust  control  law  is  being  implemented.  Fixed  secondary 
parameter  estimation  is  the  method  incorporated  in  this 
thes i s . 

Three  filters  are  incorporated  into  the  system  shown  in 
Figure  1-1  to  smooth  variations  in  the  parameter  estimates. 
The  first  is  a  digital  band  pass  filter  which  filters  ydk) 
and  u(k).  The  low  frequency  components  of  the  input-output 
signals  must  be  removed  to  reduce  parameter  estimate  bias, 
while  high  frequency  signal  components  are  removed  to  smooth 
input  excitation  (27).  The  filter  is  implemented  as  a 
sixth-order  butterworth  digital  bandpass  filter  (2,28).  A 
second  non-linear  filter  is  added  to  the  weighting  coeffi¬ 
cients  within  the  multiple  model  algorithm.  The  purpose  of 
this  filter  is  to  smooth  sudden  changes  in  the  model  proba- 


bilities  that  could  result  in  destabilizing  rapid  changes  in 


the  primary  parameter  estimates.  Therefore,  the  rate  at 
which  the  weighting  coefficients  can  change  in  a  given 
sampling  period  is  controlled  by  a  rate  limiting  filter  (2). 
The  third  filter  is  a  low  pass  filter  that  filters  the 
primary  parameter  estimates  before  they  are  implemented  by 
the  control  law. 

4 . 4  Summary 

This  chapter  presents  the  multiple  model  estimation 
algortihm  that  was  used  in  this  thesis.  The  equations  used 
to  evaluate  the  weighting  coefficients  for  the  correctness 
of  candidate  models  are  presented  as  well  as  the  overall 
implementation  of  the  multiple  model  adaptive  estimator. 

The  next  chapter  presents  important  factors  to  be  considered 
when  selecting  models  for  the  multiple  mod^l  estimation 
algorithm  in  the  context  of  achieving  uesireo  tracking  per¬ 
formance  over  a  desired  parameter  space. 


5.1  Introduction 

Because  the  aircraft  can  maneuver  through  a  continuum 
of  different  flight  conditions,  an  infinite  number  of  state 
space  system  models  of  the  aircraft  are  possible.  The  large 
number  of  system  models  that  would  be  required  to  character¬ 
ize  all  sets  of  potential  parameters  could  not  be  implemen¬ 
ted  in  a  realistic  control  system  due  to  the  computational 
burden  that  would  be  imposed  as  well  as  having  a  limited 
amount  of  computer  memory  available.  To  alleviate  this 
situation,  the  multiple  model  algorithm  with  a  discrete  set 
of  system  models  to  represent  the  region  of  selected  parame¬ 
ter  space  can  be  implemented.  (The  parameter  space  for  this 
thesis  is  defined  in  Chapter  2  and  is  shown  in  Figure  2-4.) 
Selecting  a  discrete  number  of  models  to  represent  a  con¬ 
tinuous  parameter  space  means  that  the  actual  values  of  the 
parameters  do  not  correspond  exactly  to  one  of  the  models 
but  are  close  to  one.  Therefore,  the  models  must  be  ’close’ 
enough  to  one  another  to  ensure  that  the  selected  model  is 
indeed  representative  of  the  actual  flight  condition. 

However,  the  system  models  cannot  be  ’too  close'  or 


difficulties  occur  with  the  multiple  model  algorithm.  If  the 
system  models  are  such  that  the  difference  between  measure¬ 
ment  data  and  the  model  predicted  values  is  not  significant 
for  more  than  one  model,  then  selection  of  the  ’proper’ 
model  cannot  take  place.  As  shown  in  Chapter  4,  the  multi¬ 
ple  model  algorithm  generates  a  residual  for  each  of  the 
system  models  which  is  then  passed  on  to  be  processed  by  the 
hypothesis  conditional  probability  computation  (see  Figure 
4-1).  The  algorithm  performance  depends  on  the  fact  that 
system  models  that  are  not  representative  of  the  current 
flight  condition  have  large  residuals,  while  those  models 
that  are  near  the  true  flight  condition  have  small  resid¬ 
uals.  Therefore,  it  can  be  seen  that  a  careful  choice  of 
system  models  must  be  made.  The  selected  models  must  be 
close,  but  not  too  close.  This  chapter  presents  the  cri¬ 
teria  upon  which  model  placement  was  based  to  allow  the 
multiple  model  algorithm  to  select  the  model  nearest  the 
actual  flight  condition. 

5 . 2 _ Performance  Boundaries 

The  performance  of  an  actual  physical  system  can  only 
be  approximated  by  an  appropriate  mathematical  model.  The 
degree  that  a  mathematical  model  represents  an  actual  system 


to  depends  on  the  model  uncertainties  that  are  present  at  a 


given  instant  in  time.  These  model  uncertainties  can  arise 
in  several  ways,  including  parameter  variations  or  incor¬ 
rectly  modeled  dynamics.  Both  stability  and  satisfactory 
performance  are  of  prime  concern  in  the  face  of  these  model¬ 
ling  uncertainties  (29). 

In  designing  a  feedback  control  system,  a  model  (re¬ 
ferred  to  as  the  nominal  model)  must  be  selected  from  the 
set  of  available  models  to  represent  the  plant’s  behavior. 
When  a  design  has  been  performed  with  a  specific  nominal 
model,  it  is  obvious  that  modelling  errors  are  present  when 
there  is  any  deviation  of  the  plant  model  from  the  nominal 
design  model.  A  system  is  labeled  as  robust  with  respect  to 
these  modelling  errors  when  the  feedback  control  system 
remains  stable  when  the  nominal  plant  model  is  replaced  by 
another  plant  without  altering  the  controller  (21). 

Banda  and  Ridgely  (29)  point  out  that  there  are  two 
types  of  robustness;  stability  robustness  and  performance 
robustness.  Stability  robustness  occurs  when  the  closed- 
loop  system  remains  stable  in  the  presence  of  plant  varia¬ 
tions.  Performance  robustness  is  present  when  the  perfor¬ 


mance  of  the  closed-loop  system  is  acceptable  in  the  face  of 


such  variations.  This  chapter  is  primarily  concerned  with 
the  issue  of  performance  robustness. 

Obviously,  if  a  single  control  system  possesses  enough 
performance  robustness,  the  entire  flight  envelope  shown  in 
Figure  2-4  could  be  approximated  with  a  single  nominal 
model.  In  any  case,  a  given  nominal  model  has  a  region  of 
performance  robustness  that  is  associated  with  it  as  shown 
in  Figure  5-1.  This  performance  robustness  region  indicates 
that  satisfactory  system  performance  is  achieved  when  the 
nominal  plant  is  replaced  by  any  plant  model  within  the 
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Figure  5-1.  Region  of  Performance  Robustness 


region  encompassed  by  the  performance  boundary  (see  Figure 


5-1  ) 


Recall  that  an  aircraft  state  space  model  can  be  repre¬ 


sented  by  a  autoregressive  difference  equation 


y(kT)  =  T  <kT)©(kT)  +  e(kT) 


(  5-1  ) 


where  6  contains  the  parameters  of  the  difference  equation, 


In  Chapter  3,  it  is  shown  that  the  control  law  gains  are 


given  by 


K :  =  H~  1  (  T  )£ 


(  5-2  ) 


— 2  =P*1 


(5-3) 


where  the  step  response  matrix  is 


H(T)  =  B, 


(  5-4  ) 


and  can  be  identified  as  a  subset  of  the  parameter  vector 


0.  Each  flight  condititon  has  a  unique  parameter  vector 


associated  with  it  and  hence  a  unique  step  response  matrix, 


Therefore,  it  can  be,  seen  from  Equations  (5-2)  and  (5-3), 


that  each  model  has  a  specific  set  of  gains  associated  with 


it  (assuming  the  tuning  matrix  £  and  the  parameter  p  remain 


cons  tant  )  . 


i 


‘i 
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In  determining  a  performance  boundary  for  a  specific 
nominal  model,  the  control  law  gains  remain  fixed  at  the 
values  that  are  determined  by  the  nominal  model’s  parameter 
vector.  By  replacing  the  nominal  plant  model  with  other 
models  and  driving  the  closed-loop  system  with  the  commanded 
inputs  outlined  in  Chapter  2,  the  performance  boundary  can 
be  established  for  a  given  performance  criteria.  Hence, 
plant  models  that  are  outside  the  performance  boundary  fail 
to  meet  the  desired  performance  criteria,  while  those  models 
inside  the  boundary  yield  acceptable  closed-loop  perfor¬ 


mance  . 


The  performance  criteria  used  in  this  thesis  are  the 
criteria  found  in  the  VISTA  statement-of-work  and  are  also 
the  criteria  used  by  Pineiro  and  Berens.  The  criteria  for 
tracking  performance  require  that  the  average  response  error 
absolute  value  should  not  exceed  the  average  absolute  value 
of  ten  percent  of  the  reference  signal.  The  time-average 
response  error  absolute  value  is  given  by 


1/t  j  (cmd  input  -  output)  dt 

0.01 


(  5-5  ) 


c 


where 


t  =  time 


cmd  input  =  commanded  input 


The  time-average  value  of  ten  percent  of  the  reference 


signal  is 


1/t  j  0.1  *  cmd  input  dt 

0.01 


(  5-6  ) 


For  satisfactory  closed-loop  tracking  response  to  be 


accomplished  over  the  flight  envelope  of  interest,  the  per¬ 


formance  boundaries  of  the  required  number  of  models  must 


overlap  as  shown  in  Figure  5-2.  This  requirement  is  due  to 


the  fact  that  if  the  operating  point  lies  outside  the  per¬ 


formance  boundaries,  the  multiple  model  algorithm  assigns 


probabilites  to  the  respective  models  that  result  in  control 


law  gains  that  cause  the  closed  loop  system  to  not  meet  the 


performance  criteria.  From  Figure  5-2,  it  can  be  seen  that 


there  is  a  minimum  number  of  models  that  will  result  in  the 


performance  boundaries  totally  covering  the  desired  flight 


envelope.  This  thesis  addresses  the  formation  of  perfor¬ 


mance  boundaries,  the  number  of  models  required  to  yield 
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Figure  5-2.  Overlapping  of  Performance  Boundaries 

acceptable  tracking  performance  over  the  flight  envelope, 
and  the  amount  of  overlapping  of  the  performance  boundaries 
that  the  multiple  model  algorithm  requires. 

5.3  Sensor  Noise  Effects 

By  simple  algebraic  manipulation  of  the  autoregressive 
difference  equation,  it  was  shown  in  Chapter  4  that  the 
prediction  error  is  defined  as 

e  ( kT )  =  y ( kT )  -  TT(kT)9(kT)  (5-7) 
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Therefore,  the  prediction  error  covariance  becomes 

T 

E(e(kT)e  ( kT ) )  =  A  (5-8) 

~ e 

If  the  operating  point  of  a  fixed  gain  system  (de¬ 
scribed  above)  is  the  nominal  plant  model,  then  the  predic¬ 
tion  error  variance  (assuming  no  sensor  noise)  is  due  to  the 
modeling  error  introduced  when  representing  a  continuous 
system  with  an  autoregressive  difference  equation.  When 
sensor  noise  is  present,  then  the  prediction  error  variance 
is  composed  of  the  linear  approximation  errors  just  men¬ 
tioned  and  the  error  introduced  by  addition  of  the  sensor 
noise. 

As  shown  in  Chapter  4,  the  multiple  model  algorithm 
assigns  a  conditional  probability  to  each  candidate  model 
based  on  the  calculated  value  of  each  model’s  prediction 
error  covariance.  When  selecting  a  model  as  the  best  model 
to  represent  the  current  plant  dynamics  from  a  set  of  candi¬ 
date  models,  the  prediction  error  covariance  is  composed  of 
three  error  terms.  These  error  terms  are  the  errors  intro¬ 
duced  by  sensor  noise  and  linear  approximation  (as  for  the 
fixed  gain  case),  as  well  as  a  modelling  error  term  induced 
by  the  model  not  being  the  best  fitting  model.  Therefore, 
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for  the  multiple  model  case,  the  prediction  error  covariance 
becomes 

A.  =  A  +  A  +  A  ( 5-9  ) 

le  e  n  1m 

where 

A.  is  the  total  error  covariance  associated 

_ie 

with  the  ith  candidate  model 

A  is  the  error  covariance  introduced 
e 

by  linear  approximation 

A  is  the  error  covariance  introduced 

by  the  addition  of  sensor  noise 

A  is  the  error  covariance  introduced 
l  m 

by  a  model  not  being  the  best  fitting  model 

It  can  be  seen  from  Equation  (5-9)  that,  in  order  for 
each  candidate  model  to  have  a  distinct  prediction  error 
covariance,  the  third  term  on  the  right  hand  side  of  the 
equal  sign  must  be  the  dominant  term.  If  the  error  asso¬ 
ciated  with  the  sensor  noise  becomes  the  dominant  term,  it 
will  mask  the  modelling  error  for  each  of  the  candidate 
models.  In  essence,  a  large  sensor  noise  causes  the  •• 
siduals  generated  for  each  candidate  model  to  b>  t?  . 
and  not  distinct.  Therefore,  proper  model 
occur . 

The  above  discussion  results  ;  r. 


6 


ty  which  must  be  satisfied  in  order  for  the  multiple  model 
algorithm  to  be  able  to  distinguiush  one  model  from  another. 


A  +  A.  >  >  A 
~e  "im  “h 


(5-10) 


Assuming  that  the  linear  approximation  error  term  is  small, 
equation  (5-10)  becomes 


A.  >>  A 
"im  ~ n 


(5-11) 


Equation  (5-11)  shows  that,  for  a  given  level  of  sensor 
noise,  there  is  a  minimum  spacing  between  models  that  is 
necessary  for  the  multiple  model  algorithm  to  determine  the 
proper  model.  This  is  shown  graphically  in  Figure  5-3. 
Figure  5-3  shows  that  if  a  second  model  is  placed  in  the 
region  where  the  sensor  noise  would  mask  the  residuals,  then 
proper  model  selection  cannot  occur.  Therefore,  there  is  a 
region  that  is  not  acceptable  for  placement  of  models  subse¬ 
quent  to  the  placement  of  an  initial  model  (determined  by 
sensor  noise  level).  On  the  other  hand,  when  models  are 
placed  such  that  the  modelling  errors  dominate  the  sensor 
noise,  proper  model  selection  can  be  made.  Additionally, 
Equation  (5-11)  shows  that  the  farther  apart  the  candidate 


moat's  are  spaced,  the  more  noise  resistant  the  multiple 
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1.  Region  where  get  substantial  model¬ 
ing  error.  Acceptable  region  for  place¬ 
ment  of  subsequent  model. 

2.  Region  where  sensor  noise  would  mask 
residuals.  Unacceptable  region  for  place¬ 
ment  of  subsequent  model . 


Mach  Number 


Figure  5-3.  Noise  Effects  on  Model  Selection 
model  algorithm  becomes.  This  thesis  addresses  the  issues 
just  discussed  by  investigating  how  model  selection  for  the 
multiple  model  algorithm  is  affected  by  the  addition  of 
sensor  noise. 


5 ■ 4  Summary 

This  chapter  presents  the  factors  considered  in  this 
thesis  for  selecting  models  for  the  multiple  model  algo¬ 


rithm.  Specifically,  the  concept  of  establishing  perfor- 


mance  boundaries  using  fixed  gain  simulations  is  discussed, 
showing  how  overlapping  boundaries  can  be  used  to  establish 
desired  tracking  performance  over  a  given  set  of  flight 
conditions.  Sensor  noise  effects  are  also  shown  to  be  a 
factor  to  be  considered  when  selecting  a  discrete  number  of 
aircraft  models  for  the  multiple  model  algorithm. 
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Simulation  Setup  and  Results 


6.1  Simulation  Setup 

The  parameter  adaptive  simulation  software  used  in  this 
thesis  was  developed  by  Pineiro  and  Berens.  Figure  6-1 
shows  the  top  level  description  of  the  parameter  adaptive 
system  simulation  package  conducted  with  the  aid  of  MATRIX 

x 

control  design  software  and  its  System  Build  simulation 
tool.  The  block  in  Figure  6-1  that  is  labeled  ADAPT2  imple¬ 
ments  the  multiple  model  parameter  estimation  algorithm  and 
also  performs  the  design  calculations  of  updating  the  con¬ 
trol  law  gains  and  •  The  block  labled  CTRL  accepts  the 
updated  control  law  gains  and  implements  the  control  law 
equations  that  are  presented  in  Equations  (3-47)  and  (3-48) 
with  a  sampling  time,  T,  of  0.01  second.  The  block  labled 
A/C  implements  the  longitudinal  linearized  dynamics  of  the 
AFTI/F-16  as  presented  in  Chapter  2.  The  A/C  block  also 
implements  actuator  dynamics  and  includes  nonlinearities 
such  as  rate  and  position  limits.  The  block  labled  SEN2 
allows  sensor  noise  to  be  added  to  the  longitudinal  states 
as  desired.  For  a  more  detailed  presentation  of  each  of  the 
blocks  in  Figure  6-1  consult  Beren’s  thesis  (2). 
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Figure  6-1.  Top  Level  Description  of  Parameter  Adaptive  Simulation  Package 
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The  reference  tracking  signals  used  are  those  that  were 
implemented  by  Pineiro  and  Berens.  The  inputs  were  obtained 
from  real-time,  nonlinear  simulation  with  elevators  and 
flaperons  as  control  surfaces  and  with  flight  path  angle  and 
pitch  rate  as  outputs.  The  reference  tracking  signals  are 
shown  on  Figures  6-2  and  6-3. 

In  preparation  to  run  the  simulations,  several  data 
files  were  created  to  input  specific  flight  condition  data 
along  with  the  necessary  control  surface  deflection  limits 
into  the  parameter  adaptive  simulation  software.  These  data 
files  consist  of  the  A  and  B  matrices  of  the  state  space 
representation  of  a  flight  condition  as  presented  in  Appen¬ 
dix  A,  along  with  the  elevator  and  flaperon  position  limits 
as  listed  in  Appendix  C.  This  flight  condition  data  was 
then  input  using  the  MATRIX^  System  Build  options  into  the 
block  labled  A/C  as  shown  on  Figure  6-1.  In  addition,  sever¬ 
al  parameter  vectors  were  created  from  the  state  space 
models  using  MATRIX^.  These  parameter  vectors  were  input  as 
needed  into  the  ADAPT2  block  (see  Figure  6-1)  as  the  fixed 
secondary  parameter  estimator  models  of  the  multiple  model 
parameter  estimator  algorithm  (see  Figure  4-2). 
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Figure  6-2.  Pitch  Rate  Input 
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Figure  6-3.  Flight  Path  Angle  Input 


6.2. 1.1  Performance  Boundary  Evaluation.  In  preparation 
for  determining  the  number  of  models  and  the  spacing  of 
those  models  to  obtain  satisfactory  tracking  performance 
over  the  desired  flight  envelope,  fixed  gain  simulations 
were  run  to  determine  the  performance  boundaries  of  a  single 
model.  Recall  that,  as  discussed  in  Chapter  5,  by  using  the 
performance  criteria  of  Equations  (5-5)  and  (5-6)  perfor¬ 
mance  boundaries  for  a  given  model  can  be  obtained  by  re¬ 
placing  the  nominal  state  space  aircraft  model  with  another 
aircraft  model  and  running  subsequent  simulations.  If  the 
performance  criteria  are  satisfied,  then  the  non-nominal 
aircraft  model  yields  satisfactory  tracking  performance  with 
the  nominal  models  control  law  gains.  This  satisfactory 
performance  dictates  that  the  off-nominal  model  is  within 
the  performance  boundary  for  the  nominal  model.  Failure  to 
meet  the  performance  criteria  indicates  that  the  off-nominal 
model  is  outside  of  the  nominal  model  performance  boundary. 

The  first  point  chosen  as  a  nominal  aircraft  model  was 
at  the  flight  condition  of  18,000  ft  altitude  and  0.45  Mach. 
The  parameter  vector  (Appendix  D)  was  calculated  and  imple- 


merited  into  parameter  adaptive  system  for  the  purpose  of 
defining  the  control  law  gains  via  the  step  response  matrix 
for  fixed  gain  simulations.  The  first  simulation  was  run 
with  the  nominal  aircraft  plant  and  the  nominal  gains.  The 
response  of  this  simulation  is  shown  on  Figures  6-4  and  6-5. 
Figure  6-4  plots  the  pitch  rate  input  versus  the  output 
response  while  Figure  6-5  shows  a  plot  of  the  performance 
index  for  this  simulation.  Figures  6-6  through  6-9  present 
the  elevator  and  flaperon  positions  and  rates  for  this 
simulation.  From  the  data  presented  for  this  simulation,  it 
can  be  seen  that  the  performance  criteria  are  well  satisfied 
and  that  tracking  response  of  the  reference  signal  is  very 
good.  It  was  noticed  early  in  the  simulations  that  the 
pitch  rate  response  was  the  critical  response  for  the  com¬ 
manded  inputs.  The  pitch  rate  response  failed  the  perfor¬ 
mance  criteria  well  before  the  flight  path  angle  response 
for  all  simulations  run.  For  this  reason,  only  the  pitch 
rate  responses  are  presented.  Subsequent  simulations  were 
made  by  replacing  the  nominal  aircraft  model  with  aircraft 
models  at  the  same  altitude  and  decreased/increased  Mach 
number.  Mach  number  was  decreased/increased  from  the  nomi¬ 
nal  point  in  increments  of  0.05  until  the  resulting  data 
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Figure  6-4.  Pitch  Rate  Response 
Nominal  Model:  18K  0 . 45M/Operating  Point:  18K  0.45M 


Figure  6-5.  Pitch  Rate  Performance  Criterion 
Nominal  Model  18K  0 . 45M/Operating  Point:  18K  0.45M 
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Figure  6-6.  Flaperon  Position 
Nominal  Model:  18K  0 . 4 5M/0perating  Point:  18K  0.45M 
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Figure  6- 

Nominal  Model:  18K 

7.  Elevator  Position 

0 . 45M/Operating  Point:  18K  0.45M 

showed  that  the  performance  criteria  had  been  violated. 
Figures  6-10  through  6-15  present  the  data  when  the  aircraft 
model  is  set  at  18,000  ft  and  0.8  Mach.  From  the  data,  it 
can  be  seen  that  the  performance  criteria  is  violated  and 
hence  the  model  is  outside  the  performance  boundary  for  the 
nominal  model  of  18,000  ft,  0.45  Mach. 

Following  the  simulations  at  18,000  ft,  the  control  law 
gains  for  the  nominal  plant  were  used  with  aircraft  plant 
models  at  altitudes  of  10,000,  14,000,  22,000,  26,000, 
30,000,  34,000  and  38,000  ft  with  Mach  number  varying  from 
0.25  to  0.9  in  increments  of  .05.  The  flight  conditions 
that  failed  to  meet  the  performance  criteria  are  found  in 
Table  6-1.  Using  the  data  presented  in  Table  6-1,  the 
performance  boundary  for  the  18,000  ft,  0.45  Mach  nominal 
point  can  be  obtained  and  is  shown  on  Figure  6-16.  From 
Figure  6-16,  it  can  be  seen  that  the  performance  boundary 
encompasses  the  dynamic  pressure  line  on  which  the  nominal 
plant  resides.  Figure  6-16  also  shows  that  the  upper  and 
lower  limits  of  the  performance  boundary  are  contingent  on 
the  differential  of  dynamic  pressure  between  the  nominal 
model  and  the  actual  flight  condition  model. 


Figure  6-10.  Pitch  Rate  Response 
Nominal  Model:  18K  0 . 4 5M/0perating  Point:  18K  0 . 8M 
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Figure  6-11.  Pitch  Rate  Performance  Criterion 
Nominal  Model:  18K  0 . 45M/Operating  Point:  18K  0 . 8M 
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Figure  6-12.  Flaperon  Position 
Nominal  Model:  18K  0 . 45M/Operating  Point:  18K  0 . 8M 


Figure  6-13.  Elevator  Position 
Nominal  Model:  18K  0 . 45M/Operating  Point:  18K  0.8M 
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Table  6-1 

Flight  Conditions  Failing  Performance  Criteria 
Nominal  Flight  Condition  =  18,000  feet  0.45  Mach 


Flight  Conditions  of 

Failure 

Altitude 

(Kft)  Mach  Number 

Dynamic  Pressure 

10 

0.25 

63.68 

10 

0.70 

499.27 

14 

0.30 

78.35 

14 

0.75 

489.69 

18 

0.30 

66.62 

18 

0.80 

473.75 

22 

0.35 

76.71 

22 

0.85 

452.45 

26 

0.40 

84.30 

26 

0.90 

426.79 

30 

0.45 

89.26 

34 

0.50 

91.61 

38 

0.55 

91.61 

Due  to  the  fact  that  the  performance  boundary  above 
seemed  to  be  dependent  on  dynamic  pressure,  the  next  flight 
condition  was  chosen  so  that  the  dynamic  pressure  was  ap¬ 
proximately  the  same  as  the  18,000  ft  0.45  Mach  point.  This 
led  to  choosing  the  next  flight  condition  at  38,000  ft,  0.7 
Mach.  Flight  conditions  that  failed  to  meet  the  performance 
criteria  for  this  nominal  flight  condition  are  given  in 
Table  6-2  and  the  resulting  performance  boundary  is  presen¬ 
ted  on  Figure  6-17.  Comparing  the  performance  boundaries  of 


Table  6-2 

Flight  Conditions  Failing  Performance  Criteria 
Nominal  Flight  Condition  =  38,000  feet  0.7  Mach 


Flight  Conditions 

of  Failure 

Altitude 

(Kft)  Mach 

Dynamic  Pressure 

10 

0.25 

63.68 

10 

0.70 

499.27 

14 

0.30 

78.35 

14 

0.75 

489.69 

18 

0.35 

90.68 

18 

0.80 

473.75 

22 

0.35 

76.71 

22 

0.85 

452.45 

26 

0.40 

84.30 

26 

0.90 

426.79 

30 

0.45 

89.26 

34 

0.50 

91.61 

38 

0.55 

91.61 

Figures  6-16  and  6-17  reveals  that  they  are  essentially  the 
same.  This  indicates  that,  from  a  performance  criteria 
standpoint,  nominal  models  with  approximately  the  same  dyna¬ 
mic  pressure  yield  the  same  performance  boundaries.  Several 
other  performance  boundaries  were  obtained  for  nominal  mo¬ 
dels  with  similar  dynamic  pressures.  These  flight  condi¬ 
tions  along  with  the  associated  dynamic  pressures  are  listed 
in  Table  6-3.  The  result  of  obtaining  performance  bounda¬ 
ries  for  the  flight  conditions  in  Table  6-3  showed  that,  for 
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Figure  6-17.  Performance  Boundary  for 
Nominal  Model  at  38K  0.7M 
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Table  6-3. 

Flight  Conditions  of  Nominal  Models  With 
Similar  Dynamic  Pressures 


Case 

Altitude  (Kft) 

Mach 

Dynamic  Pressure 

1 

22 

0.65 

264.58 

34 

0.85 

264.75 

2 

14 

0.70 

426.57 

26 

0.90 

426.79 

3 

26 

0.50 

131.73 

34 

0.60 

131.92 

4 

26 

0.75 

296.38 

34 

0.90 

296.82 

5 

30 

0.50 

110.19 

34 

0.55 

110.85 

6 

34 

0.50 

91.61 

38 

0.55 

91.61 

7 

34 

0.75 

206 . 12 

10 

0.45 

206.33 

8 

26 

0.60 

189.69 

22 

0.55 

189.43 

each  of  the  presented  cases,  models  with  similar  dynamic 
pressures  yield  the  same  set  of  performance  boundaries. 
Another  point  that  became  evident  was  that,  when  running  a 
fixed  gain  simulation,  if  the  dynamic  pressure  of  the  air¬ 
craft  model  differed  signi f icantly  from  that  of  the  nominal 
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model  the  performance  criteria  would  be  violated. 

The  fact  that  aircraft  models  with  similar  dynamic 
pressures  yield  equivalent  performance  boundaries  is  also 
indicated  by  examining  the  equations  for  the  dimensionalized 
derivatives  of  the  aircraft  state  space  model.  Recall  that 
the  numbers  generated  by  implementing  the  dimensionalized 
derivative  equations  are  inserted  directly  into  the  A  matrix 
of  the  aircraft  state  space  model.  By  examining  the  dimen¬ 
sionalized  derivative  equations,  it  can  be  seen  that  dynamic 
pressure  is  present  in  all  but  one  of  the  equations.  Not 
only  is  dynamic  pressure  present  in  the  dimensionalized 
derivative  equations,  it  is  also  the  dominant  variable  of 
the  parameters  that  are  not  constant.  Therefore,  dynamic 
pressure  can  be  seen  to  have  a  large  influence  on  the  equiv¬ 
alency  of  one  aircraft  model  with  another. 

Upon  comparing  the  fixed  gain  simulation  data  generated 
for  the  18,000  ft  altitude,  0.45  Mach  and  38,000  ft  alti¬ 
tude,  0.7  Mach  flight  conditions,  an  interesting  observation 
was  made.  Although  the  performance  boundaries  were  equiva¬ 
lent  for  both  flight  conditions,  the  responses  of  the  con¬ 
trol  surfaces  differed  for  equivalent  simulations.  As  an 
example,  consider  the  flight  condition  of  26,000  ft  altitude 
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and  0.45  Mach.  Figures  6-18  through  6-23  show  the  responses 
for  the  case  where  the  nominal  model  was  18,000  ft  altitude, 
0.45  Mach,  and  Figures  6-24  through  6-29  present  data  for 
the  case  where  the  nominal  model  was  38,000  ft  altitude  and 
0.7  Mach.  Comparison  of  the  elevator  and  flaperon  surface 
deflections  for  each  nominal  case  shows  that,  although  the 
general  responses  are  similar,  the  nominal  model  of  38,000 
ft  altitude  and  0.7  Mach  has  noisier  surface  deflection 
responses.  The  elevator  deflection  rates  are  much  larger 
for  the  38,000  ft,  0.7  Mach  case  and  the  flaperon  deflection 
rates  reach  the  rate  limits  more  of  the  time.  Because  the 
control  surfaces  respond  differently  for  the  same  flight 
conditions,  the  equivalence  of  one  nominal  model  with  ano¬ 
ther  is  based  on  more  than  just  dynamic  pressure.  If  dyna¬ 
mic  pressure  were  the  sole  measure  of  model  equivalence,  not 
only  should  the  performance  boundaries  be  the  same  for 
nominal  models  with  the  same  dynamic  pressure,  but  the 
surface  deflections  and  rates  should  be  also.  This  would 
dictate  that  the  eigenvalues  of  the  A  matrices  of  the  nomi¬ 
nal  state  space  models  and  the  control  matrices  be  equal. 

The  eigenvalues  for  the  18,000  ft,  0.45  Mach  and  38,000  ft, 
0.7  Mach  flight  conditions  are  found  in  Table  6-4.  From  the 
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Figure  6-19. 
Nominal  Model: 

Pitch  Rate  Performance  Criterion 

18K  0 . 45M/Operating  Point:26K  0.45M 
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Figure  6-20.  Flaperon  Position 
Nominal  Model:  18K  0 . 45M/Operating  Point:  26K  0.45M 
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Figure  6-21.  Elevator  Position 
Nominal  Model:  18K  0 , 45M/Operating  Point:  26K  0.45M 
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Figure  6-25.  Pitch  Rate  Performance  Criterion 
Nominal  Model:  38K  0 . 7M/0perating  Point:  26K  0.45M 
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Figure  6-26.  Flaperon  Position 
Nominal  Model:  38K  0 . 7M/0perating  Point:  26K  0.45M 


Figure  6-27.  Elevator  Position 
Nominal  Model:  38K  0 . 7M/0perating  Point:  26K  0.45M 
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Table  6-4 


r 

K  Eigenvalues  for  the  A  Matrix  of  State  Space  Model 


Flight  Condition 

Eigenvalues 

38,000  ft 

-0.01254387 

+ 

0.09246732 

i 

0.7  Mach 

-0.01254387 

- 

0.09246732 

i 

Dynamic 

1.11849823 

0.00000000 

i 

Pressure: 148.4 

-1.79023020 

+ 

0.00000000 

i 

18,000  ft 

-0.00386373 

+ 

0.06747448 

i 

0.45  Mach 

-0.00386373 

- 

0.06747448 

i 

Dynamic 

-2.05216607 

+ 

0.00000000 

i 

Pressure: 149 . 9 

1 .11328404 

+ 

0.00000000 

i 
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data  shown  on  Table  6-4,  it  is  apparent  that  the  eigenvalues 
for  the  A  matrices  of  flight  conditions  with  similar  values 
of  dynamic  pressure  are  not  equivalent.  This  conclusion  was 
further  verified  by  examining  the  eigenvalues  of  the  models 
shown  on  Table  6-3.  A  comparison  of  the  B  matrices  (see 
Appendix  A)  for  flight  conditions  with  similar  dynamic  pres¬ 
sures,  shows  they  are  not  equivalent  either. 

The  reason  for  the  differences  in  the  control  surface 
deflections  for  flight  conditions  with  similar  dynamic  pres¬ 
sures  can  be  found  in  the  aerodynamic  data.  When  examining 
the  aerodynamic  data  obtained  for  the  flight  conditions  of 


that  the  trim  values  of  angle  of  attack  differed.  Table 


6-5  shows  the  angle  of  attack  for  each  of  the  flight  condi¬ 
tions.  Figure  6-30  shows  that  the  lift  coefficient  varies 
as  the  angle  of  attack  changes.  When  the  lift  coefficient 
changes,  the  stability  derivatives  associated  with  it 
change,  and  hence  so  do  the  dimens i onal i zed  derivatives.  The 
pitching  moment  also  changes  as  the  lift  coefficient  chan¬ 
ges.  Therefore,  the  stability  derivatives  associated  with 
the  pitching  moment  will  differ  also. 

From  the  above  discussion,  it  can  be  summarized  that 
there  are  some  differences  in  the  responses  of  the  control 
surfaces  for  models  having  the  same  dynamic  pressure.  Even 
with  these  differences,  model  equivalency  is  primarily  as¬ 
sociated  with  dynamic  pressure.  This  was  shown  by  generating 

Table  6-5 


Angle  of  Attack  for  Nominal  Flight  Conditions 


Flight  Condition 

Angle  of  Attack  (deg) 

38,000  ft 

5.11515 

0 . 7  Mach 

Dynamic  Pressure: 1 4 8 . 4 

18,000  ft 

4.96835 

0.45  Mach 

Dynamic  Pressure:  1 49 . 9 

c 
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Figure  6-30.  Lift  Curve 


performance  boundaries  for  several  pairs  of  models  with 
similar  values  of  dynamic  pressure  (see  Table  6-3)  and 


observing  their  similarity. 


6.2.  1.2  Control  Surface  Activity.  The  effects  on  the  con¬ 
trol  surfaces  of  replacing  a  nominal  model  with  plants 
having  increased/decreased  Mach  number  and/or  altitude  can 
be  made  clear  by  examining  the  changes  in  the  aerodynamic 
gain.  For  a  given  simulation,  the  system  gain  is  composed  of 
the  control  law  gains  and  the  aerodynamic  gain.  The  aerody¬ 
namic  gain  is  derived  from  the  state  space  model  for  a 
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specific  flight  condition.  As  an  example,  the  aerodynamic 
gain  for  the  flight  condition  of  30,000  ft  and  1.2  Mach  is 
shown  in  the  discussion  to  follow.  Recall  that  the  open 
loop  aircraft  model  is  longitudinally  unstable  in  the  sub¬ 
sonic  region.  The  reason  the  flight  condition  chosen  was  in 
the  supersonic  region  is  that  the  open  loop  plant  is  stable. 
Therefore,  the  resulting  transfer  functions  that  were  used 
to  calculate  the  aerodynamic  gain  will  have  all  of  the  poles 
in  the  left  half  plane. 

The  aerodynamic  gains  for  the  30,000  ft, 1.2  Mach  flight 
condition  were  calculated  for  the  following  transfer  func¬ 
tions:  q/g  ,  q/g  .  ®/6  »  andQ!/6_.  The  state  aircraft  state 
e  i  e  t 

space  model  used  was 
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and  the  output  equation  was 
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The  proper  dimensionalized  stability  derivatives  were  ob¬ 
tained  and  inserted  into  the  state  space  model  above.  The 


desired  transfer  functions  were  then  obtained  and  are  pre- 


The  steady  state  aerodynamic  gains  were  calculated  for  the 
transfer  functions  of  Equations  (6-3)  through  (6-6).  The 
aerodynamic  gains  for  the  above  transfer  functions  are  pre 
sented  in  Table  6-6. 

Next,  the  steady  state  aerodynamic  gains  for  a  flight 
condition  at  a  different  altitude  were  calculated.  The 
flight  condition  was  10,000  ft  and  1.2  Mach.  The  transfer 
functions  for  this  flight  condition  were  calculated  as 


c 


Table  6-6 


Aerodynamic  Gains  for  Transfer  Functions 
of  Equations  (6-3)  through  (6-6) 


Transfer  Function 

Magnitude  of  Gain 

q/'b 

0.7845 

(Eq  6-3) 

q/6f 

0.2034 

(Eq  6-4 ) 

a/i>e 

0.7237 

( Eq  6-5  ) 

a/bf 

0.2464 

( Eq  6-6 ) 

before  and  are  presented  in  the  following  equations 

-60.5221(s  +  2.150764) 


+  3.460089s  +  80.19058 

19.4011(s  +  1.79332) 


+  3.460089s  +  80.19058 
. 259837 ( s  +  229.2444) 

+  3.460089s  +  80.19058 

.170677(s  +  112.372) 

+  3.460089s  +  80.19058 


(6-10) 


(6-11  ) 


As  before,  the  steady  state  aerodynamic  gains  were 


calculated  and  are  presented  in  Table  6-7.  By  comparing  the 
aerodynamic  gains  of  Table  6-6  and  6-7,  it  is  seen  that  the 
gains  for  the  pitch  rate  transfer  functions  are  much  larger 
for  the  altitude  of  10,000  ft.  The  gains  for  the  transfer 
functions  with  angle  of  attack  as  a  parameter  are  approxi¬ 
mately  the  same.  Because  this  thesis  uses  pitch  rate  as  a 
variable  to  be  tracked,  the  fact  that  the  aerodynamic  gains 
for  the  higher  altitude  flight  condition  are  larger  is 
significant.  Therefore,  an  increase  in  altitude  results  in 
an  decreased  aerodynamic  gain  and  hence  faster  surface  acti¬ 
vity.  Similar  results  were  obtained  when  considering  the 

Table  6-7 

Aerodynamic  Gains  for  Transfer  Functions 
of  Equations  (6-8)  through  (6-11) 

Transfer  Function  Magnitude  of  Gain 


effects  on  the  aerodynamic  gain  when  increasing  Mach  number. 
The  data  showed  that,  when  Mach  number  and/or  altitude 
increased,  the  aerodynamic  gain  decreased  accordingly. 

6.2.2  Two-Model  Configurations 

In  determining  the  number  of  models  required  to  allow 
the  system  with  the  multiple  model  algorithm  to  yield  satis¬ 
factory  tracking  performance  over  the  desired  flight  enve¬ 
lope,  several  two-model  configurations  were  evaluated.  The 
term  "two-model  configuration”  refers  to  placing  two  air¬ 
craft  models  in  the  parallel  bank  to  function  as  the  secon¬ 
dary  parameter  estimators  (see  Figure  4-1).  The  two-model 
configuration  evaluation  consisted  primarily  of  examining 
the  amount  the  performance  boundaries  were  required  to  over¬ 
lap  to  yield  proper  tracking  performance.  The  two  models 
for  each  configuration  were  chosen  on  different  dynamic 
pressure  lines  for  the  reasons  discussed  earlier.  The  two- 
model  configurations  are  listed  on  Table  6-8. 

The  performance  boundaries  were  determined  as  described 
previously  for  each  of  the  configurations  listed  in  Table  6- 
8.  Figures  6-31  through  6-35  show  the  performance  bounda¬ 
ries  for  each  of  the  two-model  configurations. 


Table  6-8 


Two  Model  Configurations 


Configuration  Nominal  Model  Flight  Conditions 

Alt  (Kft)  Mach  Dynamic  Pressure 


1 

26 

0.55 

154.39 

18 

0.80 

473.75 

2 

10 

0.50 

254 .73 

10 

0.75 

573.14 

3 

38 

0.65 

127.95 

10 

0.65 

430.50 

4 

10 

0.50 

254 . 73 

10 

0.65 

430.50 

5 

10 

0.35 

124.81 

10 

0.75 

573.14 

From  Figures  6-31  through  6-35  it  is  apparent  that  the 
flight  envelope  of  interest  can  be  encompassed  by  the  union 
of  the  performance  boundaries  of  two  nominal  models  (depen¬ 
ding  on  the  position  of  the  individual  performance  bounda¬ 
ries).  However,  in  attempting  to  cover  the  desired  flight 
envelope  with  two  models,  several  observations  were  made  as 
to  the  effects  of  model  placement  on  the  performance  of  the 
multiple  model  algorithm. 


Due  to  the  fact  that  the  multiple  model  algorithm 
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Figure  6-31.  Performance  Boundaries 
for  Configuration  1 
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Nominal  Model  at  10K  0 . 5M 
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Figure  6-32.  Performance  Boundarie 
for  Conf i gura t i on  2 
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Figure  6-34.  Performance  Boundaries 
for  Configuration  4 
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arrives  at  the  estimated  value  of  the  portion  of  the 
parameter  vector  by  a  blending  process  of  the  hypothesis 
probabilities  {see  Chapter  4),  the  amount  of  overlap  of  the 
performance  boundaries  is  important.  Figures  6-36  through 
6-44  show  the  significance  of  the  overlap  of  the  performance 
boundaries.  It  is  seen  from  Figures  6-36  through  6-44  that 
although  the  operating  point  of  each  of  the  examples  lies 
within  the  performance  boundaries  (see  Figures  6-32 
through  6-34)  the  control  law  gains  that  were  implemented 
caused  the  performance  criteria  to  be  violated. 

The  reason  that  the  performance  criteria  were  violated, 
was  due  to  the  manner  in  which  the  multiple  model  algorithm 
was  implemented.  By  looking  at  the  plots  of  the  hypothesis 
probabilities  (see  Figures  6-38,41,44),  it  can  be  seen  that 
the  multiple  model  algorithm  calculated  a  probability  of 
approximately  one  for  one  of  the  two  models,  while  the  other 
model  was  given  a  probability  of  approximately  zero.  If  the 
probabilties  were  exactly  one  and  zero  for  the  respective 
cases,  the  time  responses  for  the  adaptive  case  would  be  the 
same  as  those  generated  for  the  fixed  gain  case.  The  fact 
that  the  responses  were  different  for  the  adaptive  versus 
the  fixed  gain  system  showed  that  the  probabilities  were  not 
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Figure  6-36.  Pitch  Rate  Response 
Configuration  2/Operating  Point:  14K  0.35M 
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Figure  6-37.  Pitch  Rate  Performance  Criterion 
Configuration  2/Operating  Point:  14K  0.35M 
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Figure  6-38.  Model  Probability  Weightings 
Configuration  2/Operating  Point:  14K  0.35M 
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Figure  6-39.  Pitch  Rate  Response 
Configuration  4/Operating  Point:  14K  0 . 9M 
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Figure  6-40.  Pitch  Rate  Performance  Criterion 
Configuration  4/Operating  Point:  14K  0 . 9M 
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Figure  6-41.  Model  Probability  Weightings 
Configuration  4/Operating  Point:  14K  0 . 9M 
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Figure  6-42.  Pitch  Rate  Response 
Configuration  3/Operating  Point:  10K  0 . 9M 
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Figure  6-43.  Pitch  Rate  Performance  Criterion 
Configuration  3/Operating  Point:  10K  0 . 9M 
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The  multiple  model  algorithm  is  implemented  with  a 
lower  probability  bound  so  the  probability  of  any  given 
candidate  model  will  not  be  set  to  zero.  The  reason  for  this 
is  that  the  multiple  model  algorithm  is  updated  recursively 
and,  if  a  candidate  model’s  probability  goes  to  zero,  the 
contribution  of  that  model  will  be  eliminated  from  that 
point  forward. 

Lower  bounding  the  candidate  model  probabilities  means 
that  every  candidate  model  will  have  an  effect  on  the  con¬ 
trol  law  gains  to  some  degree.  This  means  that,  even  if  the 
operating  point  of  the  system  is  at  a  candidate  model  posi¬ 
tion  the  control  law  gains  will  differ  slightly  from  the 
gains  of  the  specific  candidate  model.  When  the  candidate 
models  are  placed  on  constant  dynamic  pressure  lines  that 
are  'too  far  apart’  and  the  operating  point  is  close  to  the 
performance  boundary  of  one  of  the  models  (see  Figure  6- 
32)  the  gains  may  become  such  that  the  performance  criteria 
is  violated.  By  referring  to  Figure  6-32,  it  can  be  seen 
that  the  operating  point  is  well  outside  the  perfomance 
boundary  of  model  2  and  on  the  boundary  of  model  1.  Al¬ 
though  model  1  has  a  high  probability  associated  with  it  as 


the  correct  model,  by  including  the  small  contribution  of 
the  control  law  gains  associated  with  model  2  the  perfor¬ 
mance  criteria  was  violated.  This  is  shown  in  the  data 
presented  on  Figures  6-36  through  6-38. 

By  running  simulations  with  different  values  of  the 
lower  bound  for  the  weighting  factors,  it  was  seen  that 
there  is  an  interplay  between  the  amount  of  overlap  required 
for  the  performance  boundaries  and  the  lower  bound  size. 

The  smaller  the  value  of  the  lower  bound  the  less  overlap 
was  required  of  the  performance  boundaries  (so  that  satis¬ 
factory  tracking  performance  would  be  obtained  for  all  the 
flight  conditions  within  the  union  of  the  performance  boun¬ 
daries  )  . 

By  using  a  lower  bound  value  for  the  weighting  factors 
of  0.01  it  was  determined  that  the  candidate  models  need  to 
be  placed  so  that  there  is  a  sufficient  amount  of  overlap 
between  the  performance  boundaries.  As  shown  on  Figure 
6-45,  when  the  boundaries  overlap  it  is  preferable  to  have 
the  boundary  of  one  model  well  inside  the  performance  bound¬ 
ary  of  the  second  model  to  minimize  the  problem  discussed 
above.  Also,  the  position  of  the  candidate  models  that  will 
cover  the  outer  boundaries  of  the  desired  flight  envelope 


#  Location  of  Nominal  Model  1 

*  Location  of  Npmin  d  Model  2 


—  Performance  Boundary 
of  Model  1 

--Performance  Boundary 
of  Model  2 


MACH 


Figure  6-45.  Performance  Boundary  Overlap 


need  to  be  placed  so  the  outer  boundaries  of  the  flight 
envelope  are  well  inside  the  specific  candidate  models’ 
performance  boundaries  (see  Figure  6-45).  Because  the  model 
placement  requirements  could  be  better  achieved  with  three 
models  instead  of  two,  several  three-model  configurations 
were  evaluated  in  an  attempt  to  achieve  the  desired  tracking 
performance  over  the  flight  envelope  of  interest. 
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6.2.3  Three-Model  Con figurations 

The  three-model  configurations  that  were  evaluated  are 
presented  in  Table  6-9.  A  graphical  presentation  of  the 
three-model  configurations  are  given  in  Figure  6-46.  From 
the  beginning  of  the  three-model  configuration  evaluation, 
it  was  clear  that  the  desired  tracking  performance  could  be 
obtained  over  the  flight  envelope  of  interest  with  three 
mode  1 s . 

The  first  three-model  configurations  that  were  eval¬ 
uated  were  configurations  6,  7,  and  8  (see  Table  6-9).  The 
models  for  configurations  6,  7,  and  8,  were  chosen  so  that 
the  values  of  dynamic  pressures  for  the  candidate  models 
would  be  as  close  as  possible  with  the  data  used  for  this 
thesis.  The  reason  for  doing  this  was  to  examine  the  effect 
of  choosing  different  nominal  models  (with  essentially  the 
same  dynamic  pressures)  on  the  multiple  model  algorithm’s 
perf  ormance . 

The  performance  boundaries  for  configurations  6,  7,  and 
8  are  presented  on  Figures  6-47,  6-48,  and  6-49.  They  show 
that  the  performance  boundaries  for  each  configuration  are  ap¬ 
proximately  the  same.  This  equivalency  of  performance 
boundaries  is  due  to  the  fact  that  each  configuration  con- 
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Table  6-9 


Three-Model  Configurations 


Configuration 

Nominal  Model  Flight  Conditions 

Alt  (Kft)  Mach  Dynamic  Pressure 

6 

38 

0.65 

127.95 

26 

0.70 

258.18 

10 

0.75 

573  .  14 

7 

30 

0.50 

110.19 

22 

0.65 

264 . 58 

14 

0.80 

557 . 18 

8 

10 

0.35 

124.81 

10 

0.50 

254 . 73 

10 

0.75 

573.14 

9 

22 

0.45 

126.81 

22 

0.65 

264.58 

22 

0.90 

507.24 

10 

10 

0.35 

124.81 

38 

0.80 

193.82 

10 

0.75 

573.14 

1 1 

10 

0.35 

124.81 

22 

0.65 

264 . 58 

10 

0.75 

573.14 

12 

10 

0 . 35 

124.81 

26 

0.70 

258.18 

10 

0.75 

573. 14 

13 

22 

0.45 

126.81 

22 

0.65 

264 . 58 

10 

0.75 

573.14 

14 

22 

0.45 

126.81 

10 

0.50 

254.73 

10 

0.75 

573.14 
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Figure  6-47.  Performance  Boundaries 
for  Configuration  6 
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Figure  6-48.  Performance  Boundaries 
for  Configuration  7 


sists  of  models  with  approximately  the  same  dynamic  pres¬ 
sures.  The  reason  for  the  equivalency  of  the  performance 
boundaries  of  models  with  the  same  dynamic  pressure  was 
discussed  earlier  in  this  chapter. 

Following  the  establishment  of  the  performance  bounda¬ 
ries  of  configurations  6,  7,  and  8,  simulations  were  run 
using  the  multiple  model  algorithm  for  every  flight  condi¬ 
tion  in  the  flight  envelope  of  interest.  The  results  of 
those  simulations  are  presented  on  Tables  6-10  through  6-12 
and  on  Figures  6-50  through  6-52. 

Figures  6-50  through  6-52  show  graphically  which  model 
or  models  the  multiple  model  algorithm  selected  as  the  model 
to  represent  the  current  flight  condition.  It  should  be 
pointed  out  that  the  performance  criteria  were  met  for  every 
flight  condition  with  the  nominal  models  of  configurations 
6,  7,  and  8.  As  an  example  of  how  the  individual  entries 
for  each  flight  condition  were  determined  for  Figures  6-50 
through  6-52,  consider  the  flight  condition  26,000  feet  and 
0.45  Mach  for  configuration  6.  By  looking  at  the  plot  of 
the  conditional  probabilities  for  the  candidate  models  (see 
Figure  6-53),  it  can  be  seen  that  the  probability  of  candi¬ 
date  model  one  (38,000  ft,  0.65  Mach)  is  above  0.9.  There- 
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Figure  6-50.  Model  Selection  Data 
for  Configuration  6 
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Figure  6-52.  Model  Selection  Data 
for  Configuration  8 
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-  -  Model  at  38K  0.65M 

--  Model  at  26K  0.70M 
• •  Model  at  10K  0.75M 


Time  (Sec) 

Figure  6-53.  Model  Probability  Weightings 
Configuration  6/Operating  Point:  26K  0.45M 


fore  a  *  is  placed  at  the  flight  condition  of  26,000  feet, 


<0 


0.45  Mach  on  Figure  6-50.  If  the  conditional  probability 
fell  below  0.9  for  a  candidate  model  during  the  running  of  a 
simulation,  then  the  two  models  with  the  highest  probability 
were  placed  at  the  position  of  the  flight  condition  of 
interest.  As  an  example,  the  conditional  probabilities  for 
the  flight  condition  of  26,000  feet,  0.65  Mach  (configuration 
8)  are  shown  on  Figure  6-54.  Figure  6-54  shows  that  the 
conditional  probability  of  model  one  falls  below  0.9  and 
that  the  models  with  the  most  significant  probabilities  are 
models  one  and  two,  respectively.  Therefore,  a  **  was 
placed  on  the  location  of  the  26,000  feet,  0.65  Mach  flight 
condition  of  Figure  6-52.  This  annotation  delineated  which 
candidate  model (s)  had  the  highest  conditional  probability 
at  the  given  flight  condition. 

By  comparing  Figures  6-50  through  6-52,  it  can  be  seen 
that  the  multiple  model  algorithm  performance  was  similar 
for  configurations  6,  7,  and  8.  This  similarity  indicated 
that  the  conditional  probability  calculated  for  a  candidate 
model  was  associated  with  dynamic  pressure.  Tables  6-10 
through  6-12  show  the  candidate  models  that  were  selected  as 
closest  to  the  actual  flight  condition  as  well  as  the  dif- 
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Figure  6-54.  Model  Probability  Weightings 
Configuration  8/Operating  Point:  26K  0.65M 


•itv 


ferences  between  the  dynamic  pressure  of  the  flight  condi¬ 


tion  and  the  candidate  models.  The  data  showed  that,  in 
over  ninety  percent  of  the  three  hundred  cases  run,  the 
candidate  model  that  received  the  highest  conditional  proba¬ 
bility  was  the  model  with  a  dynamic  pressure  closest  to  the 
dynamic  pressure  of  the  actual  flight  condition.  Again, 
this  pointed  out  the  fact  that  model  equivalency  is  largely 
dependent  on  dynamic  pressure. 

As  just  discussed,  configurations  6,  7,  and  8  all 
yielded  satisfactory  tracking  performance  (from  a  perfor¬ 
mance  criteria  standpoint)  over  the  flight  envelope  of  inte¬ 
rest.  However,  when  comparing  the  performance  index  respon¬ 
ses  of  configurations  6,  7,  and  8,  it  was  determined  that 
configuration  8  resulted  in  the  least  tracking  error  for  a 
majority  of  the  flight  conditions.  Table  6-13  presents  a 
representative  sample  of  the  performance  index  results  for 
configurations  6,  7,  and  8.  Configuration  6  had  the  largest 
tracking  error  as  determined  by  the  performance  index  for 
all  of  the  flight  conditions.  As  seen  from  Table  6-13, 
configuration  8  was  the  best  of  the  configurations  labled  as 
6,  7,  and  8  and  configuration  6  was  the  least  desirable. 

In  an  effort  to  obtain  a  three-model  configuration  that 
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Table  6-13 


C 


Performance  Index  Comparison  of  Selected 
Flight  Conditions 


Flight  Condition 

Conf iguration 

with  Minimal  Error 

Alt  (Kft) 

Mach 

6 

7 

8 

10 

0.35 

X 

10 

0.65 

X 

10 

0.90 

X 

14 

0.35 

X 

14 

0.65 

X 

14 

0.90 

X 

22 

0.40 

X 

22 

0.65 

X 

22 

0.90 

X 

26 

0.45 

X 

26 

0.65 

X 

26 

0.90 

X 

30 

0.50 

X 

30 

0.90 

X 

38 

0.60 

X 

38 

0.90 

X 

yielded  better  performance  than  configuration  8, 

several 

other  three 

-model 

conf igutrations  were 

examined . 

Specif i - 

cally,  these  were  configurations  9  through  14  as  given  on 
Table  6-9.  In  each  case,  the  overall  performance 
of  configuration  8  proved  to  be  the  superior  configuration. 

Configuration  9  had  performance  index  responses  that 
were  very  similar  to  those  of  configuration  8.  Although 
configuration  8  did  have  a  minimal  amount  of  error  most  of 
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the  time,  configuration  9  had  less  error  at  the  low  alti- 


i 


tude/high  Mach  flight  conditions.  Configurations  11,  13, 
and  14  had  performance  index  responses  that  were  less  than 
those  for  configuration  8  at  38,000  ft,  0.9  Mach.  However, 
configuration  8  had  better  responses  over  the  rest  of  the 
flight  envelope.  So,  from  a  performance  criteria  stand¬ 
point,  configuration  8  was  the  best  configuration  of  those 
investigated  with  configuration  9  yielding  similar  perfor¬ 


mance  . 


In  an  effort  to  better  the  performance  of  configuration 
8,  a  four-model  configuration  was  evaluated  with  all  of  the 
models  at  10,000  ft  and  Mach  numbers  of  0.35,  0.5,  0.75,  and 
0.9,  respectively.  However,  the  performance  of  this  config¬ 
uration  was  not  any  better  than  that  afforded  by  the  three 
models  of  configuration  8. 


6.2.4  Control  law  Gain  Adjustment.  In  his  thesis,  Barfield 
(21)  noted  that  although  the  Porter  design  technique  yields 
exceptionally  fast  aircraft  responses,  many  times  this  is 
achieved  through  large  and  extrememly  fast  surface  move¬ 
ments.  After  a  review  of  several  applications  of  the  design 
technique,  it  became  apparent  that  for  a  given  design,  the 


larger  the  control  law  gains,  the  faster  the  surface  move¬ 
ments.  Reaching  either  the  position  limits  or  the  rate 
limits  is  a  condition  that  should  be  avoided.  Reaching  a 
position  limit  prevents  the  input  command  from  being  a- 
chieved.  Reaching  rate  limits  results  in  the  system  becom¬ 
ing  nonlinear  and  can  result  in  an  increase  in  phase  lag  in 
a  feedback  system.  This  phase  lag  increase  has  been  shown 
to  cause  loss  of  control  in  aircraft  with  a  full  authority 
fly-by-wire  system  (21). 

Recall,  from  Chapter  3,  that  the  control  law  gains  are 
given  by 

K i  =  H_1(T)|  (6-12) 

K2  =  p  K  (6-13) 

where  the  diagonal  weighting  matrix,  s,  can  be  altered  to 
achieve  desired  tracking  characteristics  and  the  parameter, 
p  , is  a  constant  which  assigns  the  ratio  of  proportional  to 
integral  control  . 

For  all  of  the  simulations  that  were  accomplished  up  to 
this  part  of  this  thesis,  the  control  law  design  parameters 
that  were  implemented  were  those  used  by  Pineiro  and  Berens. 
The  control  law  design  parameters  had  the  following  values: 


£ 


0 . 4 
0.0 


0.0 

0.7 


and 


P  =  0.8 


(6-14  ) 


(  6-15  ) 


As  is  shown  earlier,  configuration  8  had  satisfactory 
performance  over  the  flight  envelope  of  interest  (from  a 
performance  index  standpoint).  However,  for  several  of  the 
flight  conditions  the  rate  limits  were  encountered  for  both 
the  elevators  and  flaperons.  In  an  effort  to  reduce  the 
control  surface  activity  for  the  multiple  model  algorithm, 
an  investigation  into  the  effects  of  altering  the  control 
law  design  parameters  was  performed. 

To  begin,  the  effects  of  altering  the  control  law- 
design  parameters  on  the  performance  of  a  fixed  gain  simula¬ 
tion  were  investigated.  The  nominal  flight  condition  chosen 
for  this  part  of  the  analysis  was  22,000  ft,  0.65  Mach.  The 
responses  for  the  fixed  gain  simulation  using  the  control 
law  design  parameters  values  as  shown  above  are  shown  Fig¬ 
ures  6-55  through  6-60.  As  can  be  seen  on  Figure  6-60,  the 
rate  limits  for  the  flaperons  were  encountered  several 
times.  For  the  reasons  given  above,  this  is  a  condition 
that  should  be  avoided. 
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Figure  6-59.  Flaperon  Rate 
Nominal  Model:  22K  0 . 6 5M/0perat i ng  Point:  22K  0.65M 
SIG1 =0 . 4  SIG2  =  0 . 7  RHO=0.8 


Time  (Sec) 

Figure  6-60.  Elevator  Rate 
Nominal  Model:  22K  0 . 65M/Operating  Point:  22K  0.65M 
S IG 1  =  0 . 4  S IG2  =  0 . 7  RHO=0.8 


Table  6-14  lists  the  simulations  that  were  run  subse- 


r 


quent  to  the  initial  22,000  ft,  0.65  Mach  simulation  and  the 
respective  values  of  the  control  law  design  parameters.  The 
entries  labled  Sigma  1  and  Sigma  2  refer  to  the  1,1  and  2,2 
elements  of  the  diagonal  weighting  matrix,  respectively 
while  the  term  rho  refers  to  the  design  parameter,  p. 

As  can  be  seen  from  Table  6-14,  simulations  1  and  2 
were  run  with  varying  values  of  the  parameter,  p,  simula¬ 
tions  3  through  6  were  performed  with  different  values  of 
the  1,1  element  of  the  weighting  matrix,  and  simulations  7 
through  9  were  accomplished  with  different  values  of  the  2,2 


element  of  the  weighting  matrix. 


Table  6-14 


Design  Parameter  Values  for  Simulations 
at  22,000  ft,  0.65  Mach 


Simulation 

Number 

Control 
Sigma  1 

Law  Design 
Sigma  2 

Parameters 

Rho 

1 

0.40 

0.70 

0.40 

2 

0.40 

0.70 

0.10 

3 

0.80 

0.70 

0.80 

4 

0.20 

0.70 

0.80 

5 

0. 10 

0.70 

0.80 

6 

0.05 

0.70 

0.80 

7 

0.40 

1 . 50 

0.80 

8 

0.40 

0.35 

0.80 

9 

0.40 

0.10 

0.80 
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Decreasing  the  parameter,  p,  resulted  in  an  increase  in 
the  tracking  error  as  shown  on  the  performance  index  respon¬ 
ses.  Figures  6-61  through  6-66  present  the  simulation  re¬ 
sults  when  the  parameter,  p,  was  set  to  0.4  (simulation  1  on 
Table  6-14).  By  comparing  the  responses  of  the  elevators 
and  f laperons  for  the  original  simulation  and  those  of 
simulation  1,  it  can  be  seen  that  there  is  little  difference 
in  control  surface  performance.  Therefore,  decreasing  the 
ratio  of  proportional  to  integral  control  results  in  an 
increase  in  the  tracking  error  as  determined  by  the  perfor¬ 
mance  index  while  resulting  in  minimal  changes  in  the  con¬ 
trol  surface  performances. 

Figures  6-67  through  6-72  present  the  responses  for 
simulation  3  on  Table  6-14.  As  shown  on  Figure  6-68,  in¬ 
creasing  the  value  of  Sigma  1  resulted  in  an  increase  in  the 
performance  error.  Also,  the  increase  in  Sigma  1  resulted 
in  the  control  surfaces  becoming  more  active  as  evidenced  by 
the  deflection  rates.  The  negative  aspect  of  this  is  that 
the  control  surfaces  spent  a  larger  amount  of  time  at  the 
rate  limit.  Conversely,  reducing  Sigma  1  as  in  simulation  5 
(see  Figures  6-73  through  6-78)  resulted  in  a  performance 
error  that  was  essentially  the  same  as  that  of  the  original 
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Figure  6-65.  Flaperon  Rate 
Nominal  Model:  22K  0 . 6 5M/0perat i ng  Point:  22K  0.65M 
SIG1=0 . 4  S IG2  =  0 . 7  RHO=0.4 
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Figure  6-66.  Elevator  Rate 
Nominal  Model:  22K  0 . 65M/Operating  Point:  22K  0.65M 
SIG1 =0 . 4  SIG2=0 . 7  RHO=0.4 
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Figure  6-67.  Pitch  Rate  Response 
Nominal  Model:  22K  0 . 65M/Operating  Point:  22K  0.65M 
SIG 1 =0 . 8  SIG2=0 . 7  RHO=0.8 
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Figure  6-68.  Pitch  Rate  Performance  Criterion 
Nominal  Model:  22K  0 . 65M/Operating  Point:  22K  0.65M 
SIG1 =0 . 8  SIG2=0 . 7  RHO=0.8 
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Figure  6-75.  Flaperon  Position 
Nominal  Model:  22K  0 . 65M/Operating  Point:  22K 
SIG1  =  0.1  SIG2  =  0 . 7  RHO  =  0 . 8 
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Figure  6-76.  Elevator  Position 
Nominal  Model:  22K  0 . 65M/Operating  Point:  22K  0.65M 
SIGlsO.l  S IG2  =  0 . 7  RHO=0 . 8 
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Figure  6-77.  Flaperon  Rate 
Nominal  Model:  22K  0 . 65M/Operating  Point:  22K  0.65M 
SIG1=0.1  SIG2=0 . 7  RHO=0 . 8 
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Figure  6-78.  Elevator  Rate 
Nominal  Model:  22K  0 . 65M/Operating  Point:  22K  0.65M 
SIG1=0.1  SIG2=0 . 7  RHO=0 . 8 
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case.  However,  the  elevator  and  flaperon  rates  were  much 


improved.  The  flaperon  rates  never  reached  rate  limit  for 
the  entire  simulation.  The  important  point  here  is  that, 
although  the  performance  index  was  not  affected  significant¬ 
ly,  decreasing  Sigma  1  resulted  in  improved  control  surface 
activity . 

Increasing  Sigma  2  resulted  in  less  error  as  shown  on 
Figures  6-79  through  6-84  for  simulation  7,  but  the  control 
surfaces  were  more  active,  as  can  be  seen  on  the  elavator 
and  flaperon  positions  and  rates.  Decreasing  Sigma  2  had 
the  effect  of  increasing  the  error  as  measured  by  the  per¬ 
formance  criteria.  This  is  shown  on  Figures  6-85  through  6- 
90  which  are  simulations  for  simulation  9  on  Table  6-14. 
Reducing  Sigma  2  did  not  have  much  effect  on  the  control 
surface  rates,  although  Figure  6-89  shows  that  the  amount  of 
flaperon  required  to  follow  the  input  increased.  The  data 
showed  that  the  ability  of  the  output  to  tightly  track  the 
input  was  affected  significantly  by  Sigma  2. 

Following  the  determination  of  the  effects  of  varying 
the  control  law  design  parameters,  Sigma  1  was  set  equal  to 
0.1  to  improve  the  control  surface  responses.  With  Sigma  1 
equal  to  0.1,  simulations  were  run  to  determine  the  perfor- 
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Figure  6-79.  Pitch  Rate  Response 
Nominal  Model:  22K  0 . 65M/Operating  Point:  22K  0.65M 
SIG1 =0 . 4  SIG2=1 . 5  RHO=0.8 
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Figure  6-80.  Pitch  Rate  Performance  Criterion 
Nominal  Model:  22K  0 . 65M/Operating  Point:  22K  0.65M 
SIGU0.4  SIG2  =  1 .5  RHOsO.8 
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Figure  6-84.  Elevator  Rate 
Nominal  Model:  22K  0 . 65M/Operating  Point:  22K  0.65M 
SIG1=0 . 4  SIG2= 1 . 5  RHO=0 . 8 
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Figure  6-85.  Pitch  Rate  Response 
Nominal  Model:  22K  0 . 65M/Operating  Point:  22K  0.65M 
SIG1=0 . 4  SIG2=0 . 1  RHO=0 . 8 
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Figure  6-86.  Pitch  Rate  Performance  Criterion 
Nominal  Model:  22K  0 . 65M/Operating  Point:  22K  0.65M 
SIG1=0.4  SIG2=0.1  RHO=0 . 8 
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Figure  6-87.  Flaperon  Position 
Nominal  Model:  22K  0 . 65M/Operating  Point:  22K  0.65M 
SIG1=0.4  SIG2=0.1  RHO=0 . 8 
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Figure  6-88.  Elevator  Position 
Nominal  Model:  22K  0 . 65M/Operating  Point:  22K  0.65M 
SIG1=0.4  SIG2=0 . 1  RHO=0 . 8 
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Figure  6-89.  Flaperon  Rate 
Nominal  Model:  22K  0 . 65M/Operating  Point:  22K  0.65M 
SIG1 =0 . 4  SIG2=0 . 1  RHOrO.8 
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Figure  6-90.  Elevator  Rate 
Nominal  Model:  22K  0 . 65M/Operating  Point:  22K  0.65M 
SIG1=0 . 4  SIG2=0.1  RHOrO . 8 
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stance  boundary  for  the  nominal  flight  condition  of  22,000 
ft,  0.65  Mach.  The  data  obtained  from  these  simulations 
showed  that  the  performance  boundary  was  not  affected  by  the 
reduction  of  the  value  of  Sigma  1.  The  performance  boundary 
for  the  decreased  Sigma  1  case  was  the  same  as  that  shown  on 
Figure  6-48.  Figure  6-48  shows  the  performance  boundary  for 
the  22,000  ft,  0.65  Mach  flight  condition  with  the  control 
law  design  parameters  as  given  in  Equations  (6-14)  and  (6- 
15).  Performance  boundaries  for  nominal  flight  conditions 
of  18,000  ft,  0.45  Mach  and  38,000  ft,  0.7  Mach  were  also 
determined  with  Sigma  1  set  equal  to  0.1.  As  was  the  case 
earlier,  the  performance  boundaries  with  Sigma  1  equal  to 
0.1  were  the  same  as  those  determined  previously. 

Due  to  the  fact  that  the  performance  boundaries  re¬ 
mained  the  same  when  decreasing  the  value  of  Sigma  1,  the 
three-model  configuration,  configuration  8,  had  performance 
boundaries  that  covered  the  flight  envelope  of  interest  as 
before  (see  Figure  6-49).  By  using  configuration  8  as  the 
models  in  the  parallel  bank  of  secondary  estimators  for  the 
multiple  model  algorithm  (with  Sigma  1=0.1),  satisfactory 
tracking  performance  over  the  flight  envelope  of  interest 
(from  a  performance  criteria  standpoint)  was  achieved  as 
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before.  However,  the  control  surface  responses  were  much 
improved  with  the  value  of  Sigma  1  set  equal  to  0 . 1 .  In 
contrast  to  the  simulations  for  configuration  8  with  the 
gains  of  Equations  (6-14)  and  (6-15),  the  rate  limits  were 
never  encountered  when  Sigma  1  was  set  equal  to  0.1.  Fig¬ 
ures  6-91  through  6-102  present  simulation  data  for  config¬ 
uration  8  with  Sigma  1  equal  to  0.1.  Note  the  responses  of 
the  control  surfaces  and  the  fact  that  they  do  not  reach  the 
rate  limits  for  either  the  elevators  or  flaperons. 

6.2.5  Sensor  Noise  Effects 

To  determine  the  effects  of  sensor  noise  on  the  perfor¬ 
mance  on  the  multiple  model  algorithm  as  the  means  of  para¬ 
meter  identification,  simulations  were  run  where  independent 
white,  gaussian  noise  was  injected  into  each  of  the  quanti¬ 
ties  of  interest.  The  noise  levels  used  in  this  thesis  are 
realistic  noise  levels  for  a  comparable  aircraft  (Grumman 
F-14  Tomcat)  and  are  presented  in  Table  6-15. 

To  gain  some  insight  into  the  effects  of  noise  on  the 
multiple  model  algorithms  performance,  the  two-model  con¬ 
figuration  shown  in  Table  6-16  was  considered.  For  the 
configuration  in  Table  6-16,  an  operating  point  of  10,000 
ft,  0.35  Mach  was  selected  which  was  one  of  the  models  in 
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Figure  6-91.  Pitch  Rate  Response 
Configuration  8/Operating  Point:  10K  0 . 9M 
SIG1=0.1  SIG2=0 . 7  RHO=0 . 8 
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Figure  6-92.  Pitch  Rate  Performance  Criterion 
Configuration  8/Operating  Point:  10K  0.9M 
SIGlsO.l  SIG2=0 . 7  RHO=0 . 8 
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Figure  6-95.  Flaperon  Rate 
Configuration  8/Operating  Point:  10K  0 . 9M 
S IG1 =0 . 1  S IG2=0 . 7  RHO=0 . 8 
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Figure  6-96.  Elevator  Rate 
Configuration  8/Operating  Point:  10K  0 . 9M 
S IG1 =0 . 1  SIG2=0.7  RHO=0 . 8 
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Figure  6-97.  Pitch  Rate  Response 
Configuration  8/Operating  Point:  38K  0 . 9M 
SIGlrO.l  S IG2  =  0 . 7  RHO=0.7 


Figure  6-98.  Pitch  Rate  Performance  Criterion 
Configuration  8/Operating  Point:  38K  0 . 9M 
S IG1 =0.1  SIG2=0 . 7  RHO=0 . 8 
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Figure  6-99.  Flaperon  Position 
Configuration  8/Operating  Point:  38K  0 . 9M 
SIG1 =0 . 1  SIG2=0 . 7  RHO=0 . 8 


Figure  6-100.  Elevator  Position 
Configuration  8/Operating  Point:  38K  0.9M 
SIG1 =0 . 1  SIG2=0 . 7  RHO=0 . 8 
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Figure  6-101.  Flaperon  Rate 
Configuration  8/Operating  Point:  38K  0.9M 
SIG1=0.1  SIG2=0 . 7  RHO=0.8 
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Figure  6-102.  Elevator  Rate 
Configuration  8/Operating  Point:  38K  0 . 9M 
SIG1 =0 . 1  SIG2=0 . 7  RHO=0 . 8 


Table  6-15 


Sensor  Noise  Data 
(units  are  degrees) 


Measured  Quantity 

Noise  Mean 

Variance 

q 

0.0 

7.438E-8 

© 

0.0 

4 . 886E-7 

a 

o 

• 

o 

3.404E-6 

the  bank  of  secondary  estimators.  As  a  baseline,  a  simula¬ 
tion  without  sensor  noise  was  run.  The  weights  associated 
with  each  candidate  model  as  well  as  the  prediction  error 
variances  are  shown  in  Figures  6-103  through  6-107.  From 
Figure  6-103,  it  can  be  seen  that  the  probabilty  of  the 
model  associated  with  the  10,000  ft,  0.35  Mach  flight  con- 

Table  6-16 

Two-Model  Configuration  for  Noise  Consideration 


Model 

Flight  Condition 

1 

10,000  ft 

0.35  Mach 

2 

10,000  ft 

0.50  Mach 
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Figure  6-103.  Model  Probability  Weightings 
Operating  Point:  10K  .35M 
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Figure  6-104.  Prediction  Error  Variance  (1,1)  for  Model  1 

Operating  Point:  10K  0.35M 
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Figure  6-105.  Prediction  Error  Variance  (2,2)  for  Model  1 

Operating  Point:  10K  0.35M 
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Figure  6-106 
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Prediction  Error  Variance  (1,1)  for  Model  2 
Operating  Point:  10K  0.35M 
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Figure  6-107 


Prediction  Error  Variance  (2,2)  for  Model  2 
Operating  Point:  10K  0.35M 
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figuration  is  very  high.  In  comparing  the  prediction  error 
variances  for  model  1  and  model  2,  Figures  6-104  through 
6-107  show  that  the  variances  associated  with  model  1  are 
smaller,  indicating  that  it  is  the  better  of  the  two  models. 

Adding  the  sensor  noise  levels  as  given  in  Table  6-15 
and  running  a  simulation  for  the  configuration  in  Table  6-16 
resulted  in  the  data  presented  on  Figures  6-108  through 
6-112.  The  performance  of  the  multiple  model  algorithm  as 
shown  on  the  probability  curve  of  Figure  6-108  is  very 
similar  to  that  of  Figure  6-103  (same  conditions  with  no 
noise).  The  prediction  error  variances  on  Figures  6-109 
through  6-112  have  responses  that  are  very  similar  to  the  no 
noise  case,  but  are  of  increased  magnitude.  Figures  6-113 
through  6-117  show  a  simulation  that  was  performed  with  the 
sensor  noise  figures  set  at  one-hundred  times  those  shown  in 
Table  6-15.  The  probability  weighting  curve  shown  on  Figure 
6-113  shows  the  effect  of  adding  too  much  sensor  noise  on 
the  performance  of  the  multiple  model  algorithm.  Instead  of 
model  1  being  selected  as  the  correct  model  with  a  probabil¬ 
ity  of  almost  one,  a  substantial  contribution  of  model  2  is 
now  present.  The  reason  for  this  can  be  seen  by  comparing 
the  prediction  error  variances  as  shown  in  Figures  6-114 
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•  Figure  6-108.  Model  Probability  Weightings 

Sensor  Noise/Operating  Point:  10K  0.35M 
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Figure  6-109.  Prediction  Error  Variance  (1,1)  for  Model  1 
Sensor  Noise/Operating  Point:  10K  0.35M 


Figure  6-110.  Prediction  Error  Variance  (2,2)  for  Model  1 
Sensor  Noise/Operating  Ponit:  10K  0.35M 
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Figure  6-111.  Prediction  Error  Variance  (1,1)  for  Model 
Sensor  Noi ae/Operat ing  Point:  10K  0.35M 
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Figure  6-112.  Prediction  Error  Variance  (2,2)  for  Model  2 
Sensor  Noise/Operating  Point:  10K  0.35M 
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Figure  6-113.  Model  Probability  Weightings 
Sensor  Noise (* 100 ) /Operating  Point:  10K  0.35M 
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Figure  6-114.  Prediction  Error  Variance  (1,1)  for  Model  1 
Sensor  Noi se (* 1 00 ) /Operating  Point:  10K  0.35M 


Figure  6-115.  Prediction  Error  Variance  (2,2)  for  Model  1 
Sensor  Noise {* 100 ) /Operating  Point:  10K  0.35M 
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Figure  6-116.  Prediction  Error  Variance  (1,1)  for  Model  2 
Sensor  Noise (* 100 ) /Operating  Point:  10K  0.35M 
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Figure  6-117.  Prediction  Error  Variance  (2,2)  for  Model  2 
Sensor  Noise (* 1 00 ) /Opera ting  Point:  10K  0.35M 


through  6-117.  By  comparing  Figures  6-114  and  6-116  (pre¬ 


diction  error  variance  for  pitch  rate  for  models  1  and  2 
respectively),  one  can  see  that  the  variances  are  very 
similar  for  the  time  of  7  through  10  seconds.  Comparing 
Figure  6-115  to  6-117  (prediction  error  variance  for  flight 
path  for  models  1  and  2,  respectively),  reveals  the  same 
effect.  This  effect  can  be  attributed  to  the  sensor  noise 
being  added  to  the  system.  The  noise  is  actually  masking 
the  residuals,  thereby  leading  to  a  similarity  in  prediction 
error  variances  and  hence  to  degraded  multiple  model  algo¬ 
rithm  performance.  Comparison  of  the  prediction  error  var¬ 
iances  of  the  no-noise  case  of  Figures  6-104  through  6-107 
with  those  of  Figures  6-114  through  6-117  shows  clearly  that 
the  addition  of  sensor  noise  can  mask  the  residuals,  and  if 
sufficiently  high  noise  is  added,  the  prediction  error  var¬ 
iances  will  all  have  the  same  relative  responses.  This 
similarity  of  response  will  result  in  degraded  performance 
by  the  multiple  model  algorithm  in  assigning  the  weighting 
factors  to  the  proper  models. 

The  second  two-model  configuration  evaluated  is  shown 
in  Table  6-17.  The  two  models  were  spaced  farther  apart 


than  those  of  the  two-model  configuration  presented  in  Table 


Table  6-17 


Two  Model  Configuration  for  Noise  Consideration 


Model 

Flight  Condition 

1 

10,000  ft 

0.35  Mach 

2 

10,000  ft 

0.75  Mach 

6-16  to  ascertain  the  effects  of  model  spacing  in  relation 
to  noise  resistance.  Figures  6-118  through  6-122  present 
simulation  data  for  the  two-model  configuration  with  the 
actual  flight  condition  set  at  10,000  ft,  0.35  Mach.  Figure 
6-118  shows  a  very  high  probability  associated  with  model  1, 
which  is  the  model  associated  with  the  10,000  ft,  0.35  Mach 
flight  condition.  Allowing  the  sensor  noise  to  take  on 
values  one-hundred  times  those  presented  in  Table  6-15  and 
running  a  simulation  resulted  in  the  data  presented  on 
Figures  6-123  through  6-127.  The  data  shows  that  the  proba¬ 
bility  weighting  curve  (Figure  6-123)  indicates  that  model  1 
is  the  correct  model,  as  was  the  case  with  no  noise.  A 
comparison  of  the  prediction  error  variances  for  models  1 
and  2  (Figures  6-124  through  6-127)  shows  that  the  variances 
associated  with  model  1  are  smaller  than  those  of  model  2. 
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—  Model  at  10K  0.35M 
--  Model  at  10K  0.75M 
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Figure  6-118.  Model  Probability  Weightings 
Operating  Point:  10K  0.35M 
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Operating  Point:  10K  0.35M 
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Figure  6-120 
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Prediction  Error  Variance  (2,2)  for  Model  1 
Operating  Point:  10K  0.35M 
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Prediction  Error  Variance  (1,1)  for  Model  2 
Operating  Point:  10K  0.35M 
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Figure  6-123. 


Model  Probability  Weightings 


Sensor  Noise (* 100 ) /Operating  Point:  10K  0.35M 
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Figure  6-124.  Prediction  Error  Variance  (1,1)  for  Model  1 
Sensor  Noise( *  100 ) /Operating  Point:  10K  0.35M 
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Figure  6-125.  Prediction  Error  Variance  (2,2)  for  Model  1 
Sensor  Noise {* 100 ) /Operating  Point:  10K  0.35M 
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Figure  6-126.  Prediction  Error  Variance  (1,1)  for  Model  2 
Sensor  Noi se (* 1 00 ) /Operating  Point:  10K  0.35M 


Figure  6-127.  Prediction  Error  Variance  (2,2)  for  Model  2 
Sensor  Noise (* 100 ) /Operating  Point:  10K  0.35M 
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Even  with  the  sensor  noise  figures  increased  by  a  factor  of 
one-hundred,  the  multiple  model  algorithm’s  performance  has 
not  been  severely  degraded,  and  the  correct  model  has  been 
selected.  This  is  in  contrast  to  the  results  obtained  for 
the  two-model  configuration  of  Table  6-16.  For  the  config¬ 
uration  of  Table  6-16,  an  increase  in  the  sensor  noise 
figures  of  one-hundred  times  resulted  in  degraded  multiple 
model  algorithm  performance.  The  reason  for  the  performance 
difference  between  the  two-model  configurations  of  Tables  6- 
16  and  6-17  is  due  to  the  spacing  of  the  models.  The  models 
of  Table  6-13  are  spaced  farther  apart  (in  dynamic  pressure) 
than  those  of  Table  6-16  and  therefore  the  difference  of  the 
two  models  as  seen  in  the  residuals  is  more  significant. 
Therefore,  the  error  variance  term  introduced  by  a  model  not 
being  the  best  fitting  model  (see  Equation  (5-9))  is  still 
the  dominant  term  and  a  difference  in  the  prediction  error 
variance  responses  is  the  result.  This  clearly  points  out 
the  fact  that,  the  farther  apart  the  models  of  the  multiple 
model  algorithm  are  placed  (from  a  dynamic  pressure  stand¬ 
point),  the  more  noise  resistance  is  possible. 

Following  the  simulations  for  the  two-model  configura¬ 
tions,  simulations  using  configuration  8  (see  Table  6-9) 


were  run.  Recall  that  configuration  8  was  the  three-model 
configuration  that  gave  the  best  performance  from  a  perfor¬ 
mance  criteria  standpoint.  The  purpose  of  these  simulations 
was  to  see  if  the  performance  of  three-model  configuration 
was  severely  degraded  in  the  presence  of  realistic  sensor 
noise.  Several  flight  conditions  were  selected  as  operating 
points,  the  first  of  which  was  26,000  ft,  and  0.9  Mach.  The 
data  presented  on  Figures  6-128  through  6-134  are  for  the 
case  where  no  sensor  noise  was  added.  The  probability 
weighting  curve  (see  Figure  6-128)  shows  that  models  2  and  3 
contribute  significantly  to  the  solution  of  finding  the 
correct  model.  Figures  6-129  through  6-134  show  the  predic¬ 
tion  error  variances  for  the  respective  models.  Figure  6- 
135  shows  the  probability  weighting  curve  for  the  sensor 
noise  values  presented  in  Table  6-15.  From  this  figure,  it 
can  be  seen  that  although  models  2  and  3  are  still  the  most 
heavily  weighted  models,  model  3  has  more  weight  than  it  did 
in  the  simulation  without  noise.  Further  increasing  the 
sensor  noise  to  one-hundred  times  the  values  in  Table  6-15 
results  in  the  data  on  Figures  6-136  through  6-142.  The 
probability  weighting  curve  (see  Figure  6-136)  is  now  dras¬ 
tically  different  from  the  case  with  no  noise,  and  model  3 


•  Figure  6-128.  Model  Probability  Weightings 

Configuration  8/Operating  Point:  26K  0.9M 
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Figure 


6-129.  Prediction  Error  Variance  (1,1)  for  Model  1 
Configuration  8/Operating  Point:  26K  0.9M 
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Figure  6-130.  Prediction  Error  Variance  (2,2)  for  Model  1 
Configuration  8/Operating  Point:  26K  0 . 9M 
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Figure  6-133.  Prediction  Error  Variance  (1,1)  for  Model  3 
Configuration  8/Operating  Point:  26K  0 . 9M 
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Figure  6-134.  Prediction  Error  Variance  (2,2)  for  Model  3 
Configuration  8/Operating  Pr-int:  26K  0.9M 
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Figure  6-135.  Model  Probability  Weightings 
onfiguration  8/Sensor  Noise/Operating  Point:  26K 
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Figure  6-137.  Prediction  Error  Variance  (1,1)  for  Model  1 
Configuration  8/Sensor  Noise (* 100 ) /Operating  Point:  26K  0 . 9M 
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Figure  6-138.  Prediction  Error  Variance  (2,2)  for  Model  1 
Configuration  8/Sensor  Noise (* 100 ) /Operating  Point:  26K  0 . 9M 
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Figure  6-139.  Prediction  Error  Variance  (1,1)  for  Model  2 
Configuration  8/Sensor  Noise (* 100 ) /Operating  Point:  26K  0 . 9M 
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Figure  6-140.  Prediction  Error  Variance  (2,2)  for  Model  2 
Configuration  8/Sensor  Noise (* 100 ) /Operating  Point:  26K  0.9M 
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Figure  6-141.  Prediction  Error  Variance  (1,1)  for  Model  3 
Configuration  8/Sensor  Noise (* 1 00 ) /Operat ing  Point:  26K  0.9M 


has  been  selected  as  the  most  correct  model.  The  prediction 
error  variances  as  shown  in  Figures  6-136  through  6-142  also 
have  much  different  responses  from  those  that  were  obtained 
from  the  simulations  without  noise.  Several  other  simula¬ 
tions  were  run  with  the  same  general  results.  Figures  6-143 
through  6-146  present  the  probability  weighting  curves  for 
increasing  values  of  sensor  noise  for  the  case  where  the 
actual  flight  condition  was  at  38,000  ft,  0.6  Mach.  As  can 
be  seen  from  the  weighting  curves  of  Figures  6-143  through 
6-146,  increasing  the  level  of  sensor  noise  results  in 
increasing  degradation  of  the  multiple  model  algorithm’s 
performance.  From  the  numerous  simulations  accomplished,  it 
was  determined  that  the  multiple  model  algorithm  was  more 
sensitive  to  sensor  noise  when  the  probability  weighting 
curves  for  no  noise  had  two  models  with  significant  influ¬ 
ence  on  determining  the  correct  model.  In  contrast,  the 
multiple  m  ^del  algorithm  is  less  sensitive  to  sensor  noise 
when  the  actual  flight  condition  is  such  that  the  probabili¬ 
ty  weighting  calculation  for  no  noise  yields  a  single  model 
with  a  probability  near  one.  The  simulation  results  for 
configuration  8  showed  that  the  multiple  model  algorithm 
yielded  satisfactory  results  for  the  sensor  noise  levels 
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Figure  6-146.  Model  Probability  Weightings 
Configuration  8/Sensor  Noi se (* 100 ) /Operating  Point:  38K  0.6M 


given  in  Table  6-15.  However,  increasing  the  noise  levels 


beyond  those  in  Table  6-15  resulted  in  serious  performance 
degradation  of  the  multiple  model  algorithm. 

6^3 _ Summa  ry 

This  chapter  presents  the  results  that  were  obtained 
during  the  accomplishment  of  this  research  effort.  Using 
several  nominal  models  as  data  points,  performance  bound¬ 
aries  are  generated  in  an  effort  to  determine  regions  of 
robustness  for  the  given  models,  and  observations  about  the 
primary  factor  in  determining  model  equivalence  are  made. 

The  effects  of  varying  the  flight  condition  on  the  control 
surface  performance  is  also  evaluated.  Next,  simulations 
with  the  multiple  model  algorithm  are  accomplished.  At 
first,  two  models  are  placed  in  the  secondary  estimator 
bank,  to  determine  the  required  amount  of  overlap  of  the 
performance  boundaries  of  the  elemental  models  to  achieve 
the  desired  tracking  performance.  Following  the  evaluation 
of  the  two-model  configurations,  several  three-model  config¬ 
urations  are  evaluated  and  are  f ound  to  to  yield  satisfacto¬ 
ry  tracking  performance  over  the  desired  flight  envelope. 

In  an  effort  to  minimize  the  amount  of  time  the  control 
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surfaces  reached  the  rate  limits,  the  design  parameters  of 


the  controller  are  adjusted.  Finally,  the  effects  of  sensor 
noise  on  the  performance  of  the  multiple  model  algorithm  are 
evaluated.  The  next  chapter  of  this  thesis  presents  the 
conclusions  drawn  from  the  results  detailed  in  this  chapter. 
Recommendations  for  further  study  are  also  presented. 


7.  Conclusions  and  Recommendations 

7.1  Introduction 

As  mentioned  several  times  throughout  this  thesis,  the 
primary  objective  of  this  research  is  to  determine  the 
effects  of  model  placement  on  the  performance  of  the  multi¬ 
ple  model  algorithm  within  an  adaptive  control  law.  Several 
hundred  simulations  were  run  which  are  the  basis  for  the 
results  presented  in  Chapter  6.  The  following  list  of 
conclusions  presents  the  highlights  of  the  results  in  Chap¬ 
ter  6.  The  recommendations  outline  further  areas  of  study 
in  this  area  of  research. 

7 . 2  Conclusions 

The  following  is  a  list  of  the  highlights  of  the  re¬ 
sults  obtained  in  the  course  of  this  research  effort. 

1,  Performance  criteria  must  be  selected  prior  to 
determining  the  performance  robustness  of  a  nominal  aircraft 
model.  The  performance  criteria  dictate  the  degree  to  which 
the  plant  parameters  can  vary  before  unacceptable  perfor¬ 
mance  is  reached  and  therefore  influence  the  performance 
boundaries  for  a  given  nominal  model. 

2.  In  determining  the  performance  boundaries  for  sev- 
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eral  different  nominal  models,  it  has  been  learned  that  the 
quantity  model  equivalence  is  most  dependent  on  is  dynamic 
pressure.  By  taking  several  different  models  on  approxi¬ 
mately  the  same  dynamic  pressure  line,  it  is  shown  that  the 
performance  boundaries  are  the  same. 

3.  Although  the  performance  boundaries  for  nominal 
models  on  the  same  dynamic  pressure  line  are  the  same,  the 
control  surface  responses  may  be  different.  It  is  shown 
that  the  reason  for  this  is  the  different  trim  values  of 
angle-of-attack .  It  is  also  shown  that  control  surface 
performance  is  influenced  by  the  aerodynamic  gains  as¬ 
sociated  with  the  actual  flight  condition.  By  decreasing 
altitude  and/or  increasing  Mach  number,  the  aerodynamic  gain 
increases,  causing  the  control  surfaces  to  respond  faster. 
This  information  is  useful  when  tuning  the  controller.  Af¬ 
ter  selection  of  the  models  for  the  multiple  model  al¬ 
gorithm,  the  controller  can  be  tuned  at  the  flight  condition 
with  the  worst  case  aerodynamic  gains. 

4.  When  trying  to  achieve  satisfactory  tracking  per¬ 
formance  over  the  flight  envelope  of  interest  with  two 
models  in  the  multiple  model  algorithm’s  secondary  bank,  it 
is  seen  that  the  amount  of  overlap  of  the  performance  bound- 


aries  is  very  important.  Due  to  the  fact  that  the  multiple 
model  adaptive  estimator  arrives  at  a  parameter  estimate  by 
weighting  the  "correctness"  of  the  models  in  the  bank,  the 
final  parameter  estimate  is  a  blend  of  those  models’  parame¬ 
ters.  Because  of  the  weighted  blending  of  parameters  of 
each  of  the  models  in  the  bank,  a  significant  amount  of 
overlap  of  the  performance  boundaries  is  required  to  yield 
satisfactory  tracking  performance.  However,  it  is  also 
shown  that  the  lower  bound  on  the  computed  probabilities 
affects  the  amount  of  overlap  required  of  the  performance 
boundaries.  The  amount  of  overlap  increases  as  the  value  of 
the  lower  bound  becomes  larger.  This  performance  boundary 
overlap  requirement  leads  to  the  use  of  three  models  as 
secondary  parameter  estimators  in  lieu  of  two. 

5.  Due  to  the  robustness  of  the  controller  to  varia¬ 
tions  in  the  plant  parameters,  the  flight  envelope  of  inter¬ 
est  can  be  modelled  with  three  nominal  models  in  the  bank  of 
secondary  estimators  in  the  multiple  model  algorithm.  The 
robustness  of  the  control  law  allows  a  large  parameter  space 
to  be  modelled  by  a  relatively  small  number  of  discrete 
models  while  yielding  excellent  tracking  performance. 

6.  In  the  analysis  of  the  three-model  configurations, 


it  is  shown  that,  regardless  of  the  altitudes  or  Mach  num¬ 
bers  of  the  nominal  models  in  the  secondary  estimator  bank 
of  the  multiple  model  algorithm,  the  model  assigned  the 
highest  weighting  is  the  model  with  a  value  of  dynamic 
pressure  closest  to  that  of  the  actual  flight  condition. 

7.  The  limitation  on  how  close  the  candidate  models 
can  be  placed  is  shown  to  be  dependent  on  the  amount  of 
senior  .ioise  added  to  the  longitudinal  states.  Sensor  noise 
tends  to  mask  the  residuals  generated  for  each  of  the  models 
in  the  secondary  estimator  bank  and  hence  degrades  the 
multiple  model  algorithm’s  performance.  It  is  shown  that 
with  greater  separation  between  dynamic  pressure  parameteri¬ 
zation  the  susceptabi 1 i ty  of  the  multiple  model  algorithm  to 
the  sensor  noise  is  decreased. 

From  the  above  list  of  conclusions,  it  can  be  seen  that 
model  placement  for  the  multiple  model  algorithm  is  depen¬ 
dent  on  many  factors.  First  of  all,  performance  criteria 
must  be  selected,  which  in  turn  determine  the  region  of 
performance  robustness  of  a  given  nominal  model.  The  perfor¬ 
mance  criteria  also  have  a  large  effect  on  the  number  of 
models  (for  the  multiple  model  algorithm's  secondary  estima¬ 
tion  bank)  that  are  necessary  to  achieve  the  desired  track- 


ing  performance  over  the  specified  flight  envelope.  Less 
demanding  performance  criteria  would  result  in  larger  re¬ 
gions  of  performance  robustness  and  hence  fewer  models  would 
be  required  to  model  the  desired  parameter  space. 

The  amount  of  overlap  of  the  performance  boundaries  is 
shown  to  require  a  tradeoff  of  conflicting  requirements. 
These  conflicting  requirements  are  the  result  of  wanting  a 
large  amount  of  overlap  for  enhanced  tracking  performance, 
yet  desiring  minimal  overlap  for  sensor  noise  resistance. 

The  results  show  that  the  minimum  amount  of  overlap  that  is 
acceptable  is  that  which  yields  satisfactory  tracking  per¬ 
formance  over  the  flight  envelope  of  interest.  The  maximum 
amount  of  overlap  is  dictated  by  the  expected  levels  of 
sensor  noise.  Therefore,  once  performance  criteria  have 
been  established  and  sensor  noise  figures  determined,  an 
analysis  as  to  how  far  apart  in  dynamic  pressure  the  secon¬ 
dary  estimator  models  need  to  be  can  be  accomplished. 

7.3  Recommendations  for  Further  Study 

This  thesis  presents  an  initial  evaluation  of  the  fac¬ 
tors  that  are  important  in  the  selection  of  models  for  the 
multiple  model  algorithm.  The  following  list  presents  areas 


of  interest  that  could  be  investigated  to  further  this 
ef  f  ort . 

1.  The  flight  envelope  of  interest  should  be  expanded 
to  include  transonic  and  supersonic  flight  conditions.  The 
purpose  of  investigating  this  expanded  flight  envelope  would 
be  to  verify  if  model  equivalence  is  based  primarily  on 
dynamic  pressure  or  on  some  other  criteria.  Model  selection 
in  this  portion  of  the  flight  envelope  might  have  different 
characteristics  from  those  addressed  in  this  thesis. 

2.  The  multiple  model  algorithm  could  be  tailored  so 
that  the  gains  implemented  in  the  controller  would  not  be  a 
weighted  gain  matrix  of  the  models  in  the  bank,  but  would  be 
the  set  of  gains  associated  with  the  model  with  the  highest 
probability.  This  would  seem  to  have  the  effect  of  not 
requiring  much  overlap  of  the  performance  boundaries  and 
therefore  would  provide  greater  sensor  noise  resistance. 

3.  Although  the  general  effects  of  sensor  noise  on 
model  selection  is  addressed  in  this  effort,  further  re¬ 
search  in  this  area  needs  to  be  accomplished  to  include  an 
extensive  Monte  Carlo  analysis. 


Appendix  A 


This  appendix  presents  the  A  and  B  matrices  of  the 
state  space  models  for  each  of  the  flight  conditions  that 
are  considered  in  the  flight  envelope  of  interest.  In 
addition,  the  aircraft  parameters  that  are  constant  for  each 
flight  condition  are  presented.  This  data  was  obtained  from 
the  Flight  Dynamics  Laboratory  at  Wright-Patterson  AFB ,  OH. 


Constant  Aircraft  Parameters 

S  (wing  reference  area  -  ft~2)  =  300.00 
c  (wing  mean  aerodynamic  chord  -  ft)  =  11.32 
ACGW  (aircraft  gross  weight  -  lbs)  =  21505.00 
LOADFAC  (load  factor  -  g)  =  1.0 
b  (wing  span  -  ft)  =  30.00 

Inertias: 

Ixx  (slug  -  ft~2)  =  14224.04 
Iyy  (slug  -  ft~2)  =  61352.40 
Izz  (slug  -  ft'2)  =  72651.75 


Flight  Condition:  10,000  ft,  0 

0.0000  0.0000  0.0000 

A  =  -32.1640  -0.0154  43.4875 

-0.0008  0.0000  -2.0397 

0.0003  -0.0005  6.9843 

Flight  Condition:  10,000  ft,  0 

0.0000  0.0000  0.0000 

A  =  -32.1625  -0.0156  39.7461 

-0.0009  0.0000  -1.8507 

0.0003  -0.0010  7.1505 

Flight  Condition:  10,000  ft,  0 

0.0000  0.0000  0.0000 

A  =  -32.1603  -0.0149  35.8901 

-0.0011  -0.0001  -1.6730 

0.0003  -0.0010  7.1562 

Flight  Condition:  10,000  ft,  0 

0.0000  0.0000  0.0000 

A  =  -32.1584  -0.0138  33.9936 

-0.0012  -0.0001  -1.5069 

0.0003  -0.0010  6.8526 

Flight  Condition:  10.000  ft,  0 

0.0000  0.0000  0.0000 

A  =  -32.1575  -0.0118  30.9741 

-0.0014  -0.0001  -1.3407 

0.0003  -0.0010  6.5776 

Flight  Condition:  10,000  ft,  0 

0.0000  0.0000  0.0000 

A  =  -32.1550  -0.0112  30.6762 

-0.0016  -0.0001  -1.1975 

0.0003  -0.0010  6.1077 


9  Mach 

1.0000  0.0000  0.0000 
-24.0510  B  =  -0.9568  19.0164 
0.9999  -0.2063  -0.3610 

-0.9769  -30.8710  -9.4022 

85  Mach 

1.0000  0.0000  0.0000 
-24.4587  B  =  -0.4927  15.6819 
1.0000  -0.3380  -0.1833 

-0.9116  -26.8660  -7.5573 

8  Mach 

1.0000  0.0000  0.0000 
-25.1231  B  =  -0.2750  12.8572 
1.0000  -0.1833  -0.3151 

-0.8477  -23.1933  -5.9530 

75  Mach 

1.0000  0.0000  0.0000 

-25.1510  B  =  0.6016  10.4851 

1.0000  -0.1662  -0.2979 

-0.7854  -19.8129  -4.6123 

7  Mach 

1.0000  0.0000  0.0000 

-24.1602  B  =  2.2460  3.8331 

1.0000  -0.1547  -0.2636 

-0.7244  -16.7132  -3.1283 

65  Mach 

1.0000  0.0000  0.0000 

-24.0485  B  =  3.5408  -5.0936 

1.0000  -0.1432  -0.2406 

-0.6648  -13.9687  -1.6329 


e 


Flight  Condition:  10,000  ft,  0.6  Mach 


* 


© 


* 


0.0000 

0.0000 

0.0000  1.0000 

0.0000 

i 

A  =  -32. 1506 

-0.0108 

32.2231  -24.6667 

B 

4.3316 

-i 

-0.0019 

-0.0001 

-1.0657  1.0000 

-0.1317 

- 

0.0003 

-0.0004 

5.4030  -0.6067 

- 

11 .6081 

- 

Flight  Condition:  10,000  ft,  0.55  Mach 

0.0000 

0.0000 

0.0000  1.0000 

0.0000 

A  =  -32.1429 

-0.0099 

32.7156  -26.0481 

B 

4.4748 

- 

-0.0024 

-0.0001 

-0.9855  1.0000 

-0.1261 

- 

0.0003 

0.0002 

4.3259  -0.5577 

-9.7689 

- 

Flight  Condition:  10,000  ft,  0.5  Mach 

0.0000 

0.0000 

0.0000  1.0000 

0.0000 

A  =  -32.  1305 

-0.0102 

32.2804  -28.0069 

B 

= 

4.3889 

- 

-0.0031 

-0.0001 

-0.8995  1.0000 

-0.1146 

- 

0.0004 

0.0002 

3.5523  -0.5124 

-8.0844 

- 

Flight  Condition:  10,  i 

000  ft,  0.45  Mach 

0.0000 

0.0000 

0.0000  1.0000 

0.0000 

A  =  -32  .  1092 

-0.0100 

31.8221  -30.7710 

B 

= 

3.7471 

- 

-0.0042 

-0.0001 

-0.8079  1.0000 

-0.1089 

- 

0.0004 

0.0001 

3.0825  -0.4675 

-6.5832 

- 

Flight  Condition:  10,000  ft,  0.4  Mach 

0.0000 

0.0000 

0.0000  1.0000 

0.0000 

A  =  -32.0708 

-0.0102 

27.0092  -34.5258 

B 

= 

2.6872 

- 

-0.0060 

-0.0002 

-0.6933  1.0000 

-0.0974 

- 

0.0005 

0.0001 

2.7158  -0.4228 

-5.1967 

- 

Flight  Condition:  10, 

000  ft,  0.35  Mach 

0.0000 

0.0000 

0.0000  1.0000 

0.0000 

A  -  -31.9931 

-0.0109 

19.7499  -40.0088 

B 

1.9309 

- 

-0.0090 

-0.0002 

-0.5672  1.0000 

-0.0859 

- 

0.0006 

0.0001 

2.0512  -0.3905 

-3.9362 

% 


& 


0.0000 

6.2338 

0.2349 

0.9454 


0.0000 

5.0133 

0.2120 

0.7850 


0.0000 
4.  1826 
0.1891 
0.5730 


0.0000 

3.3174 

0.1662 

0.3724 


0.0000 
2.1371 
0. 1432 
0.1833 


0.0000 
3.3690 
0. 1031 
0.0344 
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e 

Flight  Condition:  10, 

000  ft,  0.3  Mach 

0.0000 

0.0000 

0.0000  1.0000 

0.0000  0.0000 

• 

A  =  -31.8063 

-0.0127 

19.4232  -48.9139 

B 

=  1.6157  -3.3289 

-0.0150 

-0.0003 

-0.4526  1.0000 

-0.0745  -0.0630 

0.0008 

0.0001 

1.8793  -0.3862 

-2.9794  0.3724 

• 

Flight  Condition:  10, 

000  ft,  0.25  Mach 

0.0000 

0.0000 

0.0000  1.0000 

0.0000  0.0000 

A  =  -31.2389 

0.0075 

11.6253  -65.0974 

B 

=  -0.2005  -2.1257 

-0.0286 

-0.0004 

-3.5626  1.0000 

-0.0630  -0.0573 

• 

0.0003 

0.0001 

0.6131  -0.4665 

-2.0684  0.2349 

Flight  Condition:  14,000  ft,  0.9  Mach 

0.0000 

0.0000 

0.0000  1.0000 

0.0000  0.0000 

A  =  -32. 1621 

-0.0129 

37.7866  -25.9832 

B 

=  -0.9397  16.5356 

* 

-0.0009 

0.0000 

-1.7819  1.0000 

-0.1776  -0.3266 

0.0003 

-0.0005 

5.6207  -0.8498 

-26.6024  -8.1932 

Flight  Condition:  14, 

000  ft,  0.85  Mach 

« 

0.0000 

0.0000 

0.0000  1.0000 

0.0000  0.0000 

A  =  -32 .  1599 

-0.0132 

34.5150  -26.6752 

B 

=  -0.6303  13.6708 

-0.0011 

-0.0001 

-1.6157  1.0000 

-0.1662  -0.3037 

0.0003 

-0.0009 

5.8270  -0.7927 

-23.1246  -6.5776 

* 

Flight  Condition:  14, 

000  ft,  0.8  Mach 

0.0000 

0.0000 

0.0000  1.0000 

0.0000  0.0000 

A  =  -32.1569 

-0.0126 

34.6410  -27.6587 

B 

=  -0.5386  11.2529 

-0.0012 

-0.0001 

-1.4610  1.0000 

-0.1604  -0.2807 

<• 

0.0003 

-0.0009 

5.7353  -0.7369 

-19.9447  -5.1738 

Flight  Condition:  14, 

000  ft,  0.75  Mach 

0.0000 

0.0000 

0.0000  1.0000 

0.0000  0.0000 

A  =  -32 . 1549 

-0.0115 

32.4638  -27.4111 

B 

=  6.4171  9.2246 

-0.0014 

-0.0001 

-1.3178  1.0000 

-0.1489  -0.2636 

0.0003 

-0.0009 

5.4947  -0.6826 

-16.9997  -4.0164 

« 
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m 

Flight  Condition:  14,000  ft,  0.7  Mach 


0.0000  0.0000  0.0000 

A  =  -32.1532  -0.0099  29.7479 

-0.0016  -0.0001  -1.1745 

0.0003  0.0010  5.2196 

Flight  Condition:  14,000  ft,  0 

0.0000  0.0000  0.0000 

A  =  -32.1490  -0.0102  30.6876 

-0.0018  -0.0001  -1.0509 

0.0003  -0.0010  4.8210 

Flight  Condition:  14,000  ft,  0 

0.0000  0.0000  0.0000 

A  =  -32.1426  -0.0098  31.8444 

-0.0022  -0.0001  -0.9408 

0.0003  -0.0004  4.2991 

Flight  Condition:  14,000  ft,  0 

0.0000  0.0000  0.0000 

A  =  -32.1322  -0.0087  31.8936 

-0.0028  -0.0001  -0.8662 
0.0003  0.0002  3.5602 

Flight  Condition:  14,000  ft,  0 

0.0000  0.0000  0.0000 

A  =  -32.1151  -0.0088  31.5910 

-0.0037  -0.0001  -0.7893 

0.0004  0.0001  3.0968 

Flight  Condition:  14,000  ft,  0 

0.0000  0.0000  0.0000 

A  =  -32.0859  -0.0085  28.8189 

-0.0050  -0.0002  -0.6937 

0.0004  0.0001  2.8036 


1.0000  0.0000  0.0000 

-26.7275  B  =  1.9881  0.7391 

1 .0000  -0. 1375  -0.2291 

-0.6294  -14.3182  -2.4981 

65  Mach 

1.0000  0.0000  0.0000 

-27.1841  B  =  3.2771  -6.7577 

1.0000  -0.1255  -0.2087 

-0.5777  -11.9849  -1.2483 

6  Mach 

1.0000  0.0000  0.0000 

-28.1420  B  =  3.8231  -4.9670 

1.0000  -0.1164  -0.2058 

-0.5277  -9.9623  -0.8613 

55  Mach 

1.0000  0.0000  0.0000 

-29.7446  B  =  3.9308  -4.2254 

1.0000  -0.1089  -0.1853 

-0.4883  -8.3808  -0.6284 

5  Mach 

1.0000  0.0000  0.0000 

-32.1056  B  =  3.8867  -3.4343 

1.0000  -0.1007  -0.1649 

-0.4492  -6.9347  -0.4166 

45  Mach 

1.0000  0.0000  0.0000 

-35.3574  B  =  3.2049  -2.1644 

1.0000  -0.0926  -0.1440 

-0.4105  -5.6439  -0.1982 


Flight  Condition:  14,000  ft,  0.4  Mach 


A 


0.0000 

-32.0320 

-0.0071 

0.0005 


0.0000 

-0.0090 

-0.0002 

0.0001 


0.0000 

22.3943 

-0.5802 

2.3182 


1 .0000 
-39.8968 
1 .0000 
-0.3763 


Flight  Condition:  14,000  ft,  0.35  Mach 


A 


0.0000 
-31 .9175 
-0.0109 
0.0006 


0.0000 

-0.0104 

-0.0002 

0.0001 


0.0000 

17.6824 

-0.4728 

1.8844 


1.0000 
-46 . 9432 
1.0000 
-0.3560 


Flight  Condition:  14,000  ft,  0.3  Mach 


A 


0.0000 

-31.6421 

-0.0183 

0.0007 


0.0000 

-0.0061 

-0.0003 

0.0001 


0.0000 
17 . 8659 
-0.3965 
1 . 3232 


1.0000 

-57.9781 

1.0000 

-0.3805 


Flight  Condition:  14,000  ft,  Q.25  Mach 


A 


0.0000 

-30.8164 

-0.0349 

-0.0004 


0.0000 

0.0196 

-0.0004 

0.0000 


0.0000 

8.2127 

-0.3112 

0.1568 


1.0000 
-77 .3577 
1.0000 
-0.4605 


Flight  Condition:  18,000  ft,  0.9  Mach 


A 


0.0000 

-32.1594 

-0.0010 

0.0003 


0.0000 

-0.0107 

0.0000 

-0.0004 


0.0000 
33. 1373 
-1 . 5498 
4 .4287 


1.0000 
-28 . 3554 
1.0000 
-0.7359 


Flight  Condition:  18,000  ft,  0.85  Mach 


A 


0.0000 

-32.1564 

-0.0012 

0.0003 


0.0000 

-0.0110 

-0.0001 

-0.0009 


0.0000 
33.7813 
-1 . 4085 
4.6260 


1.0000 
-29 . 3747 
1.0000 
-0.6862 


0.0000 

B  =  2.2655  - 

-0.0834  - 
-4.4398  - 


0.0000 

B  =  1.9270  - 

-0.0730  - 
-3.4005 


0.0000 

B  =  0.6524  - 

-0.0632  - 
-2.5582 


0.0000 

B  =  -0.9529  - 
-0.0558  - 
-1.7615 


0.0000 

B  =  -0.9453  1 
-0.1563  - 
-22.8082  - 


0.0000 

B  =  -0.7707  1 
-0.1469  - 
-19.8074  - 


0.0000 
1 .8722 
0. 1248 
0. 1450 


0.0000 

4.4673 

0.0672 

0.2782 


0.0000 

2.4686 

0.0569 

0.3444 


0.0000 

1.4814 

0.0425 

0.1750 


0.0000 
4 . 3090 
0.2896 
7.1031 


0.0000 

1.8716 

0.2698 

5.7001 


Flight  Condition:  18,000  ft,  0.8  Mach 


0.0000  0.0000  0.0000  1.0000 

■32.1523  -0.0106  33.8536  -30.7268  B 

-0.0014  -0.0001  -1.2760  1.0000 

0.0003  -0.0008  4.5196  -0.6377 


0.0000  0.0000 

-0.7800  9.8209 

-0.1376  -0.2503 
■17.0651  -4.4756 


Flight  Condition:  18,000  ft,  0.75  Mach 


0.0000  0.0000  0.0000  1.0000 

■32.1502  -0.0095  31.0658  -30.1219  B 

-0.0016  -0.0001  -1.1479  1.0000 

0.0003  -0.0009  4.3139  -0.5906 


0.0000  0.0000 

0.6827  8.0728 

-0.1269  -0.2323 
-14.5146  -3.4819 


Flight  Condition:  18,000  ft,  0.7  Mach 


0.0000  0.0000  0.0000  1.0000 

•32.1471  -0.0085  29.1475  -29.9090  B 

-0.0018  -0.0001  -1.0273  1.0000 

0.0003  -0.0010  4.0470  -0.5446 


0.0000  0.0000 
1.8901  -2.0696 
-0.1168  -0. 1970 
-12.2165  -1 .9520 


Flight  Condition:  18,000  ft,  0.65  Mach 


0.0000  0.0000  0.0000  1.0000 

-32.1408  -0.0093  30.5283  -30.8667  B 

-0.0022  -0.0001  -0.9214  1.0000 

0.0003  -0.0010  3.7710  -0.5006 


0.0000  0.0000 
3.0186  -5.9469 
-0.1091  -0.1829 
-10.2335  -1 .0895 


Flight  Condition:  18,000  ft,  0.6  Mach 


0.0000  0.0000  0.0000  1.0000 

•32.1318  -0.0088  31.2822  -32.1163  B 

-0.0026  -0.0001  -0.8267  1.0000 

0.0003  -0.0004  3.4283  -0.4603 


0.0000  0.0000 
3.3522  -4.1756 
-0. 1011  -0. 1787 
-8.5088  -0.6662 


Flight  Condition:  18,000  ft,  0.55  Mach 


0.0000  0.0000  0.0000  1.0000 

-32.1173  -0.0075  31.1300  -34.1163  B 

-0.0033  -0.0001  -0.7600  1.0000 

0.0003  0.0001  2.9861  -0.4264 


0.0000  0.0000 
3.4816  -3.4524 
■0.0945  -0.1602 
-7.1575  -0.4477 


Flight  Condition:  18,000  ft,  0.5  Mach 


0.0000  0.0000  0.0000  1.0000 

A  =  -32.0936  -0.0075  29.8492  -36.9519  B 

-0.0044  -0.0001  -0.6817  1.0000 

0.0004  0.0001  2.7699  -0.3929 


0.0000  0.0000 
-3.4330  -2.4262 
-0.0874  -0.1415 
-5.9225  -0.2430 


Flight  Condition:  18,000  ft,  0.45  Mach 

0.0000  0.0000  0.0000  1.0000 

A  =  -32.0531  -0.0073  24.1220  -40.7828  B 

-0.0059  -0.0002  -0.5797  1.0000 

0.0004  0.0001  2.4959  -0.3597 

Flight  Condition:  18,000  ft,  0.4  Mach 

0.0000  0.0000  0.0000  1.0000 

A  =  -31.9747  -0.0082  18.7640  -46.5509  B 

-0.0086  -0.0002  -0.4866  1.0000 

0.0005  0.0001  1.8995  -0.3399 

Flight  Condition:  18,000  ft,  0.35  Mach 

0.0000  0.0000  0.0000  1.0000 

A  =  -31.8054  -0.0099  19.3550  -55.4223  B 

-0.0133  -0.0002  -0.3987  1.0000 

0.0U06  0.0000  1. 9253d  -0.3394 

Flight  Condition:  18,000  ft,  0.3  Mach 

0.0000  0.0000  0.0000  1.0000 

A  =  -31.3589  0.0030  10.9466  -70.7275  B 

-0.0230  -0.0003  -0.3144  1.0000 

0.0004  0.0001  0.7266  -0.3968 

Flight  Condition:  22,000  ft,  0.9  Mach 

0.0000  0.0000  0.0000  1.0000 

A  =  -32.1556  -0.0089  32.7182  -31.2667  B 

-0.0012  0.0000  -1.3438  1.0000 

0.0003  -0.0004  3.4854  -0.6342 


0.0000  0.0000 

=  2.6617  -1.8464 

-0.0803  -0.1238 
-4.8122  -0.1299 


0.0000  0.0000 

=  1.9813  -3.1842 

-0.0720  -0.0895 
-3.7773  0.0107 


0.0000  0.0000 

=  1.5535  -3.0116 

-0.0632  -0.0564 
-2.9367  0.3274 


0.0000  0.0000 
=  -0.1042  -2.2708 
-0.0552  -0.0518 
-2.1851  0.2266 


0.0000  0.0000 
=  -0.9673  12.3245 
-0.1355  -0.2575 
-19.4559  -6.1317 


Flight  Condition:  22,000  ft,  0.85  Mach 


0.0000  0.0000  0.0000 

A  =  -32.1516  -0.0092  33.1820 

-0.0014  -0.0001  -1.2226 

0.0003  -0.0008  3.6189 

Flight  Condition:  22,000  ft,  0 

0.0000  0.0000  0.0000 

A  =  -32.1465  -0.0088  32.6936 

-0.0016  -0.0001  -1.1084 

0.0003  -0.0008  3.4979 

Flight  Condition:  22,000  ft,  0 

0.0000  0.0000  0.0000 

A  =  -32.1439  -0.0079  29.7815 

-0.0018  -0.0001  -0.9977 

0.0003  -0.0009  3.2901 

Flight  Condition:  22,000  ft,  0 

0.0000  0.0000  0.0000 

A  =  -32.1385  -0.0079  29.2537 

-0.0021  -0.0001  -0.8954 

0.0003  -0.0009  3.1065 

Flight  Condition:  22,000  ft,  0 

0.0000  0.0000  0.0000 

A  =  -32.1295  -0.0085  30.1848 

-0.0025  -0.0001  -0.8049 

0.0003  -0.0009  2.9456 

Flight  Condition:  22,000  ft,  0 

0.0000  0.0000  0.0000 

A  =  -32.1167  -0.0077  30.6728 

-0.0031  -0.0001  -0.7235 

0.0003  -0.0004  2.7733 


1.0000  0.0000  0.0000 

-32.6659  B  =  -0.9029  10.2662 

1.0000  -0.1272  -0.2392 

-0.5912  -16.8790  -4.9160 

8  Mach 

1.0000  0.0000  0.0000 

-34.0209  B  =  -0.7153  8.5536 

1.0000  -0.1188  -0.2213 

-0.5493  -14.5178  -3.8630 

75  Mach 

1.0000  0.0000  0.0000 

-33.3828  B  =  0.7297  3.9226 

1.0000  -0.1091  -0.1963 

-0.5087  -12.3246  -2.7491 

7  Mach 

1.0000  0.0000  0.0000 

-33.8576  B  =  1.9675  -5.0544 

1.0000  -0.1012  -0.1655 

-0.4698  -10.3758  -1.4094 

65  Mach 

1.0000  0.0000  0.0000 

-35.1522  B  =  2.8210  -5.1979 

1 .0000  -0.0944  -0. 1580 

-0.4336  -8.6946  -0.8203 

6  Mach 

1.0000  0.0000  0.0000 

-36.8395  B  =  2.9755  -3.3906 

1.0000  -0.0875  -0.1537 

-0.4003  -7.2322  -0.4626 


Flight  Condition:  22, 

000  f t ,  0.55  Mach 

0.0000 

0.0000 

0.0000  1.0000 

0.0000 

0 . 0000 

A  =  -32.0964 

-0.0064 

29.8824  -39.3114 

B 

=  3.0743 

-2.6019 

-0.0039 

-0.0001 

-0.6564  1.0000 

-0.0817 

-0. 1371 

0.0003 

0.0001 

2.6459  -0.3714 

-6.0831 

-0.2751 

Flight  Condition:  22, 

000  ft,  0.5  Mach 

0 . 0000 

0 . 0000 

0.0000  1.0000 

0.0000 

0.0000 

A  =  -32.0631 

-0 . 0063 

25.0218  -42.7143 

B 

=  2.9967 

-1  .  7863 

-0.0052 

-0.0001 

-0.5651  1.0000 

-0.0756 

-0. 1206 

0.0004 

0.0001 

2.5765  -0.3428 

-5.0334 

-0.1229 

Flight  Condition:  22,000  ft,  0.45  Mach 

0.0000 

0.0000 

0.0000  1.0000 

0.0000 

0.0000 

A  =  -32.0028 

-0.0067 

19.2694  -47.7705 

B 

=  2.3106 

-1.6663 

-0.0072 

-0.0002 

-0.4761  1.0000 

-0.0692 

-0.1055 

0.0004 

0.0001 

2.0689  -0.3229 

-4.0799 

-0. 1306 

Flight  Condition:  22, 

000  ft,  0.4  Mach 

0.0000 

0.0000 

0.0000  1.0000 

0.0000 

0.0000 

A  =  -31.8847 

-0.0080 

16.6361  -55.2176 

B 

=  1.9790 

-4 . 2280 

-0.0105 

-0.0002 

-0.3987  1.0000 

-0.0618 

-0.0573 

0.0005 

0.0000 

1.8429  -0.3112 

-3.2305 

0.2286 

Flight  Condition:  22,000  ft,  0.35  Mach 

0 . 0000 

0.0000 

0.0000  1.0000 

0.0000 

0.0000 

A  =  -31 .6324 

-0.0055 

16.7535  -66.1551 

B 

=  0.6164 

-2.4477 

-0.0163 

-0.0003 

-0.3430  1.0000 

-0.0544 

-0.0507 

0 . 0005 

0.0001 

1.2263  -0.3337 

-2.4932 

0.2732 

Flight  Condition:  22,000  ft,  0.3  Mach 

0.0000 

0.0000 

0.0000  1.0000 

0.0000 

0.0000 

A  =  -30.9268 

0.0149 

9.3812  -86.2052 

B 

=  -0.8816 

-1  .  3931 

-0.0287 

-0.0003 

-0.2812  1.0000 

-0.0489 

-0.0387 

-0.0001 

0.0000 

0.2887  -0.4023 

-1 . 8403 

0. 1853 
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Flight  Condition:  26,000  ft,  0.9  Mach 


0.0000  0.0000  0.0000  1.0000 

A  =  -32.1504  -0.0074  32.3066  -34.8809  B 

-0.0014  0.0000  -1.1604  1.0000 

0.0003  -0.0004  2.7299  -0.5435 

Flight  Condition:  26,000  ft,  0.85  Mach 

0.0000  0.0000  0.0000  1.0000 

A  =  -32.1447  -0.0077  32.6038  -36.7208  B 

-0.0016  -0.0001  -1.0571  1.0000 

0.0003  -0.0008  2.7933  -0.5065 

Flight  Condition:  26,000  ft,  0.8  Mach 

0.0000  0.0000  0.0000  1.0000 

A  =  -32.1395  -0.0074  31.1195  -37.4928  B 

-0.0018  -0.0001  -0.9584  1.0000 

0.0003  -0.0008  2.6541  -0.4705 

Flight  Condition:  26,000  ft,  0.75  Mach 

0.0000  0.0000  0.0000  1.0000 

A  =  -32.1353  -0.0064  28.4955  -37.2294  B 

-0.0021  -0.0001  -0.8636  1.0000 

0.0003  -0.0009  2.4161  -0.4362 

Flight  Condition:  26,000  ft,  0.7  Mach 

0.0000  0.0000  0.0000  1.0000 

A  =  -32.1264  -0.0071  29.1362  -38.5424  B 

-0.0025  -0.0001  -0.7770  1.0000 

0.0003  -0.0009  2.3830  -0.4040 

Flight  Condition:  26,000  ft,  0.65  Mach 


0.0000 
■32  .  1136 
-0.0030 
0.0003 


0.0000  0.0000  1.0000 

-0.0077  29.7819  -40.3133  B 

-0.0001  -0.6998  1.0000 

-0.0009  2.3370  -0.3742 


0.0000  0.0000 
-1.0182  10.5448 
-0.1167  -0.2272 
■16.4857  -5.2653 


0.0000  0.0000 

-1.0436  8.8273 

-0.1095  -0.2104 
-14.2904  -4.2162 


0.0000  0.0000 

-0.3898  6.8315 

-0.1016  -0.1936 
12.2624  -3.2875 


0.0000  0.0000 
0.7821  -0.1848 
-0.0932  -0. 1621 
10.3980  -2.0645 


0.0000  0.0000 
2.0781  -5.9435 
-0.0872  -0.1388 
-8.7604  -0.9883 


0.0000  0.0000 
2.6433  -4.2897 
-0.0813  -0.1354 
-7.3461  -0.6062 


■  »>,y 


Flight  Condition:  26,000  ft,  0.6  Mach 


0.0000 

0.0000 

A  =  -32.0952 

-0.0068 

-0.0037 

-0.0001 

0.0003 

-0.0004 

Flight  Condition:  26, 

0 . 0000 

0.0000 

A  =  -32.0663 

-0.0054 

-0.0047 

-0.0001 

0.0003 

0.0001 

Flight  Condition:  26, 

0.0000 

0.0000 

A  =  -32.0164 

-0.0057 

-0 . 0063 

-0.0001 

0 . 0004 

0 . 0001 

Flight  Condition:  26, 

0.0000 

0.0000 

A  =  -3  1 .9272 

-0.0064 

-0.0087 

-0.0002 

0.0004 

0.0001 

Flight  Condition:  26, 

0.0000 

0.0000 

A  =  -31.7515 

-0.0072 

-0.0128 

-0.0002 

0 . 0005 

0.0001 

Flight  Condition:  26, 

0.0000 

0.0000 

A  =  -31 . 3002 

0.0039 

-0.0210 

-0.0002 

0.0003 

0.0000 

0.0000  1.0000 
29.3384  -42.5014  B 
-0.6217  1.0000 

2.7507  -0.3472 

000  ft,  0.55  Mach 

0.0000  1.0000 
24.8845  -45.5421  B 
-0.5397  1.0000 

2.5580  -0.3227 

000  ft,  0.5  Mach 

0.0000  1.0000 
19.6408  -50.0936  B 
-0.4608  1.0000 

2.1749  -0.3043 

000  ft,  0.45  Mach 

0.0000  1.0000 
16.5774  -56.4384  B 
-0.3935  1.0000 

1.8059  -0.2912 

000  ft,  0.4  Mach 

0.0000  1.0000 
19.8437  -65.6687  B 
-0.3387  1.0000 

1.8617  -0.2993 

000  ft,  0.35  Mach 

0.0000  1.0000 
8.4819  -82.7277  B 
-0.2573  1.0000 

0.6050  -0.3517 


0.0000  0.0000 
2.6237  -2.6088 
-0.0753  -0.1308 
-6.1146  -0.2784 


0.0000  0.0000 
2.6719  -1.6807 
-0.0704  -0.1156 
-5.1425  -0.1073 


0.0000  0.0000 
2.5507  -1.6174 
-0.0651  -0.1021 
-4.2555  -0. 1150 


0.0000  0.0000 
2.0763  -3.6976 
-0.0594  -0.0654 
-3.4448  0.0732 


0.0000  0.0000 
1 . 1662  -2.6180 
-0.0528  -0.0485 
-2.7267  0.2641 


0.0000  0.0000 
-0.1514  -2.5498 
-0.0474  -0.0417 
-2.1251  0.2159 


Flight  Condition:  26,000  ft,  0.3  Mach 


S' 


0.0000  0.0000  0.0000  1.0000 

A  =  -30.2511  0.0360  -2.5161-105.4813 

-0.0361  -0.0003  -0.1702  1.0000 

-0.0005  0.0000  -0.3671  -0.3849 

Flight  Condition:  30,000  ft,  0.9  Mach 


0.0000  0.0000 
B  =  -0.7843  -2.4663 
-0.0449  -0.0252 
1.5641  0.2375 


0.0000  0.0000  0.0000  1.0000  0.0000  0.0000 

A  =  -32.1428  -0.0061  31.8627  -39.4235  B  =  -1.1166  8.9296 

-0.0016  -0.0001  -0.9970  1.0000  -0.0997  -0.1976 

0.0003  -0.0004  2.1557  -0.4624  -13.8388  -4.4916 


Flight  Condition:  30,000  ft,  0.85  Mach 

0.0000  0.0000  0.0000  1.0000 

A  =  -32.1351  -0.0064  31.9249  -41.5880 

-0.0019  -0.0001  -0.9095  1.0000 

0.0003  -0.0007  2.1452  -0.4300 

Flight  Condition:  30,000  ft,  0.8  Mach 

0.0000  0.0000  0.0000  1.0000 

A  =  -32.1298  -0.0057  29.5791  -41.7366 

-0.0021  -0.0001  -0.8244  1.0000 

0.0003  -0.0007  1.9695  -0.4002 

Flight  Condition:  30,000  ft,  0.75  Mach 


0.0000  0.0000 

B  =  -1  . 1622  6 . 8048 

-0.0935  -0.  1813 
-11.9926  -3.5649 


0.0000  0.0000 

B  =  -0.0915  5.2296 

-0.0862  -0.1654 
-10.2681  -2.7791 


0.0000  0.0000  0.0000  1.0000 

A  =  -32.1221  -0.0056  28.2165  -42.3777 

-0.0024  -0.0001  -0.7437  1.0000 

0.0003  -0.0008  1.8191  -0.3719 


0.0000  0.0000 

B  =  1.1123  -3.6678 

-0.0797  -0.1288 
-8.7032  -1  .4175 


Flight  Condition:  30,000  ft,  0.7  Mach 


0.0000  0.0000  0.0000  1.0000 

A  =  -32.1089  -0.0065  29.1151  -44.3135 

-0.0029  -0.0001  -0.6688  1.0000 

0.0003  -0.0008  1.8816  -0.3456 
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0.0000  0.0000 

B  =  2 . 1226  -4 . 9686 

-0.0748  -0.1178 
-7.3431  -0.7534 


Flight  Condition:  30,000  ft,  0.65  Mach 


0.0000  0.0000  0.0000 

A  =  -32.0904  -0.0071  27.4961 

-0.0036  -0.0001  -0.5917 

0.0003  -0.0009  2.1986 

Flight  Condition:  30,000  ft,  0 

0.0000  0.0000  0.0000 

A  =  -32.0640  -0.0060  23.9335 

-0.0045  -0.0001  -0.5044 

0.0003  -0.0004  2.5022 

Flight  Condition:  30,000  ft,  0 

0.0000  0.0000  0.0000 

A  =  -32.0197  -0.0049  19.1986 

-0.0058  -0.0001  -0.4346 

0.0003  0.0001  2.2250 

Flight  Condition:  30,000  ft,  0 

0.0000  0.0000  0.0000 

A  =  -31.9453  -0.0055  16.3750 

-0.0077  -0.0001  -0.3781 

0.0004  0.0001  1.7976 

Flight  Condition:  30,000  ft,  0 

0.0000  0.0000  0.0000 

A  =  -31.8042  -0.0063  18.5182 

-0.0109  -0.0002  -0.3214 

0.0004  0.0000  2.0308 

Flight  Condition:  30,000  ft,  0 

0.0000  0.0000  0.0000 

A  =  -31.5127  -0.0017  12.3981 

-0.0163  -0.0002  -0.2733 

0.0004  0.0000  0.8712 


1.0000  0.0000  0.0000 

-46.6131  B  =  2.4433  -3.3079 

1.0000  -0.0697  -0.1142 

-0.3217  -6. 1665  -0.3962 

6  Mach 

1.0000  0.0000  0.0000 

-49.3741  B  =  2.2557  -1.5574 

1.0000  -0.0646  -0.1090 

-0.3004  -5.1389  -0.0785 

55  Mach 

1.0000  0.0000  0.0000 

-53.6084  B  =  2.2337  -1.5512 

1.0000  -0.0604  -0.0970 

-0.2848  -4.3218  -0.0947 

5  Mach 

1.0000  0.0000  0.0000 

-59.3404  B  =  2.1582  -2.3112 

1.0000  -0.0559  -0.0773 

-0.2731  -3.5762  -0.0735 

45  Mach 

1.0000  0.0000  0.0000 

-67.9412  B  =  1.3656  -2.8570 

1 . 0000  -0.0503  -0.0460 

-0.2759  -2.8834  0.2386 

4  Mach 

1.0000  0.0000  0.0000 

-80.8201  B  =  0.1373  -2.5923 

1 .0000  -0.0448  -0.0431 

-0.3038  -2.2822  0.1731 


Flight  Condition:  30,000  ft,  0 

0.0000  0.0000  0.0000 

A  =  -30.8031  0.0187  8.5235 

-0.0267  -0.0002  -0.2365 

-0.0002  0.0000  0.1110 

Flight  Condition:  34,000  ft,  0 

0.0000  0.0000  0.0000 

A  =  -32.1324  -0.0049  31.3064 

-0.0019  -0.0001  -0.8530 

0.0002  -0.0003  1.6794 

Flight  Condition:  34,000  ft,  0 

0.0000  0.0000  0.0000 

A  =  -32.1226  -0.0049  30.8306 

-0.0022  -0.0001  -0.7789 

0.0002  -0.0006  1.6010 

Flight  Condition:  34,000  ft,  0 

0.0000  0.0000  0.0000 

A  =  -32.1161  -0.0043  28.1732 

-0.0025  -0.0001  -0.7061 

0.0002  -0.0007  1.3805 

Flight  Condition:  34,000  ft,  0 

0.0000  0.0000  0.0000 

A  =  -32.1027  -0.0052  27.7174 

-0.0029  -0.0001  -0.6320 

0.0002  -0.0008  1.5937 

Flight  Condition:  34,000  ft,  0 

0.0000  0.0000  0.0000 

A  =  -32.0835  -0.0060  25.3429 

-0.0035  -0.0001  -0.5556 

0.0003  -0.0008  1.9931 


35  Mach 

1.0000  0.0000  0.0000 

102.0108  B  =  -0.8702  -1.7816 

1.0000  -0.0417  -0.0303 

-0.3536  1.7818  0.1746 

9  Mach 

1.0000  0.0000  0.0000 

-44.7614  B  =  -1.1430  6.7633 

1.0000  -0.0847  -0.1694 

-0.3892  -11.5448  -3.7862 

85  Mach 

1.0000  0.0000  0.0000 

-46.9536  B  =  -0.9896  5.6691 

1.0000  -0.0791  -0.1555 

-0.3626  -9.9930  -3.0202 

8  Mach 

1.0000  0.0000  0.0000 

-46.9351  B  =  0.1797  0.7674 

1.0000  -0.0726  -0.1297 

-0.3383  -8.5408  -2.0317 

75  Mach 

1.0000  0.0000  0.0000 

-48.8071  B  =  1.4527  -5.4077 

1.0000  -0.0681  -0.1018 

-0.3154  -7.2383  -0.9031 

7  Mach 

1.0000  0.0000  0.0000 

-51.3329  B  =  2.1302  -3.6767 

1 .0000  -0.0638  -0.0986 

-0.2944  -6.1176  -0.4988 


Flight  Condition:  34,000  ft,  0.65  Mach 


0.0000  0.0000  0.0000 

A  =  -32.0568  -0.0065  22.2663 

-0.0043  -0.0001  -0.4779 

0.0003  -0.0009  2.2699 

Flight  Condition:  34,000  ft,  0 

0.0000  0.0000  0.0000 

A  =  -32.0148  -0.0057  17.9322 

-0.0054  -0.0001  -0.4001 

0.0003  -0.0004  2.2327 

Flight  Condition:  34,000  ft,  0 

0.0000  0.0000  0.0000 

A  =  -31.9478  -0.0047  15.7087 

-0.0071  -0.0001  -0.3535 

0.0003  0.0001  1.8135 

Flight  Condition:  34,000  ft,  0 

0.0000  0.0000  0.0000 

A  =  -31.8282  -0.0056  16.7474 

-0.0096  -0.0001  -0.3003 

0.0004  0.0000  2.0240 


1.0000  0.0000  0.0000 

-54.2504  B  =  2.2139  -2.3546 

1.0000  -0.0595  -0.0951 

-0.2756  -5.1454  -0.2089 

6  Mach 

1.0000  0.0000  0.0000 

-58.3606  B  =  1.8413  -1.4707 

1.0000  -0.0553  -0.0907 

-0.2647  -4.2934  -0.0763 

55  Mach 

1.0000  0.0000  0.0000 

-63.7817  B  =  1.8574  -1.3914 

1.0000  -0.0517  -0.0810 

-0.2546  -3.6097  -0.1241 

5  Mach 

1.0000  0.0000  0.0000 

-71.7205  B  =  1.3982  -3.5083 

1.0000  -0.0472  -0.0436 

-0.2551  -2.9656  0.1677 
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Flight  Condition:  34,000  ft,  0.45  Mach 


0.0000  0.0000  0.0000  1.0000 

A  =  -31.6185  -0.0046  14.3289  -81.8577  B 

-0.0135  -0.0002  -0.2700  1.0000 

0.0004  0.0000  1.0496  -0.2694 

Flight  Condition:  34,000  ft,  0.4  Mach 

0.0000  0.0000  0.0000  1.0000 

A  =  -31.0439  0.0105  8.0406-103.9423  B 

-0.0216  -0.0002  -0.2205  1.0000 

0.0000  0.0000  0.3937  -0.3236 


0.0000  0.0000 
0.5153  -2.5631 
-0.0424  -0.0408 
-2.3903  0.1522 


0.0000  0.0000 
-0.7184  -2.4958 
-0.0393  -0.0297 
-1.9481  0.2287 


j  i 


Flight  Condition:  38,000  ft,  0.9  Mach 


0.0000  0.0000  0.0000  1.0000 

A  =  -32.1172  -0.0038  30.8408  -51.7604  B  = 

-0.0022  -0.0001  -0.7213  1.0000 

0.0002  -0.0003  1.2864  -0.3225 

Flight  Condition:  38,000  ft,  0.85  Mach 

0.0000  0.0000  0.0000  1.0000 

A  =  -32.1070  -0.0037  28.9780  -53.1253  B  = 

-0.0025  -0.0001  -0.6545  1.0000 

0.0002  -0.0006  1.2741  -0.3015 

Flight  Condition:  38,000  ft,  0.8  Mach 

0.0000  0.0000  0.0000  1.0000 

A  =  -32.0951  -0.0037  26.2610  -54.2606  B  = 

-0.0029  -0.0001  -0.5870  1.0000 

0.0002  -0.0007  1.4028  -0.2820 

Flight  Condition:  38,000  ft,  0.75  Mach 

0.0000  0.0000  0.0000  1.0000 

A  =  -32.0744  -0.0048  23.0363  -57.1571  B  = 

-0.0035  -0.0001  -0.5114  1.0000 

0.002  -0.0007  1.8436  -0.2638 

Flight  Condition:  38,000  ft,  0.7  Mach 

0.0000  0.0000  0.0000  1.0000 

A  =  -32.0459  -0.0057  20.1578  -60.4982  B  = 

-0.0042  -0.0001  -0.4428  1.0000 

0.0002  -0.0008  2.0441  -0.2484 

Flight  Condition:  38,0000  ft,  0.65  Mach 

0.0000  0.0000  0.0000  1.0000 

A  =  -32.0028  -0.0065  16.3861  -64.9370  B  = 

-0.0053  -0.0001  -0.3751  1.0000 

0.0003  -0.0009  1.9789  -0.2394 


0.0000  0.0000 

-1.1699  5.5244 

-0.0709  -0.1430 
-9.5743  -3.1572 


0.0000  0.0000 

-0.4986  4.7349 

-0.0657  -0.1320 
-8.2649  -2.5621 


0.0000  0.0000 
0.6600  -3.3294 
-0.0608  -0.0950 
-7.0516  -1.2971 


0.0000  0.0000 
1.6924  -3.8681 
-0.0574  -0 . 0835 
-5.9839  -0.5967 


0.0000  0.0000 
2.0538  -3.0343 
-0.0538  -0.0816 
-5.0675  -0.3653 


0.0000  0.0000 
1.9161  -2. 1631 
-0.0503  -0.0788 
-4.2700  -0.2208 


Flight  Condition:  38,000  ft,  0 

0.0000  0.0000  0.0000 

A  =  -31.9389  -0.0057  14.5232 

-0.0067  -0.0001  -0.3194 

0.0003  -0.0005  1.8149 

Flight  Condition:  38,000  ft,  0 

0.0000  0.0000  0.0000 

A  =  -31.8311  -0.0049  15.5895 

-0.0088  -0.0001  -0.2740 

0.0003  0.0000  2.0114 

Flight  Condition:  38,000  ft,  0 

0.0000  0.0000  0.0000 

A  =  -31.6553  -0.0048  14.8821 

-0.0119  -0.0001  -0.2523 

0.0003  0.0000  1.2438 

Flight  Condition:  38,000  ft,  0 

0.0000  0.0000  0.0000 

A  =  -31 . 1779  0.0064  6. 1007 

-0.0182  -0.0002  -0.1880 
0.0001  0.0000  0.4764 


6  Mach 

1.0000  0.0000  0.0000 

-70.2645  B  =  1 .4981  -1.3132 

1.0000  -0.0467  -0.0747 

-0.2341  -3.5656  -0.1188 

55  Mach 

1.0000  0.0000  0.0000 

-77.8026  B  =  1.2700  -3.6987 

1.0000  -0.0433  -0.0398 

-0.2346  -2.9784  0.1275 

5  Mach 

1.0000  0.0000  0.0000 

-87.C302  B  =  0.5640  -2.5569 

1.0000  -0.0391  -0.0376 

-0.2439  -2.4339  0.1274 

45  Mach 

1.0000  0.0000  0.0000 

108.6833  B  =  -0.4276  -2.9478 

1.0000  -0.0365  -0.0273 

-0.2895  -2.0557  0.2277 


Appendix  B 


This  appendix  presents  the  mathematical  development  of 
obtaining  the  Nth  order  autoregressive  vector  difference 
equation  that  results  in  a  set  of  parameters  representing  an 
aircraft  model  at  a  given  flight  condition.  This  appendix 
was  taken  largely  from  material  presented  by  Pineiro  (1). 
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The  difference  equation  model  used  in  this  thesis  for 
the  open-loop  longitudinal  dynamics  of  the  host  airplane  is 
based  on  the  discrete  state  and  output  relationships  given 
in  Equations  (2-8)  and  (2-9)  (Equations  (2-8)  and  (2-9) 
repeated  here  for  reading  continuity): 

x{ (k  +  1 ) T }  =  $  x(kT)  +  $  u(kT)  (A-l ) 

y(kT)  =  C  x(kT)  (A-2) 

By  taking  the  Z  transform  of  Equations  (A-l)  and  (A-2)  and 
obtaining  a  transfer  function  model,  the  desired  input- 
output  relationship  can  be  obtained. 

*j{x(k  +  1))  =  zl£(z)  =$X(z)  +  fu(z)  ( A-  3 ) 

^(yfk)}  =  y(z)  =  CX  ( z )  ( A-4  ) 

Rearranging  equation  (A-3)  yields 

zIX(z)  -  $X(z)  =  ¥U(z)  ( A- 5 ) 

Therefore , 

[zl  -  $]X(z)  =  ¥y(z)  ( A-6  ) 

X ( z )  =  [zl  -  frtufz)  ( A- 7 ) 

By  substituting  Equation  (A-2)  into  Equation  (A-7)  and  per- 
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forming  some  matrix  manipulation  the  following  expression 


for  Y(z)  is  obtained 


Y  ( z )  =  C[zl  -  <1>  ]  *U(z) 

=  G  (  z  )  U  (  z  ) 


( A-8  ) 
(  A — 9  ) 


The  m  x  1  vector  Y(z),  the  m  x  m  matrix  G{z),  and  the  m  x  1 
vector  U(z)  are  given  by 


Y(  z) 


Y  (z) 
Y2(z) 


Y3(z) 


( A- 1 0 ) 


G(z) 


u'2' 

2 1  Z 

C12U‘ 

... 

...  G  ( z) 
lm 

...  G0  (z) 
2m 

•  •  •  •  •  • 

»1(21 

G  „lz) 
m2 

...  G  (z) 
mm 

(  A-  1  1  ) 


U(z)  = 


U  (z) 
U2U) 


U  (z) 
m 


The  elements  of  G(z),  denoted  G.  (z),  are  the  transfer 

i  J 

functions  relating  the  output  Y^  to  the  control  input  V 


(A- 12 ) 


and  are  of  the  form 


v;»»rv 


rn  *'«  it 


G.  .(  z) 
i  J 


(w  <  n)  ( A- 1 3 ) 


w  w  —  1 

bz  +  b„z  +  ..  +  b  z  +  b 
12  w  w+1 


zn  +  a.zn-1  +  ..  +  a  ,z  +  a 
1  n- 1  n 


By  dividing  the  numerator  and  denominator  of  each  transfer 
function  of  G(z)  by  z™ ,  the  transfer  function  matrix  is 
transformed  into  the  delay  operator  form  given  by 


,  w-n  ,  w-n-1  ,  -n 

bz  +  b  z  +..+b  ,  z 

1  2  w+ 1 

G  .  .  (  z  )  =  -  (  w  <  n  )  (  A  -  1  4  ) 

i  J 

1  +  a z-1  +  ..  +  a  z-n+1  +  a  z-n 
1  n- 1  n 

By  grouping  the  coefficients  with  the  same  amount  of 

delav  in  every  transfer  function  G  (z)  in  Glz),  and  taking 

i  j 

the  inverse  Z  transform,  the  Nth  order  autoregressive  dif¬ 
ference  equation  of  the  open-loop  plant  is  obtained 


Z(kT)  =  B  u{(k-l)T)  -  A  2 { ( k- 1 ) T }  + 


+  B  u { <  k - N ) T  3  -  A  v { ( k-N ) T }  +  e(kT)  ( A- 1 5  > 

n  ~ n 


y(kT)  =  T  (kT)  6  +  e(kT) 


( A-  1 6  ) 


where  e(kT)  is  an  equation  error  term  assumed  to  be  a  zero- 

mean  Gaussian  white  noise  vector  added  to  account  for  roodel- 
■v'T  mxl 

ling  errors,  !  (kT)€  R  is  a  matrix  of  past  values  of 


mxm 

{£<kT)  and  u(kT)) ,  and  the  matrices  A^€  R  ( i = 1 , 2 , . . . , N ) , 

mxm  n  1 

B  £  R  (i=l,2,...N)  and  the  vector  v  6  R  are  the  parame- 
i 

ters  of  the  Nth  order  difference  equation. 


f  J 
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Appendix  C 
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This  Appendix  presents  the  position  limits  for  the 
elevators  and  flaperons  for  each  flight  condition.  As  dis¬ 
cussed  in  Chapter  2,  the  reason  for  the  different  position 
limits  is  due  to  the  fact  that  the  trim  values  for  the 
elevator  and  flaperon  positions  differ  for  each  flight  con- 
d i t i on . 


Flight  Condition  Elevator  Position  Flaperon  Position 
Alt  (Kft)  Mach  Limits  Limits 


10 

.  90 

27.3171, 

-22.6829 

22.0000, 

-21.0000 

10 

•  85 

27.1407, 

-22.8259 

22.0000, 

-21.0000 

10 

.80 

27.0100, 

-22.9900 

22.0000, 

-21.0000 

10 

.75 

26.9343, 

-23.0657 

21.5916, 

-21.4084 

10 

.  70 

26.9769, 

-23.0231 

20.5543, 

-22 .4457 

10 

.65 

26.9687  , 

-23.0313 

20.5019, 

-22.4981 

10 

.60 

26.9011  , 

-23.0988 

18 . 4067  , 

-24 . 5933 

10 

.  55 

27.0884, 

-22.9116 

17.2362, 

-25.7638 

10 

.50 

27.2365, 

-22.7635 

15.9609  , 

-27.0390 

10 

.45 

27 .3375, 

-22.6625 

14.4776 , 

-28.5224 

10 

.40 

27 . 3085, 

-22.6915 

12.6417, 

-30.3583 

10 

.  35 

27.1404 , 

-22 . 8597 

10.1102, 

-32 . 8898 

10 

.  30 

26 . 1426 , 

-23 .8574 

6.0729, 

-36.9271 

10 

.25 

25 . 1583 , 

-24.8417 

5.0000, 

-38.0000 

14 

.  90 

27.2899, 

-22.7100 

22.0000, 

-21.0000 

14 

.85 

27 .1375, 

-22.8625 

22.0000, 

-21 .0000 

14 

.80 

26.9598, 

-23.0402 

22.0000, 

-21 .0000 

14 

.75 

26.9439  , 

-23.0561 

21 . 3222 , 

-21.6778 

14 

.70 

26.9699  , 

-23.0301 

20.2337, 

-22.7664 

14 

.65 

26.9715, 

-23.0285 

19.0859, 

-23.0141 

14 

.60 

26.8681  , 

-23.1319 

17.9156, 

-25.0844 

14 

.55 

27.0560, 

-22.9440 

16.6666, 

-26.3333 

14 

.  50 

27.2062, 

-22.7939 

15.2683, 

-27.7317 

14 

.45 

27 . 2496, 

-22 .7504 

13.6179, 

-29.6321 

14 

.40 

27 . 1438, 

-22.8562 

1 1 .5134  , 

-31 . 4867 

14 

.  35 

26 . 7824  , 

-23 .2176 

8.4601  , 

-34 . 5399 

14 

.  30 

25.3368, 

-24 . 6633 

5.0000, 

-38.0000 

14 

.25 

25.2301, 

-24 . 7699 

5.0000, 

-38.0000 

18 

.90 

27.2598, 

-22 .7402 

22.0000 , 

-21.0000 

18 

.85 

27.0957, 

-22.9043 

22.0000, 

-21.0000 

18 

.80 

26.9013, 

-23 . 0987 

22.0000, 

-21.0000 

18 

.  75 

26 .9586  , 

-23.0414 

20.9971 , 

-22.0029 

18 

.  70 

26.9793, 

-23.0207 

19.8341 , 

-23.1659 

18 

.65 

26 . 9669  , 

-23.0331 

18 . 5947 , 

-24.4053 

18 

.60 

26 .8530, 

-23 . 1470 

17 . 3436 , 

-25.6564 

18 

.55 

27.0466, 

-22.9534 

15.9830, 

-27.0170 

18 

.50 

27 .1503, 

-22.8497 

14.4367, 

-28.5633 

18 

.45 

27 . 1075  , 

-22.6925 

12.5850, 

-30.4150 

18 

.  40 

26.9715, 

-23.0285 

10.0985, 

-32.9015 

18 

.  35 

26 . 1469 , 

-23 . 8531 

6.4130, 

-36 . 5870 

18 

.  30 

25.0717, 

-24.9282 

5.0000, 

-38.0000 

22 

.90 

27 . 2262  , 

-22 .7738 

22.0000, 

-21.0000 

22 

.85 

27.0477, 

-22.9532 

22.0000, 

-21.0000 

22 

.80 

26.8873, 

-23. 1127 

21.8248, 

-21 . 1753 

22 

.75 

26 . 9801 , 

-23.0199 

20.6028, 

-22.3972 

22 

.70 

27.0164 , 

-22.9836 

19.3353, 

-23 . 6647 

22 

.65 

27.0194, 

-22 . 9806 

18.0135, 

-24.9865 

22 

.60 

26 . 8667  , 

-23 . 1333 

16.6537, 

-26 . 3463 

22 

.  55 

27.0119, 

-22 . 9880 

15. 1578, 

-27.8422 

22 

.  50 

27.0199, 

-22.9801 

13.5678, 

-29 . 5678 

22 

.45 

26.9190, 

-23.0810 

1 1 . 2427 , 

-31.7573 

22 

.  40 

26 . 5847  , 

-23.4153 

8.2382, 

-33 . 7618 

22 

.  35 

25 . 3709  , 

-24.6792 

5.0000, 

-38.0000 

22 

.  30 

25 . 1478  , 

-24.8523 

5.0000, 

-38.0000 

26 

.90 

27 . 1877 , 

-22 .8123 

22.0000, 

-21 .0000 

26 

.  85 

26.9906, 

-23.0094 

22.0000, 

-21 .0000 

26 

.80 

26.9366, 

-23.0634 

21.4198, 

-21 . 5802 

26 

.75 

27.0291, 

-22.9709 

20.1317, 

-22.8683 

26 

.  70 

27 . 1280 , 

-22 . 8720 

18.7323, 

-24 . 2677 

26 

.65 

27.0930, 

-22.9170 

17.3037, 

-25.6963 

26 

.60 

26.8501  , 

-23.1499 

15.8139, 

-27.1861 

26 

.  55 

26.8977, 

-23 . 1023 

14.1526, 

-28.8474 

26 

.  50 

26.8249, 

-23.1751 

12.1327 , 

-30.8673 

26 

.45 

26.7082, 

-23.2918 

9.5537  , 

-33.4464 

26 

.  40 

25.7882, 

-24.2118 

5.9549, 

-37.0451 

26 

.  35 

25.0313, 

-24.9687 

5.0000, 

-38.0000 

26 

.  30 

25.1150, 

-24.8850 

5.0000, 

-38.0000 

30 

.  90 

27. 1405, 

-22.8595 

22.0000, 

-21.0000 

30 

.85 

26 . 9340, 

-23.0660 

22.0000, 

-21 .0000 

30 

.80 

27.0269, 

-22.9731 

20.9211 , 

-22.0789 

30 

.  75 

27 . 1622, 

-22.8378 

19.5016, 

-23 .4984 

30 

.  70 

27 . 2342  , 

-22 . 7658 

17.9817, 

-25.0183 

30 

.65 

27.1207  , 

-22.8793 

16 . 4251  , 

-26.5749 

30 

.60 

26 .7404  , 

-23 . 2596 

14 . 7789  , 

-28.2211 

30 

.55 

26.6819, 

-23.3181 

12.8383, 

-30.1617 

30 

.  50 

26 .6237  , 

-23 . 3767 

10.4812, 

-32.5189 

30 

.45 

26.0799, 

-23.9201 

7.2889, 

-35.7111 

30 

.40 

25 . 1127  , 

-24.8873 

5.0000, 

-38.0000 

30 

.  35 

24.9556, 

-25.0444 

5.0000, 

-38.0000 

34 

.90 

27.0970, 

-22.9030 

22.0000, 

-21.0000 

34 

.85 

26.9523, 

-23.0477 

21 .7643, 

-21 .2357 

34 

.80 

27.1475, 

-22.8525 

20.3042, 

-22.6952 

34 

.75 

27 . 3222  , 

-22.6778 

18.7056, 

-24.2944 

34 

.70 

27 . 3322  , 

-22.6678 

17.0547, 

-25.9453 
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This  Appendix  presents  the  parameter  vectors  for  the 
flight  conditions  specifically  addressed  in  this  thesis. 
Recall  that  the  parameter  vector  is  a  set  of  parameters 
associated  with  a  given  flight  condition  that  are  used  in 
the  autoregressive  difference  equation  form  of  the  aircraft 
model  . 
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Flight  Condition:  10,000  ft,  0.35  Mach 


£ 


* 


% 


A1  = 

-3 . 990568 

ai 

= 

0.000842 

-0.039268 

0.001033 

0.000310 

A2  = 

5.971522 

“2 

= 

-0.002524 

0. 117563 

-0 . 003096 

-0.000950 

A3  = 

-3.971342 

B3 

= 

0.002520 

-0.117320 

0.003092 

0.000969 

A4  -- 

0.990387 

"4 

- 

-0.000838 

0.039025 

-0.001029 

-0.000329 

Flight 

Condition : 

10,000 

ft,  0 

.40 

Mach 

A.  = 

-3 . 989059 

Ei 

= 

0.000957 

-0.051883 

0.001442 

-0.001841 

A2  = 

5.966940 

“2 

= 

-0.002871 

0.155263 

-0.004321 

0.005471 

A3  = 

-3 . 966701 

-3 

= 

0.002866 

-0.154875 

0.004315 

-0.005419 

II 

< 

0.988821 

-4 

- 

-0.000952 

0.051495 

-0.001436 

0.001789 

Flight 

Condition : 

10,000 

ft,  0 

.50 

Mach 

A1  = 

-3.986201 

ai 

= 

0.001149 

-0.080647 

0.001891 

-0.005754 

A2  = 

5.958298 

“2 

= 

-0.003451 

0.241174 

-0.005664 

0.017146 

A3  = 

-3.957994 

63 

= 

0.003444 

-0.240404 

0.005655 

-0.017027 

It 

< 

0.985896 

-4 

- 

-0.001143 

0.079878 

-0.001881 

0.005636 

252 
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Flight 

Condition : 

10,000 

ft,  0 

.65 

Mach 

A1  = 

-3.981994 

= 

0.001427 

-0. 139296 

0.002386 

-0.016365 

A2  = 

5 . 945443 

®2 

= 

-0.004298 

0.416132 

-0.007146 

0.048754 

A  3  = 

-3.944927 

= 

0.004279 

-0.414375 

0.007130 

-0.048414 

> 

II 

0.981466 

= 

-0.001417 

0.137539 

-0.002370 

0.016024 

Flight 

Condition : 

10,000 

ft,  0. 

75 

Mach 

h  = 

-3.977787 

“l 

= 

0.001672 

-0. 197434 

0.002931 

-0.046029 

II 

CM 

< 

5.932799 

~2 

= 

-0.005032 

0.589214 

-0.008782 

0.137196 

A3  = 

-3.932239 

S3 

= 

0.005019 

-0.586125 

0.008757 

-0. 136303 

A4  = 

0 .977226 

= 

-0.001658 

0. 194345 

-0.002908 

0.045136 

F 1 ight 

Condition : 

10,000 

ft,  0 

.80 

Mach 

A1  = 

-3.975528 

B1 

= 

0.001793 

-0.231040 

0.003139 

-0.059402 

A2  = 

5.926021 

S2 

= 

-0.005404 

0.689126 

-0.009403 

0.176989 

A3  = 

-3.925456 

= 

0.005388 

-0.685132 

0.009376 

-0.175772 

A 


4 


B 


0.974963 


-0.001778 

0.227046 


-0.003112 

0.058185 


Flight  Condition: 


10,000  ft,  0.90  Mach 


A1 = 

-3.970230 

B1  = 

0.002019 

-0.307304 

0.003569 

-0.093715 

A2  = 

5.911380 

&2  = 

-0.006103 

0.915522 

-0.010695 

0.278993 

II 

CO 

< 

-3.910892 

"3  = 

0.006083 

-0.909132 

0.010660 

-0.276841 

A,  = 

0.970134 

B.  = 

-0.002000 

-0.003534 

4 

4 

0.300914 

0.091562 

Flight 

Condition : 

14,000 

ft,  0.80 

Mach 

A1  = 

-3.978575 

Bi  = 

0 . 001567 

-0. 198759 

0.002800 

-0.051639 

5 .935265 

= 

-0.004720 

-0.008388 

2 

L 

0.593277 

0. 154001 

A,  s 

-3.934805 

B,  = 

0.004709 

0.008367 

3 

O 

-0.590276 

-0.153085 

II 

< 

0.978115 

“4  = 

-0.001555 

0.195758 

-0.002779 

0.050722 

Flight 

Condi tion : 

18,000 

ft,  0.45 

Mach 

A1  = 

-3 . 990800 

B1  = 

0.000800 

-0.048047 

0.001234 

-0.001312 

A2  = 

5.972188 

~2  = 

-0.002400 

0.143841 

-0.003699 

0.003898 

A3  = 

-3.971961 

~3  = 

0.002397 

-0. 143541 

0.003694 

-0.003859 

A.  = 

0.905786 

B,  = 

-0.000797 

-0.001229 

4 

4 

0.047746 

0.001274 

Flight  Condition:  18,000  ft,  0.80  Mach 


A  = 

-3.981308 

B 

=  0.001363 

0.002486 

1 

~1 

-0.170151 

-0.044673 

A2  = 

5.943556 

“2 

=  -0.004105 

0.508218 

-0.007449 

0. 133337 

A  = 

-3.943189 

B 

=  0.004096 

0.007433 

3 

“3 

-0.505982 

-0.132653 

A  = 

0.980940 

B 

=  -0.001355 

-0.002470 

4 

“4 

0.167915 

0.043990 

Flight 

Condition : 

22,000  ft,  0 

.45  Mach 

A1  = 

-3.992165 

Ei 

=  0.000690 

-0.040741 

0.001052 

-0.001314 

A2  = 

5.976305 

=  -0.002071 

0.122015 

-0.003153 

0.003916 

A3  = 

-3.976116 

=  0.002068 

-0.121805 

0.003149 

-0.003888 

A  = 

0.991975 

=  -0.000688 

-0.001048 

4 

_4 

0.040531 

0.001286 

Flight  Condition:  22,000  ft,  0.65  Mach 


A1  = 

-3.987867 

Ei  = 

0.000938 

-0.086776 

0.001573 

-0.008208 

II 

CM 

< 

5.963342 

~2 

-0.002819 

0.259598 

-0.004714 

0.024513 

A3  = 

-3.963082 

h  -- 

0.002815 

-0.258867 

0.004706 

-0.024401 

> 

It 

0.987607 

h  -- 

-0.000934 

0.086045 

-0.001566 

0.008096 

255 

Flight 

Condition : 

22,000 

ft,  0.90 

Mach 

A1  = 

-3.992314 

B1  = 

0.001351 

-0. 194202 

A2  = 

5.976828 

B2  = 

-0.004055 

0.581808 

A3  = 

-3.976713 

~3  = 

0.004050 

-0.581010 

> 

II 

0.992199 

“4  - 

-0.001345 

0 . 193404 

Flight 

Condi tion : 

26 , 000 

ft,  0.70 

Mach 

A1  = 

-3.988397 

E1  = 

0.000867 

-0.087440 

A2  = 

5.964985 

b2  * 

-0.002607 

0.261618 

A  3  = 

-3.964777 

“3  = 

0.002603 

-0.260914 

A4  = 

0.988189 

S,  = 

-0.000864 

0.086737 

Flight 

Condition : 

30 , 000 

ft,  0.50 

Mach 

A1  = 

-3 . 993622 

0.000557 

-0.035718 

A2  = 

5 .980699 

= 

-0.001673 

0.107009 

> 

00 

II 

-3.980531 

0.001671 

-0.106863 

II 

rr 

< 

0.993454 

~4  = 

-0.000556 

0.035572 

0.002570 

-0.061218 

-0.007702 
0. 183375 

0.007692 
-0. 183095 

-0.002559 

0.060938 


0.001382 

-0.009880 

-0.004142 

0.029530 

0.004137 

-0.029421 

-0.001376 

0.009770 


0.000771 

-0.000740 

-0.002313 

0.002205 

0.002311 

-0.002189 

-0.000769 

0.000724 
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Flight 

Cond i t i on : 

34,000 

ft,  0 

.60 

Mach 

A1  = 

-3.988733 

*1 

= 

0 . 000786 

-0 . 099757 

0.001548 

-0.030160 

A2  = 

5.966067 

~2 

= 

-0.002364 

0.298482 

-0.004641 

0.090220 

A3  = 

-3.965934 

“3 

= 

0.002362 

-0.297690 

0.004636 

-0.089960 

A4  = 

0.988600 

~4 

= 

-0.000783 

0.098966 

-0.001542 

0.029899 

Flight 

Condition : 

34,000 

ft,  0 

.85 

Mach 

A1  = 

-3.993531 

*1 

= 

0 . 000551 

-0.042884 

0.000905 

-0.000771 

II 

CS) 

< 

5.980380 

• 

“2 

= 

-0.001654 

0.128468 

-0.002714 

0.002292 

A3  = 

-3.980166 

~3 

= 

0.001652 

-0.128282 

0.002711 

-0.002268 

A4  = 

0.993317 

*4 

= 

-0.000549 

0.042698 

-0.000903 

0.000748 

Flight 

Condition : 

38,000 

ft,  0 

.60 

Mach 

A1  = 

-3.994599 

&1 

- 

0.000466 

-0.035619 

0.000745 

-0.001193 

A2  = 

5 . 983623 

~2 

= 

-0.001398 

0.106733 

-0.002236 

0.003562 

II 

CO 

< 

-3.983447 

~3 

= 

0.001397 

-0.106609 

0.002234 

-0.003545 

> 

II 

0.994423 

~4 

= 

-0.000465 

0.035495 

-0.000744 

0.001176 

Flight  Condition:  38,000  ft,  0.65  Mach 


A1 

= 

-3.994002 

0.000501 

0.000786 

-0.002213 

-0.042655 

A„ 

_ 

5.981814 

= 

-0.001504 

-0.002357 

2 

“2 

0. 127793 

0.006616 

A  „ 

_ 

-3.981620 

Bn  = 

0.001503 

0.002355 

3 

“3 

-0. 127620 

-0.006592 

A , 

_ 

0.993808 

B,  = 

-0.000500 

-0.000784 

4 

~4 

0.042482 

0.002189 

Flight 

Cond i t i on : 

38,000 

ft,  0.70 

Mach 

A1 

-3.993253 

*1  = 

0 . 000536 

0.000814 

-0.050619 

-0.003657 

A  . 

_ 

5.979562 

B„  = 

-0.001611 

-0.002441 

2 

“2 

0. 151622 

0.010938 

A 

_ 

-3.979365 

Bn  ' 

0.001609 

0.002439 

3 

~3 

-0. 151384 

-0.010905 

A. 

_ 

0 .993056 

B.  = 

-0.000535 

-0.000812 

4 

-4 

0.050382 

0.003624 

Flight 

Condition : 

38 , 000 

ft,  0.80 

Mach 

A 

_ 

-3.991436 

B,  = 

0.000605 

0.000946 

1 

-1 

-0.070422 

-0.012959 

A 

_ 

5.974184 

Bn  = 

-0.001818 

-0.002837 

2 

-2 

0.210843 

0.038788 

A0 

_ 

-3.974057 

Bn  = 

0.001816 

0.002835 

3 

-3 

-0.210419 

-0.038698 

A. 

0.991310 

B.  = 

-0.000603 

-0.000943 

4 

"4 

0.069998 

0.012869 
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Flight  Condition:  38,000  ft,  0.90  Mach 


A  =  -3.989684 


A  =  5.968947 

2 


A  =  -3.968841 
3 


A  =  0.989578 

4 


B  =  0.000705 

1  -0.095594 


~2  = 


"3 


-0.002120 

0.286084 

0 . 002118 
-0 . 285385 

-0 . 000703 
0.094895 


0.001424 

■0.031531 

■0.004271 

0.094348 

0.004267 

-0.094101 

-0.001420 

0.031285 
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